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Quarter-Wave Corrugated Transformer for 
Broadband Matching of a Corrugated Feed 


By C. DRAGONE® and W. E. LEGG* 
(Manuscript received June 3, 1983) 


This paper describes a broadband technique for obtaining the transforma- 
tion TE,, — HE; at the input of a corrugated feed. Usually, this transfor- 
mation is obtained by operating the feed in the vicinity of a frequency, w, 
corresponding to one of the zeroes of the input surface reactance since the 
feed input reflection, p,, vanishes at w;. Here we use a matching transformer 
to produce p, ~ 0 at two frequencies simultaneously, i.e., at w; and a much 
lower frequency, wo. Using two sections of corrugated waveguides, with a total 
of eight teeth and wy = 0.63 w;, the measured reflection, p;, is about —40 dB 
in the vicinity of both frequencies wo and w;, and it remains less than —30 dB 
over the entire frequency range 0.60 w:; < w < 1.04 w,. This makes the 
corrugated feed suitable for simultaneous operation at the 4- and 6-GHz 
terrestrial microwave radio bands. 


I. INTRODUCTION 


Corrugated feeds are usually characterized by a relatively large input 
reflection,’ * which vanishes only at certain frequencies corresponding 
to the zeroes of the surface reactance due to the input corrugations. 
In this article we describe a simple technique for reducing considerably 
this reflection over a wide range of frequencies. The reflection in 
question arises at the input of the feed as shown in Fig. 1, where a 
corrugated waveguide is directly connected to an uncorrugated wave- 
guide of circular cross sections. The reflection, p;, due to the discon- 
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Fig. 1—Junction between two circular waveguides, one with corrugations of depth, d. 
The reflection, p;, is given accurately by eq. (1). 





tinuity in surface reactance between the two waveguides, is given 
accurately by the formula” 


_ bo-8 
Bo + B’ 


relating p;.to the propagation constants (8 and 6) of the dominant 
modes (TE; and HE;;) of the two waveguides. The two propagation 
constants coincide only at. certain isolated frequencies, corresponding 
to the zeroes of the surface reactance of the input corrugations. The 
feed is normally operated in the vicinity of the first zero, which is the 
frequency, w;, determined by the condition 


A, (w1) 
9 ? 
where d is the depth of the corrugations and },(w) is the wavelength 
for the radial waves excited in the grooves at an input frequency, w. 
Using eq. (1), one can readily determine, as in Refs. 2 and 3, the 
variation of p,; with frequency, in the vicinity of @. 
One finds a rapid increase in p;, with |w — w,| causing, under the 





Pt (1) 


c= (2) 
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Fig. 2—The two waveguides of Fig. 1 are matched by using N sections of waveguides 
with corrugations of the same depth, d, but different s;/h;. 


conditions of Ref. 2,* reflections greater than —30 dB for w < 0.832 
w;. Thus, the junction of Fig. 1 is unsuitable for applications requiring 
negligible | p,|? at widely spaced frequencies, such as, for instance, 4 
and 6 GHz, as required in terrestrial radio systems.’ 

In this article we place a matching transformer between the two 
waveguides of Fig. 1 to cause p; = 0 also in the vicinity of a frequency 
wo appreciably lower than w;. The transformer, shown in Fig. 2, 
consists of two sections of corrugated waveguide having the same d 
but different gaps s; between successive teeth in each section. Since 
all corrugations have the same depth, d, the total reflection, p;, 
vanishes for w = w;. To obtain p; = 0 in the vicinity of wo, the lengths, 
4, of the various sections and the values of s; are chosen by the same 
procedure commonly used for multisection quarter-wave matching 
transformers.* By properly choosing the parameters 4 and s;, one can 
cause p; to have N zeroes in the vicinity of wo. Here we make these 
zeroes coincide with wo to obtain a maximally flat characteristic in the 


* Notice a misprint in Ref. 2. The frequency w, = 17.5 GHz, given correctly on page 
878, is given incorrectly as 19 GHz on pp. 810 and 876. 

*The TD and TH bands correspond to the intervals (3.7, 4.2) and (5.925, 6.325) in 
GHz. 
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vicinity of w = wo. As a result, the reflection, p;, is very small at 
frequencies close to wo and w;. Furthermore, we shall see that p; is also 
small in the entire interval (w9, w). 

Section III describes a transformer with inside radius a with a 
measured reflection of less than —30 dB for values of ka in the interval 
(2.437, 4.22), and less than —35 dB in (2.437, 2.766) and (3.69, 4.22). 
This performance is obtained using two sections with ordinary teeth, 
as in Fig. 2. Moreover, it can be further improved by increasing the 
number of sections or by using special teeth, as in Ref. 5. 

Equation (1) was derived in eq. (2) by a perturbation analysis that 
assumes that the difference 8) — 6 is small. Recently, p, has been 
derived® without this restriction for the special case where one of the 
two waveguides is uncorrugated, as in Fig. 1, and a more general 
derivation is given in Ref. 7. By letting 8) — 6 — 0 in the results of 
Refs. 6 and 7, we obtain eq. (1). For the purpose of this analysis, 
however, eq. (1) is more than adequate since all reflections, p;, are 
small. 

Reflections of less than —30 dB were reported in Ref. 8 over a band 
2.7 < ka < 3.8, using only five slots. A mode-matching technique was 
used to accurately determine the transformer scattering matrix. Here, 
however, we are interested in a ratio w;/wo, which is much greater than 
3.8/2.7. Experimental and theoretical results on the problem of mini- 
mizing the input reflection of a feed using slots of varying depth were 
reported in Ref. 9, and excellent performance was obtained using ring- 
loaded slots, which greatly reduced excitation of the HE,. mode. We 
also mention the very low cross-polarization levels obtained in Ref. 10 
using a curved-aperture corrugated horn. 


Il. APPROXIMATE DERIVATION OF 4, s; 


In this section we minimize the total reflection, p;, in Fig. 2 in the 
vicinity of wo. We assume that only the HE;, mode propagates for 
z> 0 and that the TE,; mode is incident from the left. Let M; be the 
number of corrugations in the ith section, and let \,; be the wavelength 
for the HE,; mode. The reflection, p;, at the ith junction will be 
determined approximately using eq. (1), which is strictly valid only if 
pi is small and 


M; > 1, (3) 


implying h; < \,;, where h; is the teeth separation in the ith section. 
The propagation constant, 8, of a mode is related to its transverse 
wave number, o, 


Ba = v(ka)*? — u’, (4) 
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where u = oa and k is the free-space propagation constant. For the 
TEn1 mode, 


U = Up = 1.8411. (5) 


For the HE,; mode in a corrugated waveguide, u is determined by s/h 
and d, and it can be calculated as in Refs. 1 and 3. 

Assume the reflections, p;, in Fig. 2 are small, and let the frequency 
dependence of p; and 6;/k be neglected. Then, for a maximally flat 
passband characteristic,* 


_1_™! 
Pe> 8 (N= Til ©” 

where N is the number of sections and 
pr = & pi. (7) 


The reflection, p;, is related® to the propagation constants §;-, and 8; 
of the waveguides at the ith junction, 

Bi-1 — Bi: 

pi = — > (8) 

Bi-1 + B; 
and, therefore, using eqs. (4), (6), and (8) one can determine Wy, Up, 
etc., if WU) and uy+; are given. Once the u; are known, the values of 
s;/h; can be determined approximately as in Refs. 1 and 3. Notice the 
length, 4 of the ith section must be chosen so that 


Nei 
= a ? at Wo. (9) 


(6) 


In the following section the corrugated waveguide at the input of 
the feed (see Fig. 1) will be characterized by 


b h 
—= 1.789, —= 1.256, (10) 
a Ss 


which can be shown to give 


ka = 4.08, at w =a). (11) 
Then, if one chooses 
wo = .63 w), (12) 
one finds 
Un = Uz = 2.114, for w= wo. (13) 


Therefore, eqs. (4), (6), and (8) give 
u, = 1.928, Ug = 2.056, for w = wo. 
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This requires 


tee 16.8, ee 3.15, (14) 


81 S2 


taking into account that up = 1.841 in the input guide. 


Ill. EXPERIMENTAL RESULTS 


The values of eq. (14) will not produce exactly a maximally flat 
characteristic because of the frequency dependence of p; and £;/k. 
Furthermore, in the experiment each section was realized using only 
four corrugations, in which case the values calculated for p; and u; in 
the preceding section for w = wo are not expected to be accurate since 
condition (3) is violated (a consequence of this is pointed out in the 
appendix of Ref. 3). Thus, taking into account the difficulty involved 
in accurately determining the dependence of p; on h;/s;, the optimum 
values of h;/s; were determined experimentally, and they are 


h h 
+= 15.23, — = 2.492, 
Si S2 
assuming the conditions in (10). 

The experiment was carried out choosing 


a = 0.651”, d = 0.514”, (15) 


which implies wo = 7.18 GHz and w,; = 11.455-GHz. The two sections 
were constructed using eight rings assembled as indicated in Fig. 2. 
The output corrugated waveguide was realized by the technique of 
Ref. 2. The input circular waveguide was connected to a rectangular 
waveguide using a long transition, and the input reflection, p;, was 
measured as in Fig. 3, using a commercial directional coupler. The 
measured reflection, shown in Figs. 4(a), (b), and (c), includes a 
component due to small reflections from the directional coupler and 
the transition. This component, of about —45 dB, is responsible for 
the fast ripples in Fig. 4. Also shown in Fig. 4 are the values of p;, 
calculated using eq. (8), for h;/s; given by eq. (8), assuming condition 


TRANSITION 
Pt FROM CJ TOO 


MATCHING 
SECTION 


Lege ieomescet 


Fig. 3—Apparatus used to measure the reflection, p:. 
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10.5 11.0 @W, 11.5 12.0 
GIGAHERTZ 
Fig. 4—Reflection, p,, measured with and without transformer from (a) ka = 2.4 to 
3.0, (b) ka = 3.0 to 3.7, and (c) ka = 3.7 to 4.2. 
(3) but without neglecting the frequency dependence of 6;/k and p;. 
There is considerable discrepancy for w = wo, as expected in view of 
the relatively large values of h;/,;. Also shown in Fig. 4 for comparison 
is the reflection p, measured without matching section. 
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Fig. 5—The reflection at the throat of a horn can be matched by a suitable quarter- 
wave transformer. 


IV. CONCLUSIONS 


The input reflection, p;, in Fig. 4 is less than about —30 dB for 
values of ka in the band (2.43, 4.22) and less than —35 dB over the 
two bands (2.43, 2.766) and (3.69, 4.22). Over most of these two bands, 
|p:|? is less than —40 dB. The transformer is simple to realize, 
consisting of eight rings assembled as shown in Fig. 2. 

Over part of the above frequency range, the input waveguide can be 
shown to be multimoding because of the TM; mode, which propagates 
for ka > 3.8317. Also, the output waveguide is multimoding over 
approximately the same frequency range. However, one can show, 
using the approximate formulae of Ref. 3, that the transformer greatly 
reduces generation of this mode, which should therefore be negligible. 
For this reason, and also because of the difficulty in precise measure- 
ments of the conversion coefficients, the mode in question was ignored 
here. 

Finally, in a feed, the input reflection includes a contribution from 
the throat of the horn.®”° This contribution is determined by the angle 
y in Fig. 5, and it can be eliminated by using a suitable matching 
transformer as shown in Fig. 5. This is pointed out in Ref. 9. 
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Moment Calculations by Digital Filters 


By Z. L. BUDRIKIS* and M. HATAMIANt 
(Manuscript received May 25, 1983) 


We present a simple recursive algorithm for computing moments of two- 
dimensional integer arrays. It uses only additions and can be implemented for 
high-speed and real-time computation at video rates. We describe Comple- 
mentary Metal-Oxide Semiconductor (CMOS), Very Large-Scale Integrated 
(VLSI) implementation of the algorithm in a single chip that can calculate 
the 16 moments yp’ (i, j = 0, 1, 2, 3) (i.e., up to the sixth-order moment) on 
512 X 512 array of 8-bit integers in real time (at video rate). Such a chip can 
have potential applications in image processing, graphics, and robotics. The 
basic building block of the system is a single-pole digital filter that is imple- 
mented by recursive addition. The complexities involved in designing the chip, 
as well as its area, are significantly reduced by taking advantage of the fact 
that the column samples of the data array can be processed at a much slower 
rate than the row samples. An estimate of the chip area obtained from the 
layout design of the individual cells is given. 


I. INTRODUCTION 


Moments are familiar from statistics and mechanics, and are finding 
applications in other areas, among them video processing.” The mean, 
or first moment, is a particularly robust estimator of the ‘center’ of a 
distribution. Similarly, the centroid is a good single indication of the 
location of an extended object. We want to use it for pinpointing the 
position of a light pen.° 

A possible brake on calculating moments, particularly in real-time 
applications, may well have been the computational burden that they 


* University of Western Australia. ‘AT&T Bell Laboratories. 
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appear to entail. It would seem that the calculation calls for exponen- 
tiations and multiplications, as well as additions. However, given real- 
world measurements that are finite in precision, range, and number, 
and are at regular intervals, the calculation of moments can involve 
much less effort than at first seems necessary. 

One scheme has already been reported* in which the multiplying 
coefficients are generated without exponentiation and in which the 
multiplications are distributed over the individual bits of the measure- 
ments and are therefore accomplished by single AND gates. An even 
simpler approach is possible, requiring less computation by an order 
of magnitude, and is reported here. 

Our approach comes about from recognizing that Infinite Impulse 
Response (IIR) digital filters with impulse responses h(n) = u(n), or 
nu(n) or n?u(n) or n'u(n), with u(n) the unit step sequences, will 
calculate respectively the zeroth-, first-, second-, and ith-order mo- 
ments of their input sequence. It is a straightforward matter to extend 
this to arrays in the plane and higher dimensionality, since the 
computations along the different dimensions are separable from each 
other. 

The poles of the required filters can be very easily realized by first- 
order sections, and these are nothing more than recursions with unity 
feedback. Thus, if only the poles, and not zeros, of these filters are 
realized, the computation between successive sample points is just one 
addition and can be instrumented to take place in real time on picture 
data from standard television. The basic building block of the system, 
then, will be a first-order stage, and the resulting filter will have a 
highly regular and recursive structure that is extremely attractive for 
Very Large-Scale Integrated (VLSI) implementation. 

As will be seen in Section III, the consequence of discarding the 
zeros of these filters is that the output of, say, the mth stage, rather 
than being the mth-order moment itself, is a linear combination of 
first through mth-order moments. This, however, does not create any 
major difficulty because the coefficients of this linear combination are 
exactly known, and therefore starting from zeroth-order moment, at 
each stage the effect of previous moments can be removed by subtrac- 
tion in a post-processing step. 

We describe a design for a custom Integrated Circuit (IC) in Com- 
plementary Metal-Oxide Semiconductor (CMOS) technology that can 
calculate the 16 moments, p'” (i, j] = 0, 1, 2, 3), of a 512 x 512, 8 bits/ 
pixel image in real time (i.e., at conventional video rate). 


lil. THE TASK 


Given the array x(n, m)(n = 0,1, --- N, m=0,1, --- M) of integer 
values, the moments about m = 0, n = 0 are 
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Caen 
po! = n'm’x(n, m) 


n=0 m=0 


n=0 m=% 


Nn .[M. 
>n ( Ly mx(n, m) (1) 


Only the innermost computations, corresponding to along line in 
television, have to be at speed. The outer computations are at lower 
rate by a factor of M. The choice of n = 0, m = 0 as the pivot for the 
moments is arbitrary. Well-known expressions relate the moments 
about one pivot to those about another. The pivot point about which 
we chose to do the computations is n = N, m= M, ie., 


fas. _M 
py = = (N — n)’ z (M — m)x(n, m). (2) 


il. THE ALGORITHM 


We first consider the one-dimensional case. Consider a sequence 
x(n) n=0,1, --- N; if this sequence is applied to the input of a digital 
filter with impulse response h(n) = n‘u(n), then from the well-known 
convolution theorem the output will be 


HONS Sah) 
k=0 


N . 
> x(k)(n — kt (3) 
k=0 
If this output is evaluated at n = N, then we will have 
N ; 
y(N) = Rs x(k)(N — )’, (4) 


which is the ith-order moment of x(n) about the point n = N. Figure 
1 shows the filter for i = 0. This is a single-pole filter with the transfer 
function 1/z — 1 in z-transform domain. Its implementation is trivial, 
just an adder with one sample delay and a feedback, or basically an 
accumulator. For higher-order moments (i = 1, 2, ---) the filter will 
have both poles and zeros. In general, for the ith-order moment we 
have i zeros and a pole of order i + 1 located at z = 1. The first and 


x(nT) y(nT) 





y(NT)=p° 


Fig. 1—Single-pole filter for generating zeroth-order moment. 
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Table |—Impulse response and transfer function of moment filters 











All-Pole 
Moment Impulse Transfer All-Pole Impulse 
Order, Response, Transfer Function, Function, Response, 
I h(n) H(z) H(z) h(n) 
1 1 
0 u(n) a ar u(n) 
z 1 
ae Sees -1 : 
1 nu(n) @—- 1? @-1? (n-1)u(n-1) 
2 z(z + 1) mer ean 1/2 (n — 2)°u(n — 2) 
on) (2-1) (2-1 +1/2(n-2)u(n—-2) 
2 1/6 (n — 3)'u(n — 8) 
3 n?u(n) a eS) E 7 + 1/2 (n — 3)"u(n — 8) 
(z— 1) (2 — 1) + 1/3 (n — 3)u(n — 3) 


second columns of Table I, respectively, show the impulse response 
and transfer function of the moment filters for i = 0, 1, 2, 3. If we 
discard the zeros of the transfer function, then for the ith-order 
moment we obtain the all-pole filter 


_ 1 
(z = Ly 


This filter can be very easily implemented by cascading i + 1 first- 
order stages similar to the one shown in Fig. 1. Obviously, the output 
of the filter without zeros will no longer be just the ith-order moment 
of its input signal except for 1 = 0. We examine this next. 

Starting from i = 0 with the z-transform pair, 


H(z) = (5) 


Hole) = = © hiln) = un), 


and, using the differentiation property of the z transform, we can 
progressively derive the impulse responses h;(n) corresponding to 
H(z) = 1/(z — 1)'*? for all values of i. The results for i = 0, 1, 2, 3 are 
shown in the fourth column of Table I. To see the effect of zero 
elimination, let us consider the case for i = 3 as an example. From 
Table I we have: 


ha(n) = ; (n — 3)3u(n — 3) + ; (n — 3)2u(n — 3) 


+ ’ (n — 3)u(n — 3). (6) 


Convolving this impulse response with the input x(n) and evaluating 
the output at n = N + 8 will result in: 
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N 1N 
yalN + 3) = —& 2 xlkN — ky? + 9 x(k)(N — k)? 


Oot 


1 N 
+— J x(k)(N — k) 
3 k=0 


ae 1 
=ap tiptoe 7 
gu tage tae, (7) 
which is a linear combination of first-, second-, and third-order mo- 
ments. Similarly, for an arbitrary value of i, the output of the filter 
evaluated at n = N + i will be a linear combination of first through 
ith-order moments, i.e.: 


yi(N +i) = D heute (8) 
=1 


Conversely, the moments yp’ (i = 0, 1, -- -) can be expressed as a linear 
combination of the filter outputs y;(N+7); in matrix notation, 


> 


p= CY, (9) 
where 
we yo(N) 
~_ tH ~ | (N+ 1) 
B= Tey *= Vy + 2) 


The coefficient matrix C can be easily derived from the impulse 
responses h;(n). For i = 5, this matrix is: 


0 0 0 0 
1 0 0 0 
= 2 0 0 


1 6 6 0 (10) 


—-1 14 -—36 24 
1 -30 150 -240 120 


oooocr 
oooo eo 


Figure 2 shows the complete filter structure for generating zeroth- 
through third-order moments. The delay introduced on the output of 
each stage is necessary for synchronizing the output before performing 
the matrix operation. 

So far in this section we have been discussing the moment calcula- 
tion problem for a one-dimensional sequence x(n); extending the 
algorithm to two or higher dimensions is a straightforward matter. 
This is true because the computations along the different directions 
are separable. Consider an array of integers x(n, m) (n = 0,1, --- N, 
m= 0,1, --- M), which can, for example, be the digitized samples 
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Fig. 2—Moment filter for generating zeroth- through third-order moments. 


of a video image. As mentioned in Section II, the moments p™ [i.e., 
(i + j)th order] of this array can be.computed from 


i N. . M . 
py = > (N—-n)’ > (M — m)’x(n, m) 
n=0 m=0 


2 (N= n)y(n). (11) 


The inner summation, y;(n), represents the jth-order moment of the 
nth row (“along-line moment”) which, as discussed previously, can be 
computed by filtering the samples on that row with a filter with the 
impulse response m/u(m) (row filter) and evaluating the output at 
m=M. 

Now p” is equal to the ith-order moment of the sequence y;(n) and 
can be computed using a filter with impulse response n‘u(n) (column 
filter). Figure 3 shows the filter block diagram for generating p’’. The 
input to the column filter implementing eq. (11) is updated every 
M samples of the input signal, corresponding to one row of the array 
x(n, m); the row filter is also reset at this time. If the array x(n, m) 
represents the digitized samples of a video image, then for real-time 
operation, the sampling rate will be equal to the pixel frequency for 
the row filter and line frequency for the column filter. As we will see 
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in Section IV, this fact greatly reduces the complexity of implemen- 
tation. 

The column filters can be realized with single-pole stages in exactly 
the same manner as was previously discussed in this section. There 
are two alternatives. One is to use the filter structure shown in Fig. 2, 
which includes the dematrixing operation, for the column filters as 
well and make it a building block for implementing the two-dimen- 
sional filter. The block diagram of Fig. 4 illustrates this approach. In 
this case the outputs of the two-dimensional filter will be the true 
moments yp’, The other alternative is to use the original single-pole 
stage shown in Fig. 1 as the building block for constructing the two- 
dimensional filter. In this case the outputs of the filter will be a linear 
combination of the moments; to obtain the actual moments p, a 
dematrixing operation has to be done at the end of the process (i.e., 
after the last row of the input array has passed through the filter). 
Figure 5 shows the complete two-dimensional structure implementing 
this latter filter configuration. From the implementation point of view 
this structure is more attractive than the former realization. (More on 
this in the following section.) The relationship between the array 
x(n, m) in Fig. 5 and the input signal is determined by the sampling 
strategy. For example, if f(t) is a sequentially scanned signal to be 
filtered, then 


m=M n=N 


x(nm) hAtm)=m/U(m) O O A(n)=n/U(n) Oo CO pli 


Fig. 3—Block diagram of the filter for generating p’. 


ROW FILTER 


x(nm) 











COLUMN FILTER 





COLUMN FILTER 





COLUMN FILTER 
COLUMN FILTER 


Fig. 4—Moment filter for generating yp (i,j = 0, 1, 2, 3) using the filter shown in Fig. 
2 as a building block. 
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x(n,m) 









(N+1)T" 


(N+2)T' 
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Fig. 5—Moment filter for generating y" (i, j = 0, 1, 2, 3) using the single-pole stage 
as a building block; T’ = MT. 


x(n, m) = f(mT + nT’), 
where T is the pixel period and T’” is the line period. 


IV. VLSI IMPLEMENTATION 


In this section we describe VLSI design for implementation of our 
two-dimensional filtering algorithm, which can calculate the moments 
uw (i = 0, 1, 2, 3,7 = 0, 1, 2, 3) of a 512 X 512, 8 bits/pixel image in 
real time (i.e., at conventional video rate). We consider a maximum 
chip size of 6 mm X 6 mm. To meet the speed requirement the picture 
element (pel) processing time should be less than 125 ns. 

We examine the two filter structures shown in Figs. 4 and 5 for 
possible VLSI implementation. The boxes labeled ‘row filter’ and 
‘column filter’ in Fig. 4 could be made identical, and each could be 
realized by one chip. Then the two-dimensional filter would be con- 
structed by putting five of these chips together. This may at first seem 
a reasonable approach, but in fact, it is not: the speed requirements 
are set by the calculations that are performed in the ‘row filter’, while 
the size of the accumulators is dictated by the ‘column filters’. The 
result would be very large accumulators (up to 64 bits) operating at 
the high rate of at least 8 MHz. Moreover, that implementation would 
preclude any possibility of having all calculations done on a single 
chip. 

With the structure of Fig. 5, all mismatches of speed and size can 
be avoided. What is more, the whole system can be put on a single 
chip with an area occupancy only slightly greater than needed for a 
single filter of Fig. 4. The only shortcoming of the system of Fig. 5 
compared to that of Fig. 4 is that it is the poles-only realization and 
requires external dematrixing. That, however, is only minor, for those 
calculations involve only a few multiplications and additions, and can 
be readily performed by the host processor that would use the com- 
puted moments. 

The basic building block of the filter in Fig. 5 is an accumulator 
functioning as a single-pole digital filter. The maximum word lengths 
that may appear as outputs and hence determine minimum accumu- 
lator sizes increase with order. For 8-bit input data samples in a 512 
X 512 array, the accumulations along the row section may progressively 
reach 17, 25, 33, and 41 bits, and those in the last column section, 
progressively 48, 55, 62, and 68 bits. 

To achieve the required speed in the row section, we have decided 
to use look-ahead carry adders. The column sections could be imple- 
mented with simple ripple carry adders occupying much less area and 
obviously making it possible to accommodate all circuits on the avail- 
able chip area. In fact, we do even better than that by using serial 
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additions. Conversion to serial streams when passing data from the 
top row into the columns is in any case indicated by wiring consider- 
ations. 

If given parallel operations in both row and columns, there would 
be inordinate numbers of wires that would have to pass down vertically, 
presenting a difficult routing problem. To alleviate this, we decided to 
use a parallel-to-serial convertor on the output of each stage of the 
row filter and transfer serially the outputs to the column heads. Note 
that the along-row and down-to-column transfers take place at differ- 
ent times, the latter only once every line during the horizontal- 
blanking interval, and that the horizontal connections remain in 
parallel without any sacrifice in speed. To achieve the dual feed-out, 
the output registers of the along-row accumulators are designed as 
parallel in parallel/serial out shift registers. 

Once given serial streams into the columns, there is no point in 
converting back to parallel operation or using ripple carry adders: the 
accumulations can be done equally well serially, using for each stage 
a one-bit full adder and a shift register. The block diagram of Fig. 6 
illustrates this final design. The upper portion of this figure shows the 
four stages of the row filter. The lower portion is divided into four 
sections that correspond to the four column filters and receive their 
input from the serial output of the shift registers of the row-filter 
accumulators. 


[" tSTSTAGE ~ 2NDSTAGE  3RDSTAGE  4THSTAGE | aaWelbrER 
cc ira oc oc "4 
E 2 ra 
a o a a] | 
INPUT 17-BIT |2 25-BIT |9 33-BIT |9 41-BIT |2 
DATA ADDER | cc ADDER | cx ADDER | & AppDER|«| | 
b i - - 
| a bo = we] | 
| B a w o l 
1S =e! 7s ee ee I 
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FOUR SERIAL STAGES 
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ce | 
SS 
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1-BIT 64- BIT 
FULL ADDER SHIFT REGISTER 


_ Fig. 6—Block diagram of a VLSI design for single-chip implementation of the two- 
dimensional moment-generating filter. 
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Fig. 7—Floor plan and input/output pins for the design shown in Fig. 6. 


Each section consists of four stages stacked on top of each other, 
each containing a one-bit full adder (the dark area) and a 64-bit shift 
register composing a 64-bit serial accumulator. The feedback path for 
the accumulator is not shown in the figure. 

The design of the chip is not complete yet, but a conservative 
estimate of the size has been obtained from the design of the individual 
cells. The design is for CMOS technology with a two-micrometer 
design rule. Figure 7 shows the result of the size estimate in a floor 
plan block diagram drawn in proportion to the 6 mm X 6 mm boundary. 
Each section of the floor plan is drawn 15-percent larger in both 
dimensions than the layout size estimated from individual cell designs. 
As can be seen from Fig. 7, there is still plenty of silicon area left for 
adding the timing and control section and, if desired, one or two more 
column stages for generating higher-order moments. The anticipated 
input/output pins for this chip are also shown in Fig. 7. The chip 
operation is synchronized with the input image data through three 
synch pulses: End of Pixel (EOP), End of Line (EOL), and End of 
Frame (EOF). After the output data ready signal is asserted, the host 
processor should read the 16 outputs, y®°, y®!, - - - y*%, in serial format. 
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Vv. CONCLUSION 


We have described a simple and powerful algorithm for computing 
moments of a two-dimensional array of integers using digital filters. 
The algorithm involves only addition operation and is very suitable 
for real-time implementation. A two-dimensional digital filter with a 
separable impulse response, h(n, m) = n'u(n)m/u(m), can generate the 
(i + j)th order moment of its input. Realization of an all-pole version 
of such a filter results in a filter with a highly regular structure capable 
of producing all the moments, from zeroth to (i + j)th order, of a two- 
dimensional input data array. We have shown that the consequence 
of eliminating the zeros of the filter and using an all-pole version is 
that the outputs will be a linear combination of the moments rather 
than the direct moments themselves. But this does not create any 
difficulty because the constant coefficients of the linear combination 
are known and a simple dematrixing operation at the end of the 
process will recover the actual moments. The basic building block of 
our moment filter is a single-pole digital filter that can be easily 
implemented by an accumulator. We have proposed a VLSI design for 
single chip implementation of this moment calculation algorithm. The 
preliminary results obtained from the layout design of the individual 
cells for 2-micrometer CMOS technology indicate that on a single 6 
mm X 6 mm chip, we can certainly implement a filter for generating 
the 16 moments p (i, j = 0, 1, 2, 3) of a 512 x 512, 8 bits/pixel image 
in real time. 
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Wideband Operation of Nonlinear Solid-State 
Power Amplifiers—Comparisons of Calculations 
and Measurements 
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Calculations of intermodulation (IM) noise and input/output power-trans- 
fer characteristics for a number of different nonlinear Solid-State Power 
Amplifiers (SSPAs) are found to be in very good agreement with measure- 
ments. The calculations are based on the measured AM/AM and AM/PM 
characteristics, and include results of a computer simulation, and analytical 
results based on modeling the SSPA as an ideal envelope limiter. The results 
demonstrate that the calculations predict wideband performance characteris- 
tics of nonlinear amplifiers for arbitrary input signal configurations and 
operating conditions, as well as provide a basis for comparison against which 
anomalous behavior can be identified. Earlier results, in which anomalous 
behavior of an SSPA was observed to result in significant performance 
degradation, are discussed in the light of these new results. A simple analytical 
criterion is developed for evaluating measurements to identify such degrada- 
tion. The effect of AM/PM conversion is investigated and found to cause 
significant degradation in IM performance only at low power levels for typical 
SSPAs. In the vicinity of saturation, and above, the contribution to IM noise 
due to AM/PM conversion is found to be very small in comparison with that 
due to the saturating AM/AM characteristic. 


I. INTRODUCTION 


Two important considerations that arise in the design of nonlinear 
power amplifiers such as those used in satellite and terrestrial radio 
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systems are the input/output power-transfer efficiency, and the 
amount of intermodulation (IM) noise produced as a function of input 
and output power levels. The former of these is related to what is 
termed the power-added efficiency of the device and the latter is a 
measure of its range of linearity. An earlier paper presented some 
results of analyses of these aspects of a nonlinear multistage Solid- 
State Power Amplifier (SSPA) under wideband excitation.’ The re- 
sults revealed a significant discrepancy between the calculated and 
measured power-transfer characteristics of the SSPA. In particular, 
measured output power levels in the vicinity of saturation were found 
to require approximately 3 dB more input power above predicted 
values based on AM/AM and AM/PM measurements. Ideally, the 
measured single-tone AM/AM and AM/PM characteristics of a non- 
linear amplifier uniquely specify how it will perform in operation with 
wideband signals. The results of Ref. 1, however, have shown that 
anomalous cases—in which the single-tone measurements are not 
sufficient for predicting measured wideband behavior—can occur, with 
significant degradation in performance. It is therefore necessary to 
have a basis for comparison against which measurements of wideband 
power-transfer efficiency, IM performance, and possibly other per- 
formance criteria can be evaluated. In following up these initial results, 
the techniques originally developed in Ref. 1 have been extended and 
applied to a number of additional nonlinear SSPAs with newly mea- 
sured performance data.” As we will see from the results presented 
below, calculations of this kind provide such a basis for any arbitrary 
input signal configuration and, essentially, any set of operating con- 
ditions. 

Numerical results are obtained by means of computer simulation 
and an analytical closed-form expression for the wideband power- 
transfer characteristics. In the analysis, the multicarrier signal is 
assumed to be Gaussian, which is a very good approximation for the 
signals of interest here. In all cases it is assumed that the devices 
under consideration are memoryless or instantaneous. The initial 
results in Ref. 1 were obtained for a multistage SSPA, referred to as 
MS-1. The additional SSPAs for which comparisons of calculations 
and measurements have been made include a single-stage SSPA, 
referred to as SS-1, of which six such stages were employed in MS-1, 
another multistage SSPA, referred to as MS-2, which also employed 
six SS-1’s, and three additional multistage SSPAs. 

In all cases, measured and calculated results presented below for 
these SSPAs are in very close agreement, differing for the most part 
by fractions of a decibel. This demonstrates that the methods pre- 
sented herein can be used to accurately determine what the operating 
characteristics of nonlinear SSPAs under wideband excitation should 
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be, and that they are sufficiently reliable to reveal anomalous behavior 
and to help pinpoint the sources of such anomalies. Specifically, the 
anomalous behavior of MS-1 was identified by comparison of the 
measured power-transfer characteristics with calculated values. Fur- 
thermore, the good agreement between measurements and calculations 
for SS-1 and MS-2 shows that, in general, the transistors employed in 
MS-1 can be expected to perform as predicted both singly as well as 
in a multistage configuration. Thus, the anomalous behavior observed 
with MS-1 has been clearly identified as being a peculiarity of the 
particular device. Possible reasons for this include one or more mal- 
functioning individual stages, and peculiarities in the overall design 
configuration, including possible saturation of early stages. 

In the above discussion we have used the term “anomalous behavior” 
to indicate an observed discrepancy between measurements and pre- 
dictions based on AM/AM and AM/PM characteristics. It is of interest 
to consider possible connections between predictable performance and 
optimal design. Kaye, George, and Eric have conjectured, but not 
proved, that, from the point of view of minimizing IM noise as a 
function of output power, the optimum nonlinear power amplifier has 
AM/AM characteristics that exhibit exact linearity up to the satura- 
tion level, with constant power output thereafter, and zero AM/PM 
conversion.® These characteristics are referred to as those of an Ideal 
Envelope Limiter (IEL). Reference 4 proves part of Kaye, George, and 
Eric’s conjecture, namely that AM/PM conversion can only degrade, 
never improve, IM performance for multicarrier input signals. Since, 
in practice, some degree of AM/PM conversion can be expected, it is 
of interest to investigate how the shape of the AM/PM characteristic 
affects the extent of the degradation. 

For SSPA’s, the AM/PM measurements examined thus far indicate 
that the phase shifts at low power levels are relatively small (i.e., a 
few degrees), and increase very gradually up to saturation. At this 
point they abruptly begin to increase much more rapidly, reaching 
perhaps 25 to 30 degrees a few decibels above saturation; an example 
of such a characteristic is presented below. The larger and more rapidly 
increasing phase shifts at and above saturation, however, are less 
important to the overall IM performance than the phase shifts at 
lower power levels (e.g., ~6 to 10 dB back-off). This has been deter- 
mined by calculating the IM performance of the SSPA with and 
without AM/PM conversion, which is easily done in the computer 
simulation. It is found near saturation that the AM/AM characteristic 
dominates the IM production to the extent that the difference with 
and without AM/PM conversion is of the order of 1 dB or less. On 
the other hand, for low power levels, because the AM/AM character- 
istic seen by the signal is effectively much more linear, relatively small 
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phase variations can make a very significant difference in IM perform- 
ance—up to 6 dB is observed here. This indicates that, assuming that 
some degree of AM/PM conversion cannot be avoided, the design goal 
would be to have it as small as possible below saturation, while in the 
vicinity of saturation and above, much greater leeway in the shape of 
the characteristics could be tolerated. This could be of practical 
significance in the design of networks for compensating for AM/PM 
conversion. That is, the results indicate that it is more important to 
compensate for the small, slowly varying phase shifts below saturation 
than for the large, rapidly varying phase shifts at and above saturation, 
where the compensation is presumably more difficult and where re- 
quirements on the performance of phase-compensation networks could 
be relaxed with little degradation in IM performance. 

With regard to the ideal AM/AM characteristics, to date no viola- 
tions or counter examples to Kaye, George, and Eric’s conjecture seem 
to have been found. AM/PM conversion has no effect on input/output 
total power-transfer characteristics.* Furthermore, based on compari- 
sons of measurements with eq. (6), below, it is found that, excluding 
anomalous cases, SSPAs exhibit power-transfer characteristics that 
are nearly identical to those of an IEL, which under Kaye, George, 
and Eric’s conjecture can be considered optimal. Thus, it should be 
possible for the AM/PM to be adjusted or compensated as necessary 
for minimum IM noise without affecting or sacrificing ideal power- 
transfer efficiency. With the use of eq. (6), a simple analytical criterion 
is developed for evaluating the extent to which SSPA wideband power- 
transfer measurements meet this ideal performance. 

Section II summarizes the analytical results relevant to this work, 
including the optimal power-transfer characteristics of an IEL. Section 
III describes the computer simulation, Section IV compares measure- 
ments and calculations, and a summary of results and conclusions is 
presented in Section V. 


Il. SUMMARY OF ANALYTICAL RESULTS 


The multicarrier signals of interest here can, to a very good approx- 
imation, be represented as bandpass Gaussian noise in the form 


n(t) = x(t)cos wot + y(t)sin wot, (1) 


where w is the filter center frequency. Equation (1) is the narrowband 
representation of Gaussian noise, in which the low-pass functions x(t) 
and y(t) are required to be slowly varying with respect to wo. In the 
cases of interest (e.g., satellite communications), the bandwidth of x(t) 
and y(t) is nominally 40 MHz and fp is 4 GHz, satisfying this require- 
ment by a very large margin. Equation (1) can also be written as 
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n(t) = r(t)cos(wot + 6(¢)), (2) 


where r(t) = (x?(t) + y?(t))”? is the envelope, and 6(t) = tan7'y(t)/ 
x(t) is the phase. Since x(t) and y(t) are Gaussian, r(t) has a Rayleigh 
distribution: 


= eM l20? (3) 
oO 
where o? = E(x?(t)) = E(y7(t)) is the input-signal power. 

The output of the nonlinearity will include an infinite number of 
harmonics of wo. We are interested only in the principal zone, including 
only the first harmonic and a neighboring frequency range sufficiently 
wide to encompass all the IM noise produced (~200 MHz) (see Section 
III). For an instantaneous system, the principal-zone output can be 
written as: 


A(r)cos(wot + 6(t) + v(r)), (4) 


where A(r) and v(r) vary in time as r(t), and, as functions of r, 
represent the AM/AM and AM/PM characteristics of the nonlinear- 
ity. 

As Ref. 4 shows, the total output power in the principal zone is just 
1/2(A(r)), independent of v(r). For all the SSPAs studied, the 
AM/AM characteristics very closely resemble those of an IEL and 
excluding anomalous cases, the power-transfer characteristics are also 
found to be essentially identical to those of an IEL. To calculate the 
output power for an JEL, referring to Fig. 1 we write: 


(3 aun) = J E(A%r)/2)P(rydr + | E(A?(r)/2)P(r)dr, (5) 


where r, is the value of the input envelope r necessary to drive the 
limiter into saturation and P(r) is given by (38). In the saturated region, 
r>r,, we have A(r) = rm, the maximum saturated value of the output 
envelope, and in the linear region, r < r,, we have A(r) = kr, where 
the constant k is r,,/rs. Under these conditions the integrals in (5) are 
easily evaluated and, denoting the output power, 1/2(A7(r)), as Po, 
and P, = 1/2r?, Pm = 1/2r2, it is found that 


1 —e. 3/ Pi 


P./P; : 


where in (6) we have redesignated the input power as P;. 

Equation (6) provides a convenient means for quickly evaluating 
SSPA wideband power-transfer measurements in order to determine 
whether the device is performing as it should. As diagramed in Fig. 1, 
the IEL equivalent to a given SSPA is defined by extending the linear 
and saturated portions of the AM/AM characteristic to the point of 


Po = Pm (6) 
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Fig. 1—SSPA AM/AM characteristics with equivalent IEL and IEL power-transfer 
characteristics. 


intersection that defines P, and P,,. Using these values of P, and Py, 
the theoretically correct values of Po as given by (6) can be easily 
plotted by noting that, starting with a value of P; 6 dB above P, and 
stepping down in 3-dB increments to 6 dB below P,, the corresponding 
values of output power are very nearly 0.5, 1, 2, 3.6, and 6 dB below 
Pm; the exact values are given in the table insert in Fig. 1. Using this 
method, measurements of power-transfer characteristics for a number 
of SSPAs are compared in Section IV to theoretical values and found 
to be in almost exact agreement, which demonstrates the usefulness 
of this approach. It should be noted that this method is not valid for 
traveling-wave tubes, which do not have IEL-like AM/AM character- 
istics. 


Ili. COMPUTER SIMULATION 


Figure 2 presents a block diagram of the computer simulation used 
in these analyses, and Ref. 1 gives a description of the calculation 
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Fig. 2—Block diagram of simulation. 


procedures. In the laboratory measurements, the input waveform 
consisted of 24 equal-power, filtered, Quadrature Phase-Shift Keying 
(QPSK) signals spaced on either side of an unmodulated carrier at 
band center (4 GHz) by the frequency increments shown in Table I 
(from Ref. 2). The power in the unmodulated tone was adjusted at the 
input to be equal to that of a single QPSK signal. 

In simulating this waveform, a complex baseband formulation was 
employed. The simulated input consists of the outputs of 48 (24 pairs 
of) independent, equal-power, Pseudo-Noise (PN) sequence genera- 
tors. The outputs of each pair of generators are formed into a complex 
number, which, in the baseband formulation, represents a QPSK 
signal, with the real part representing the in-phase component and 
the imaginary part, the quadrature component. Each QPSK signal is 
then passed through a four-pole Butterworth filter with a 3-dB (posi- 
tive-frequency) bandwidth of 0.5 MHz. This is essentially the same 
filtering employed at the output of each QPSK signal generator in the 
24-channel laboratory test set. The pulse duration employed in the 
simulation is 1.3 ws, corresponding to 0.772 X 10° symbols/s, or a DS- 
1 rate of 1.544 x 108 bits/s. 

After filtering, each QPSK signal is shifted in frequency on either 
side of zero by the same increments used in the measurements (see 
Table I). The 4-GHz carrier in the laboratory measurements is rep- 
resented in the simulation by a dc signal adjusted in power to be equal 
at the input to a single QPSK signal. The spectrum of the input signal 
employed in the simulation is shown in Fig. 3. 

This signal is then input to the nonlinear system under considera- 
tion. The nonlinearity is specified by entering tables of the measured 
AM/AM and AM/PM values into the program to which the simulation 
fits continuous curves. An example of the spectrum at the output of 
the nonlinearity, for input power in the vicinity of saturation, is 
presented in Fig. 4. For evaluating IM performance, the laboratory 


Table I—Frequencies used in simulation 


Frequencies Frequencies 
Channel Channel (MHz) 
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Fig. 3—Input signal spectrum showing 24 T-1 signal plus unmodulated tone at band 
center. 
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Fig. 4—MS-2 output spectrum with amplifier saturated. 


NONLINEAR AMPLIFIERS 239 


measurements employed a narrowband filter in the notch at 4-GHz 
for measuring the output unmodulated carrier power, and a number 
of narrowband filters located in the notch between the 4-GHz carrier 
and the neighboring QPSK signals for measuring IM power; the 
individual IM-power measurements were averaged. Correspondingly, 
the simulation employs a narrowband filter at zero frequency for 
measuring unmodulated output dc power, and a filter for measuring 
IM power centered between zero and -1.5 MHz. 

For determining power-transfer characteristics, the total output 
power, including all the IM as well as signal power, is the measured 
quantity. Figure 4 shows that the saturated output power becomes 
negligible outside of +100 MHz, indicating that terms higher than 
fifth order are negligible in describing the nonlinearity. In any case, 
we are safe in assuming that a filter bandwidth of, say, 200 MHz will 
be sufficient for measuring total output power. This procedure was 
followed in the computer calculations and in the laboratory measure- 
ments in which the bandwidth of the detector used to measure total 
output power was nominally 300 MHz. 


IV. COMPARISONS OF MEASUREMENTS AND CALCULATIONS 


In plotting the laboratory data, the procedure adopted was to nor- 
malize the measurement by dividing the average measured IM power 
by the noise bandwidth of the filter, thereby obtaining IM power per 
1-Hz frequency interval. Denoting this as J and the output unmodu- 
lated carrier power as C, the results were plotted as C/I, in dB-Hz, vs. 
output power, in dBm. In practice, the quantity of interest is the ratio 
of signal power to IM power in any given band, S/Njy. Since the 
signal and IM spectra are essentially flat, this is equal to the ratio of 
the signal and IM spectral densities, and for a large number of equal- 
power carriers equally spaced by AHz we have: 


S C 


Nac TA” (7) 


For these data, eq. (7) can be used with A set equal to the nominal 
intercarrier spacing of 1.5 MHz. 

Results of the computer calculations are also plotted in terms of 
C/I. In measuring IM power, an eight-pole Butterworth filter with 400 
kHz, 3-dB bandwidth was used; the measured noise bandwidth in this 
case was 420 kHz. The bandwidth was chosen to be as wide as possible, 
in order to allow for averaging, and also to be able to fit into the notch 
in between the unmodulated dc tone and the QPSK signal at —1.5 
MHz. In measuring the power in the dc tone at the output of the 
nonlinearity, an eight-pole Butterworth filter with a 20-kHz, 3-dB 
bandwidth was used, which was similar to that used in the laboratory 
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measurements. The following subsections discuss the results for the 
different devices studied. 


4.1 Single-stage SSPA, SS-1 


This single-stage SSPA used the same transistor as that used in the 
multistage MS-1 and MS-2 SSPAs. In this case, AM/PM character- 
istics for the particular amplifier that was studied were not available. 
A set of AM/PM characteristics typical of this type of single-stage 
amplifier is presented in Fig. 5, in which the phase shifts are seen to 
be very small. For this reason, as will be seen, unavailability of the 
exact set of AM/PM characteristics for this device had a negligible 
effect on the results. Calculations of C/J vs. output power are presented 
in Fig. 6 along with the measured values. For the three largest values 
of output power, which correspond to input saturation, 3- and 6-dB 
input back-off (~35-, ~34-, and ~32-dBm output power), the calcu- 
lated values differ from the measurements by the order of tenths of a 
dB. For the lowest point—13-dB input back-off, ~26-dBm output 
power—the measurement included significant thermal noise, which 
contributed to the difference of 2.5 dB from the calculated value. For 
the point at 10-dB input back-off (29-dBm output power), the differ- 
ence is 1.5 dB, which may also include some nonnegligible effects of 
thermal noise. 

In considering C/I performance, it is of interest to determine how 
much of the IM noise is due to AM/PM conversion, which is easily 
done in the simulation by setting AM/PM equal to zero. In this case 


4 


PHASE SHIFT IN DEGREES 


4 
1415 #16 #417 «18 #19 #20 21 «#22 23 #24 «#25 #26 «27 «#28 «#29 «430 (31 
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Fig. 5—Typical AM/PM characteristics for single-stage SSPA. 
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Fig. 6—C/I for single-stage SSPA. 


the AM/AM characteristics were modeled as those of the equivalent 
IEL, which is a very good approximation. The results for the IEL are 
also presented in Fig. 6 and seen to be only slightly different from 
those of the single-stage SSPA with AM/PM conversion, which dem- 
onstrates that phase shifts of the magnitude of those in Fig. 5, below 
saturation, are essentially negligible, and that the sharp increase in 
slope and larger phase values at and above saturation are also negli- 
gible in comparison with the effect of the saturating AM/AM charac- 
teristic. 

The AM/AM characteristics used in the calculation of C/I in Fig. 6 
are presented in Fig. 7, along with measurements and calculations of 
total power in vs. power out, and a plot of (6). The computer calcula- 
tions are seen to fall exactly on the curve of measured values, with the 
exception of the point at saturation where the difference is ~0.2 dB. 
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Fig. 7—AM/AM measurements and total power-transfer characteristics for single- 
stage SSPA. 


The IEL power-transfer characteristics of (6) are also seen to be 
within tenths of a decibel of the measurements and the computer 
calculations. 


4.2 Multistage SSPA, MS-2 


Figure 8 presents calculated and measured values of C/I. For input 
saturation, 3-dB input back-off and 6-dB input back-off (~36-, 
~35-, and ~34-dBm output power), the differences are fractions of a 
decibel. The largest difference is at 10-dB input back-off (~31-dBm 
output power), where the difference is ~1 dB. In this case the thermal 
noise in the measurements was negligible, as evidenced by the falloff 
in C/I below ~27-dBm output power. The computer calculations have 
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Fig. 8—C/I for MS-2. 


reproduced the irregular shape of the measured C/I curve in this 
region. 

The source of this irregular shape was determined to be the some- 
what irregular AM/PM characteristics for MS-2 shown in Fig. 9. As 
shown, the phase shift increases and then decreases somewhat in 
between —27- and —21-dBm input power. When the AM/PM was set 
equal to zero, the flattening in the C/I curve between 31 and 27 dBm 
output power disappeared, as shown by the equivalent JEL results in 
Fig. 8. 

We also note that although the phase shifts are largest and vary 
most rapidly near and above saturation (~—22-dBm input power, 
~37-dBm output power) (see Fig. 10), their effect on production of 
IM noise in Fig. 8 is of the order of 1 dB or less. However, although 
the phase shifts are much smaller and slower varying for low power 
levels, the differences between the C/I values for MS-2 with AM/PM 
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Fig. 9—AM/PM measurements for MS-2. 


and the equivalent IEL, with zero AM/PM, are much larger (e.g., ~6 
dB at ~27.5-dBm output power). The reasons for this, in terms of the 
relative importance of the AM/AM and AM/PM characteristics as a 
function of input power level, and the implications in terms of SSPA 
design and AM/PM compensation, have been discussed in Section I. 

The AM/AM characteristics for MS-2 are presented in Fig. 10, 
along with measurements and calculations of total power in vs. power 
out and a plot of (6). As for the single-stage SSPA, the differences 
between measurements, computer simulation results, and results based 
on (6) are negligible. It is of interest to note that the slope of the AM/ 
AM characteristics differs slightly from strict linearity in the gain- 
expansion region near saturation, which evidently has a negligible 
effect. 

For the remaining three multistage SSPAs that were studied, 
AM/PM measurements were not available and C/I calculations not 
possible. The comparisons were therefore restricted to total power in 
vs. power out. In these cases the laboratory measurements were 
compared with calculated values using the five-point approximation 
method discussed at the end of Section II. Differences between calcu- 
lated and measured values were found to be negligible; in the interest 
of brevity, plots of these results are not presented. For convenience 
we present in Fig. 11 the earlier result for MS-1 from Ref. 1, showing 
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Fig. 1O—AM/AM measurements and total power-transfer characteristics for MS-2. 


the significant discrepancy between calculated and measured power- 
transfer characteristics. The significant body of data in which almost 
exact agreement between measurements and calculations has been 
obtained clearly indicates that the MS-1 device represents an anom- 
alous case. 


V. SUMMARY AND CONCLUSIONS 


\ 

A number of SSPAs have been studied, for which very good agree- 
ment has been observed between calculated and measured wideband 
power-transfer characteristics and IM noise performance. The calcu- 
lations incude results of a computer simulation developed for this 
purpose, with which any arbitrary input signal and any set of operating 
conditions can be investigated, and an analytical result that describes 
the wideband power-transfer characteristics of typical SSPAs. In all 
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Fig. 11—AM/AM measurements and total power-transfer characteristics for MS-1, 
showing power-transfer anomaly. 


these calculations, the devices are characterized by their measured 
AM/AM and AM/PM characteristics; for power-transfer calculations 
AM/PM is, of course, irrelevant. 

The number of cases in which agreement has been obtained dem- 
onstrates the usefulness of calculations of this kind as means for 
predicting the performance of nonlinear power amplifiers, as well as 
for providing a tool with which possible anomalous behavior can be 
identified. Such behavior, in which the single-tone AM/AM and 
AM/PM measurements are not sufficient for characterizing the oper- 
ation of the amplifier under wideband excitation, which can result in 
severe performance degradation, can, in the case of power-transfer 
anomalies, be identified by comparing measured data with results 
predicted by (6). An example of anomalous behavior, in which a 
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significant degradation in power-transfer efficiency was identified by 
means of comparison with theoretical predictions, has been discussed. 

SSPAs in nonanomalous cases have been shown to generally exhibit 
power-transfer characteristics essentially identical to those of an ideal 
envelope limiter which, under the conjecture of Kaye, George, and 
Eric, are optimal. Part of Kaye, George, and Eric’s conjecture has been 
proved; namely, that AM/PM conversion can only degrade, never 
improve, IM performance. It has been shown that degradation due to 
AM/PM conversion is significantly more severe at low power levels 
than in the vicinity of saturation where, based on available AM/PM 
measurement data, the phase shifts are much larger and more rapidly 
varying; the reasons for this have been discussed. This could be of 
practical importance in the design of compensation networks for 
AM/PM conversion, since requirements on network performance at 
saturating power levels could be relaxed with little degradation in IM 
performance. 

Power-transfer characteristics are independent of AM/PM conver- 
sion. Ideally, this independence could permit the AM/PM character- 
istics to be adjusted or compensated to maximize IM performance 
without sacrificing ideal power-transfer efficiency. Using an analytical 
expression for power transfer in an IKL, a simple criterion has been 
developed for evaluating measurements to determine whether ideal 
power-transfer efficiency has been achieved. 
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Level Pictures 
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This paper presents a new variable-length coding algorithm for bit-rate 
reduction of digital images. The coding scheme is based on the fact that by 
taking into account the local conditional statistics of the picture, the coding 
efficiency can be improved considerably. Results of computer simulations are 
presented in terms of picture quality and required bit rate. Simulations show 
that for various Differential Pulse Code Modulation (DPCM) coded pictures, 
coding efficiency is such that a bit rate below the entropy of the DPCM 
quantizer output is attained. For a 3-bit DPCM encoded picture, a bit rate of 
1.64 bits/picture element is achieved. 


I]. INTRODUCTION 


Two popular methods for reducing bit rates of pictures are Differ- 
ential Pulse Code Modulation (DPCM)! and transform coding.? When 
these methods are applied, their reconstructed signals are not exactly 
the same as the original Pulse Code Modulation (PCM) input because 
of the quantization noise. Consequently, these encoders are often 
referred to as nonreversible encoders. Bit-rate reduction can also be 
achieved by reversible encoders. These methods are based on the fact 
that it is not necessary to assign an equal number of bits to quantized 
samples that are not equiprobable. Therefore, it may be advantageous 
to design variable-length code words in which the length of the code 


*AT&T Bell Laboratories; present affiliation Bell Communications Re- 
search, Inc. 





Copyright © 1984 AT&T. Photo reproduction for noncommercial use is permitted with- 
out payment of royalty provided that each reproduction is done without alteration and 
that the Journal reference and copyright notice are included on the first page. The title 
and abstract, but no other portions, of this paper may be copied or distributed royalty 
free by computer-based and other information-service systems without further permis- 
sion. Permission to reproduce or republish any other portion of this paper must be 
obtained from the Editor. 


249 


words are matched to the Probability Distribution Function (PDF) of 
the source by using very simple procedures. One such procedure has 
been devised by Shannon? and Fano and another by Huffman.’ An- 
other reversible coding method is run-length coding, in which a se- 
quence of identical symbols is replaced by a code indicating the length 
and symbol of the sequence. Run-length codes are useful whenever 
the long runs are highly probable. This method has been widely used 
for compression of bilevel pictures.° 

With the above-mentioned reversible coding methods, the entropy 
obtained from the PDF of a digital image gives a target for the 
minimum bit rate necessary to describe the source, ignoring the 
conditional statistics of the picture. If, by some efficient scheme, 
conditional statistics can be considered, it may be possible to lower 
this target to a figure that can be described by conditional entropy. 
One such method has recently been devised by Gharavi et al. for 
encoding the gray-level pictures using CCITT one-dimensional run- 
length codes.® In this scheme the run-length statistics of the picture 
are changed in such a way as to provide a signal that has long bit runs. 

In this paper we introduce a similar strategy by developing a scheme 
in which the picture can be directly coded by the variable-length codes. 
The code words are designed according to the local conditional statis- 
tics of the picture using the Huffman procedure. For example, the 
picture element (pel) values that are expected to occur most frequently, 
given the local PDF, are assigned to the shortest code words in 
accordance. The mathematical model and the details of the coding 
algorithm are given in the following section. 


Il. BASIC MODEL 


The n-bit PCM coded picture can be described mathematically as a 
random sequence of 2” possible outcomes with certain statistical 
properties, which are given in terms of probabilities. From these 
probabilities we can calculate the entropy of one pel, which represents 
a lowest theoretical bound for the one-dimensional memoryless model. 
The entropy of pictures can be presented as: 


Qn 


Fpe1 pa > P;LogeP;, 
i=1 


where P; is the probability of the ith outcome (ith level). In the case 
of an mth-order Markov model, the luminance level of each pel, x;, is 
considered to be dependent on the level of m previous surrounding 
pels. Therefore, the conditional entropy based on the m previous 
elements is given by: 
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where P(xj-mXi-m-1 ‘++ Xj) is the joint probability of x; --- x; and 
P(xi/Xj-1, Xi-1 +++ Xi-m) 18 the conditional probability of x;, given m 
previous pels x;-1, X1-2 --- X;-m. Equation (1) can be rewritten as 
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and, on defining Hy as a local entropy, 


gn 


Ay —-> P (xi/ Xi-1, Xi-2 ++: Xi-m)LogsP (x;/xi-1, eis <n) (3) 


We express the conditional entropy as 
H, = YY +++ Y Pim, Xiem-1 +++ Xi-1) AL. (4) 


Equation (4) shows that the conditional entropy depends not only on 
the local entropy but also on its distribution over (2”)” independent 
states. This model leads to the new coding scheme proposed here. 


2.1 Coding algorithm 


The first step is to use the m previous pels to construct a table of 
(2”)™ independent states of picture statistics. For each state, a set of 
code words is designed according to the probability distribution of the 
2” outcomes of the quantized values. These code words can be derived 
by employing optimum coding procedures, such as Huffman, or Shan- 
non and Fano, in which a shorter code is assigned to a more probable 
level. Then, for a given state, each incoming pel is coded by selecting 
an appropriate code word corresponding to its quantized value. 


Ill. DPCM CODING 


In the previous section it is apparent that the total number of states 
is heavily dependent on the number of quantization levels. For exam- 
ple, for input pictures of 256 levels corresponding to 8 bits/pel, there 
are 256” (m is the number of reference elements) states. This figure 
is rather high, and a simple reduction of the number of quantizing 
levels would lead to annoying contouring effects. Thus, to convert an 
8-bit picture into fewer bits while maintaining approximately the same 
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quantizer noise, a DPCM encoder is employed. In the DPCM encoder 
used here, a linear predictor constructs its predicted pel as the weighted 
summation of previous pels (see Fig. 1); thus 


x= 0.5A + 0.25(B + C), (5) 


where x is the prediction of the present pel. The prediction error is 
then quantized by a symmetric seven-level quantizer with the transfer 
characteristics given in Fig. 2. Applying the above coder to the 512 x 
512, 8-bit/pel picture shown in Fig. 3 results in the seven-level DPCM 
coded picture shown in Fig. 4. 


3.1 Subsampling and interpolation 


Subsampling at a factor of 2:1 gives approximately a 2:1 bit-rate 
reduction. This results in some degradation of pictures, which can be 
minimized by an efficient interpolation scheme. Figure 5 shows the 
2:1 horizontal subsampling pattern employed here. The interpolation 
is based on averaging the luminance levels of the four surrounding 
pels. The subsampled picture is DPCM encoded, in which the present 
pel D (see Fig. 5) is predicted by 


Prediction of D = 0.3B + (C + F)0.35. 
Figure 6 shows the interpolated picture where 2:1 subsampling was 
employed. 


IV. COMPUTER SIMULATIONS 


Experiments were carried out by means of computer simulation. 
Figure 7 shows the transmitter and receiver block diagrams. In our 
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Fig. 1—Configuration of pels used for prediction. 
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Fig. 2—Transfer characteristics of the quantizer. 
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Fig. 3—Original picture. 


simulation, for a set of seven-level DPCM coded pictures, tables of 49 
and 343 states were constructed by observing the quantized values of 
two and three previous elements, respectively (A, B and A, B, C of Fig. 
8). As a matter of convenience only, the table of 49 states with the 
corresponding code words is shown in Table I. For each state, the 
seven code words are obtained from the probability distribution func- 
tion of the seven outcomes, using the standard Huffman procedure. 
As Fig. 7a shows, the input signal is first DPCM encoded to generate 
a 3-bit DPCM coded picture. The coded picture is then applied to the 
conditional variable coder, which uses the line and sample delay to 
represent the values of neighbors B and A, respectively. These values, 
along with current pel value, comprise the input for the encoder Look 
Up Table (LUT). The encoder LUT uses this information to give the 
bits of the appropriate variable-length code word, as shown in Table 
I. The output is smoothed by a buffer to allow for the uniform bit rate 
in the transmission channel. In the receiver, the transmitted variable- 
length code words are received from the channel and applied to a 
buffer and then to the decoder LUT (see Fig. 7b). The conditional 
variable-length decoder LUT also receives the decoded sample and 
line-delayed signals to represent the values of neighbors B and A, 
respectively. Based on this information, the decoder determines 
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Fig. 4—Three-bit DPCM coded picture. 


Xx O xe O xe 
O x8 cA x? o) 
x O xF 0 x 
O x fe) X fe) 


X : SAMPLE SELECTED FOR TRANSMISSION 
O : SAMPLE DROPPED 


INTERPOLATION OF A = (B+C+D+E) / 4 


Fig. 5—Subsampling pattern used for interpolation and prediction. 


whether the current pel code word is the one expected to occur most 
frequently, second most frequently, and so on, as specified in Table I. 
The output of these is then applied to the DPCM decoder, which 
serves as the receiver output signal. 
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Fig. 6—Subsampled and interpolated picture. 


4.1 Results 


Table II presents the bit-rate results of the variable-length encoder 
with the corresponding entropies and conditional entropies. The re- 
sults of variable-length coding based on the statistics of the quantizer 
output (ignoring the conditional statistics) are also included in Table 
II as a reference. 

By looking at the tabulated results, it is clear that a considerable 
improvement is obtained by using this conditional variable-length 
algorithm compared with a nonconditional variable-length encoder. 
An even better improvement is obtained by increasing the number of 
reference elements from two to three (third row of Table II) at the 
expense of complexity, as the number of states is increased from 49 to 
343. In the case of 2:1 subsampling, a 2:1 bit-rate reduction is unob- 
tainable. This is due to the lower spatial correlation as a result of an 
increased distance between neighboring pels, as compared to a non- 
subsampled case. However, the most important conclusion from Table 
II is that bit rates below the entropy of the DPCM quantizer output 
are achievable. 


V. CONCLUSION 


A new concept for a variable-length coding scheme has been pre- 
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Fig. 7—Block diagram of (a) transmitter and (b) receiver. 
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A x 


SSS SS SENT LINE 


~ 
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Fig. 8—Configuration of reference elements used for conditional variable-length 
coding. 


Table |—Table of variable-length code words 


Neighbors 
State A B Code Words 
1 lL lh L; Lz Ls Ly L; Le Ls 


1 01 001 00001 000001 000001 000000 
L L L 


2 L L Ls 5 7 6 

1 le 9 ii 101. 1001 10001 100001 100000 
3 L L »w b Le Ls Ly Ig Ly 

1 ls 9 11 101 1001 10001 100001 100000 
th 2 b Li L 
5 


3 L, 5 ; Leg 
11 10 01 001 0001 00001 00000 
L L Ls; Ls L, L, Lg 4 L; 
A By a 10 01 001 0001 00001 00000 
6 L Le ls Is Ls Ly Ly Ly Ly 
4 : 1 01 0011 0010 0000 00011 00010 
7 LL ob lh Ls Lg Ls Ls La 
‘ ‘ 1 01 0011 0010 0000 00011 00010 
3 LL » lb kL Ly Ls 
0 11 101 1001 10001 100001 100000 
9 L L Li L, L; L, L; 7 6 
‘i 0 11 101 1001 10001 100001 100000 
10 L L Ls L; L, Ls Ly Lg L, 
: 0 11 101 1001 10001 100001 100000 
11 fe oe Lp Ly Li L3 L, L; Le 
? cme 10 01 001 0001 00001 00000 
12 L L Ls L; Li L. Le ly L, 
ss e 10 01 001 0001 00001 00000 
13 Le Le we ls Le Ls Ly Ly Le 
Q . 0 01 0011 0010 0000 0011 00010 
14 L L L; Ly Le Li L, L3 L; 
i 0 01 0011 0010 0000 0011 00010 
5b L L mM le Ls Ly Ls Lg Ly 
: 0 ll 101 1001 10001 100001 100000 
6 lL LD mM lb kL Ly Ls Ly Le 
; 2 0 11 101 1001 10001 100001 100000 
17 Le L Li Ls L, L; Ly Le Ly 
: : 0 11 101 1001 10001 100001 100000 
1% LL L Lp Li Lg L3 L, Ls Le 
3 0 11 101 1001 10001 100001 100000 
19 L. L L; Li Ls L, Le Ly L; 
3 5 11 10 01 001 0001 00001 00000 
20 L L Le L; Ls L, Li L, Ly 
. s 0 01 0011 0010 0000 00011 00010 
21 L L Li Le la L3 Lz L; Li 
: c 0 111 110 100 1010 10111 10110 
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State 


22 
23 
24 
25 


26 
27 
28 
29 
30 


31 
32 
33 
34 
35 


36 
37 
38 
39 
40 


41 
42 
43 
44 
45 


46 
47 
é8 
49 
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Neighbors 
A BS 
le li 
Ly Le 
L, Ls 
Ile Ly 
L, Ls 
Li Le 
L, Iy 
L; Ly 
Ls; Ly 
Ls Ls 
Ls Ls 
L; Ls 
L; Le 
L; Lr 
Le li 
Le IL: 
Ls Ls 
Le Ls 
Le Ls 
Le Le 
Le Lr 
L; la 
lL, I, 
L, Ls 
L, Ly 
L; Ls 
L; Le 
L, Ly 


Table |I—Continued 


L, 


Code Words 
Ls Lg 
1001 10001 
Ls Ls 
001 0001 
Le ly 
1001 10001 
L; L3 
1001 10001 
Ls Le 
100 1010 
Ll, Le 
001 0001 
Ls Le 
1001 10001 
L; Ly 
001 0001 
Ls Ls 
100 1010 
Ls Ly 
001 0001 
Li Ls 
001 0001 
Li Le 
001 0001 
Ly Le 
1001 10001 
Ls Lp 
1001 10001 
Ls Le 
010 000 
Ly Ls 
1011 1010 
Ls Ls 
1111 1110 
Le Li 
001 0001 
Li Le 
001 0001 
L7 Li 
0010 0000 
L, Lz 
1001 10001 
La Le 
1111 1110 
Ll Li 
010 000 
ly L3 
100 1010 
Ls li 
001 0001 
Ls Ls 
010 000 
Ls Ls 
00111 00110 
li L3 
00111 00110 


1 


Ly 
100001 
L 


7 
00001 
Ly 
100001 
Ls 
100001 


Ly 
10111 
L, 
00001 
Ii 
100001 


Lg 
00001 


1 
100001 
L; 
0011 
Ly 
1001 
L, 
1001 
Ls 
00001 
Ly 
00001 


Ly 
00011 
Ls 
100001 
L: 
1001 
Ls 
0011 
Ly 
10111 


Le 
00001 


L, 
0011 
L2 
00101 
Le 
00101 


Le 
100000 


Le 
00000 
Le 
100000 
Lg 
100000 
Li 
10110 
L; 
00000 
L; 
100000 
L, 
00000 
L; 
10110 
L, 
00000 


4 
100000 
Ls 
100000 
Ly 
0010 
Ls 
1000 
ly 
1000 
L; 
00000 
L; 
00000 
L, 
00010 
L 
100000 
Ly 
1000 
Le 
0010 


L, 
10110 


L; 
00000 


Ly 
0010 
Li 
00100 
L; 
00100 


Table II—Performance of variable-length coding algorithms for 512 x 
512 pictures 
Entropy 
of Quan-  -Condi- 


tizer Out- _ tional 
put (bits/ Entropy Bit Rate 


Variable-Length Encoder DPCM Encoder pel) (bits/pel) (bits/pel) 
Based on quantizer output Nonsubsampled 1.84 — 1.94 
Based on conditional statistics | Nonsubsampled 1.84 1.64 1.73 
ref. elements A, B 

Based on conditional statistics | Nonsubsampled 1.84 1.53 1.64 
ref. elements A, B, C 

Based on conditional statistics | Subsampled 1.05 0.92 1.01 
ref. elements A, B 

Based on conditional statistics | Subsampled 1.05 0.89 0.97 


ref. elements A, B, C 


sented in which the structure of the variable-length code word changes 
from pel to pel, depending on the local conditional statistics of the 
picture. In this scheme 3-bit DPCM words are coded by assigning a 
variable-length code word to represent each pel as a function of its 
expected frequency of occurrence, given the neighbor pel values (state). 
We have shown that for subsampled and nonsubsampled 3-bit DPCM 
coded pictures, bit rates of 0.97 bits/pel and 1.64 bits/pel, respectively, 
are achievable. We therefore conclude that by using the conditional 
variable-length encoder it is possible to obtain bit rates below the 
entropy of the quantizer output. 
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Signals Designed for Recovery After Clipping—I. 
Localization of Infinite Products 


By B. F. LOGAN, JR.* 
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It is shown for signals of the form s(t) = g(t) + cos ct, where g is band- 
limited to [—6, b],0<b<c< ©, and (—1)*s(kx/c) > 0 (e.g., if | g(t) | < 1), that 
the zeros of s in the interval (t — T, t + T) allow reconstruction of s(t) witha 
relative error less than 2e~*7(1 — e~*7)-, where \ = b — c. The reconstruction 
algorithm gives an effective way of localizing the infinite product representa- 
tion for s(t). Some practical aspects of the algorithm are discussed. The scheme 
offers advantages for the recovery of signals after nonlinear distortion, such 
as occurs, for example, in magnetic recording. 


l. INTRODUCTION 


Let us assume that a signal s(t) undergoes some nonlinear distortion, 
generally of unknown characteristic, during transmission (or perhaps 
in recording), which does not alter the zero crossings. We may suppose 
that the reliable information we receive is just sgn s(t), the clipped 
signal. We are led to ask for what class of signals s(t) can we make a 
reliable estimate of As(t) from sgn s(¢), i.e., from the zero crossings of 
s(t). 

If we assume that s(t) is a (real-valued) bandlimited signal with real 
simple zeros t,, we have the formula,’ assuming s(0) # 0, 


* AT&T Bell Laboratories. 

t This follows from a result of Titchmarsh,! who showed that (1) holds for s(t) that 
have ordinary Fourier transforms that vanish outside an interval [~a, a] but not out- 
side any translate of the interval; i.e., for bandlimited functions whose spectral end- 
points are centered about the origin. Then, (1) follows for bounded s(t) by considering 
s(t)(t — t;)""(t — t,)~', where t; and ¢, are either real or complex conjugates. 
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S50) 0 ( = t) (1) 


tp)’ 
where the product converges (conditionally) if the zeros are ordered 
such that | trai | > | t,|. Thus, from a strictly mathematical viewpoint, 
we have a “formula” for recovering s(t) within a constant multiplier 
from sgn s(t). However, for practical purposes the formula is worthless 
without additional restrictions on s(t). 

In the first place, the formula is valid only for strictly bandlimited 
s(t). If the t, are the zeros of an almost bandlimited s(t), then the 
product may not even converge, or if it does, it may not even represent 
a bounded function, and in general cannot be expected to give a 
reasonable replica of s(t). 

In the second place, even if s(t) is strictly bandlimited, the infinite 
product representation is, in general, impractical; it is very sensitive 
to (correlated) perturbations of the t, and, more important, knowledge 
of all the t, is crucial to estimating s(t), relative to s(0), for large t. In 
practice, one would like to make a good estimate, s7(t), of s(t) just 
from the knowledge of the zeros of s in the interval (t — T, t + T) for 
some fixed T and arbitrary t, (-%2 < t < ). That is, in practice we 
require a localized reconstruction formula. 

For a certain class of bandlimited s(t), we can obtain an exponentially 
localized reconstruction formula; i.e., from a knowledge of the 
zeros of s in the interval (t — T, t + T), we can make an estimate s7(t) 
such that 


| s(t) — sr(t)| S | s(t) |-er, (2) 
where 
eT = Ofe*7} T— © 


and \ > 0 is an appropriate parameter in the representation of s(t). 
This class of signals, which we denote by S(6, c), consists of signals 
of the form 


s(t) = g(t) + cos ct, (3) 


where g (real valued) is a signal whose spectrum is confined to [—, 6], 
0<b<c< ©, and such that 


(—1)*s(kx/c) > 0, = 0, +1, +2,---. (3a) 


The alternation condition (3a) ensures that s(t) has only real simple 
zeros,* one between each extremal of cos ct. Sufficient additional 


* We ask the readers to accept this fact, or provide their own proof. A more general 
result will be given in a future paper. 
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conditions on g for (3a) are 
| g(kx/c)| <1 k = 0, +1, +2, --- (3b) 
and 
}e(t)| <1 —-~<t<o, (8c) 


the latter condition being the more practical constraint. However, (3a) 
allows special functions g of the form 


g(t) = (1 + x(t))cos ct + x(t)sin ct, (4) 


where x(t) is a bandpass signal (real valued) whose spectrum is 
confined to the disjoint intervals, [\, w] and [—p, —A],O<A <u, 


A+ us 2c, (4a) 
x(t) > —1, (4b) 

and 
x(t) is the Hilbert transform of x(t). (4c) 


The spectrum of g, defined in (4), is confined to the union of the 
two intervals, [—c + A, —c + yw] and [—p + c, —A + c]. In case c > yp, 
then g, defined in (4), is bandpass with spectrum confined to the 
disjoint intervals [—b, —a], and [a, b], where 


0O<a=c-yw<b=c-.. (4d) 


In this case (c > u), we may consider s(t) to be a full carrier lower- 
sideband signal associated with the bandpass baseband signal x(t). 
Also, if g in (3) is bandpass with spectrum confined to the disjoint 
intervals [—b, —a], [a, b], and 0 <a<b<c, we can always put g in the 
form (4). 

The important parameter associated with the class S(0, c) is 


=6 =D; (5) 


In words, ) is the width of the guard band between the top frequency 
of g(t) and the frequency of cos ct. We will refer to \ as the gap 
frequency (rad/s). It is the presence of the gap that allows the expo- 
nential localization of the reconstruction of s(t), giving an error of the 
order of e~*? from the knowledge of the zeros of s in the interval (t — 
T, t+ T). We can obtain a good approximation to s(t) for AT, say, in 
the range 27 to 37; i.e., T needs to be only 2 or 3 half periods (Nyquist 
intervals) of the gap frequency to obtain 0.1-percent error. 

We should note that the signals in S(0, c) are, in effect, normalized 
by taking the coefficient of cos ct to be +1. With the signals thus 
normalized, their zeros give a complete description of the signals. 

In some applications we may think of s as a baseband signal g with 
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added bias (cos ct), and in other applications we may think of s as a 
lower-sideband signal with full carrier (cos ct). 

We first present the reconstruction formulas. Subsequent sections 
are devoted to the derivation of the formulas. Finally, some practical 
aspects of the reconstruction formulas are discussed in the concluding 
section. 


Il. THE RECONSTRUCTION FORMULAS 


Here we describe an approximation, s7(t), determined only by the 
zeros of s in the interval (¢ — T, t + T) such that for s in S(b, c), we 
have 





Oe! 
| s(t) — sr(t) | =q7_ 27 [s(t)|,  -~<t<o, (6) 
The first form of the approximation is 
sp(t) = (1/2){sgn s(t)}-exp{hr(t)}, (7) 
where 
fin(t) = Lnlt ~ ) + £ wT), (7a) 
Ly(t) = Lr(—t), (7b) 
Lr(t) = 0, \t| > qT; (7c) 
Lr(t) = - po dx, 0<t<T, (7d) 
T sinh AVT? — t? 
MO = Soh PTR” Ve 
and 








_ x1o(x) e 2 { cee \/% 
AVS sinhx sinhx Jo © o i a) 


(Here Jp is the modified Bessel function.) The function L(t) is a 
sort of truncated log |t/T'|, depending both on X and T, although in 
the notation we have suppressed the dependence on i. Using the fact 
that (0) = 1, we have 














T 
t 1 — f(x) 
= —_— d 
Lr(t) = log T + f ay 
t 
= re pas <T 
log T +F (5: x7" |t| 
= 0, |\¢| = T, (8) 
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where F(x; \T) is analytic everywhere. Graphs of F(x; 8) and L7(xT), 
8 = XT = kr/2 are shown in Fig. 1 for -1 <x <1, k= 1, 2,3, --- , 6. 
For \T = z/2, in which case the reconstruction interval (t — T, t + T) 
has a length equal to only a half-period (one Nyquist interval) of the 
gap frequency, the function F is not significantly different from zero 
over the interval (—1, 1); so L7(t) = log |t/T'| for the case AT = 7/2. 
For larger XT, the function F becomes significant. 
In order to express s7(t) in product form we define 


Vr(t) = exp Lr(é). (9) 


Here, V(t) is a sort of V-shaped function that we may write, using 


(7), 





Fig. 1—The functions Lr(xT) and F(x; 8) for 8 = kr/2, k = 1, 2,3, ---,6. 
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t t 
Vr(t) = | w(f: Ma) jt| << T 


= 4, |¢| = fie (9a) 

where 
W(x; AT) = exp{F(x; AT)}, Ea ee (9b) 

We have then, 

sr(t) = Af{sgn s(t)} Tl Vr(t — t), (10) 

where 

1 Cc 

A= 9 exp \s ar} . (10a) 


Graphs of the two functions, Vr(xT) = | x|W(x; AT) and W(x; AT), 
are shown in Fig. 2 for -1 < x < 1, AT = B = kw/2, and k = 1, 2, 3, 

., 6. 

The normalization function u(AT), defined in (7f) and appearing in 
(7a) and (10a), is graphed in Fig. 3. 

It should be noted in formulas (7) and (10), recalling that L7(t) = 0 
for |t| = T and Vr(t) = 1 for |t| = T, that the infinite sum and the 
infinite product involve for each t only those t, that satisfy |t — t,| < 
T; i.e., only those t, in the interval (t — T, t + T). So, (10) gives the 
desired localization of the infinite product (1) for signals in the class 
S(b, c). 

The error estimate in (6) is conservative, with equality possible only 
at the zeros of s. At least, we have the exponential decrease of the 
error with AT; e.g., for AT = 27, 


9 24 
a = 0.00375, 
and for AT = 37, 
y) —3r 
qos = 0.00016. 


So, for T somewhere between 2 and 3 half periods of the gap frequency 
d, we get a relative error, (sr(t) — s(t)) + s(t), amounting to no more 
than 0.1 percent. 

Figure 4 illustrates the simplest sort of test of the algorithm in 
reconstructing s(t) = cos t. Here g(t) = 0, and \ = c = 1. Four cases of 
sy(t) are shown corresponding to T' = 7/2, 7, 3/2, and 2x7. The 
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1.0 





Fig. 2—The functions Vr(xT) and W(x; 8) for 8 = kn/2, k = 1, 2,3, --- ,6. 


corresponding errors and relative errors, suitably magnified, are also 
shown. 

In the first case (T' = 2/2), shown in Fig. 4a, only one zero of cos t 
appears in the observation interval (t — T, t + T); so, we have in this 
case 


sr(t) = A-Vr(t — m/2) for O<sts< m/2. 


Here Vr(t) = Vr(t; AT) with T = 7/2, X = 1. For the relatively small 
product, \T = 7/2, the function V7(t) is essentially | t/T| for |t| < T 
(cf. Fig. 2), so that s7(t) has a quite triangular waveform. The magni- 
tude of the error (sr(t) — s(t)) in the reconstruction is maximum at 
t = 0, where it is approximately 0.22. The magnitude of the relative 
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2.5 


2.0 


por) 


0.5 


a 


Fig. 3—Normalization function p(AT). 


error, (s7(t) — s(t)) + s(t), is also maximum at ¢t = 0 with the same 
value, approximately 0.22. The estimate (6) gives an upper bound for 
the absolute relative error in this case of approximately 0.66. The 
reconstruction is rather crude in this case. 

The approximations become progressively better in the other cases. 
For example, the maximum absolute relative error in the fourth case, 
AT = 2x (see Fig. 4d), is approximately 0.001, the upper bound (6) for 
this case being approximately 0.0037. In all cases the maximum 
absolute relative error is approximately one third of the upper bound, 
given by (6). Note that in the second and fourth cases the absolute 
relative error is maximum at the zeros of cos t where the error is zero. 
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Fig. 4(a) and (b)—Reconstruction of cosine. 
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Fig. 4(c) and (d)—Reconstruction of cosine. 
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Note the similar character of the error waveforms in the first and 
third cases, where 2T is an odd multiple of z, which is distinctly 
different from that of the second and fourth cases, where 27' is an 
even multiple of 7. 


Ili. THE BASIC FORMULA 
Given the zeros {t,} of a signal s(t) in S(b, c), we define the 
following important function: 
h(t) = J(t) — ct, (11) 


where J(t) is a jump function increasing by z at each zero ¢;, of s(t), 
and 


J(0) = 0. (11a) 


The function z7'J(t) is just the zero-counting function of s(t). Since 
there is one, and only one, zero t, between each extremal of cos ct, 
we may label the zeros such that 


kr/e < tp < (k + 1)x/c, k= 0, +1, +2,.---. (12) 
It is clear then that 
h(kr/c) = 0, k = 0, +1, +2, --- (13) 
and 
—r<h(t)<7, -2w<t<», (14) 


The function h(t), which we will refer to as the fundamental function 
associated with the zeros {t,} of s(t), may be described as a sort of 
meandering sawtooth function bounded between —z and z. 

In this section we want to show that h(t) is a high-pass function 
“with no spectrum in the frequency interval (—A, \), where X is the 
gap frequency (A = c — b), and therefore has a Hilbert transform 
h(t). Furthermore, we will show that 


log] 2s(t) | = A(t) = * ff nf) 


-o t — 





dx, (15) 


and thus obtain the basic formula 
s(t) = (1/2){sgn s(t)jexp{A(d)}. (16) 


The fact that h(t) is a high-pass function allows us to make a 
good estimate h(t) of A(t) just from the knowledge of h, or equiva- 
lently the zeros of s, in the interval (¢ — T, t + T). Later, we detail 
the localization of the Hilbert transform to obtain (7) from (16). 

To obtain these results, we introduce the complex variable 7 = 
t + iu, and consider for s(t) in S(b, c) the function 
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H(r) = i log{2e*"s(r)},  u>0. (17) 
We have 
2e""s(r) = 2e”"( g(r) + cos cr) 
= 1 + 2e"g(r) + e7*7, (18) 


Now g(r) is an entire function of exponential type b < c, real and 
bounded on the real line; so (cf. Ref. 2, p. 83), 


| g(t + iu)| < (cosh bu)sup | g(é) |. (19) 


Hence, 
Qee(iMs(t + iu) =1+0(e), u>, 
where AX=c-—b>0, (20) 
and therefore, 
H(t + iu) = 0(e~”), U —> 00, (21) 


Also, since s(r) has only real zeros, H(7) is analytic in the upper 
half plane. The real part of H(r) is also bounded. It follows (Ref. 3, 
see Theorems 6.3.1 and 6.4.3) that 


h(t) = lim Re{H(t + iu)} (22) 
u-0+ 
and 
A(t) = lim Im{H(t + iu)} (23) 
u-0+ 


are high-pass functions having no spectrum in the interval (—), A), 
and h(t) is the Hilbert transform of h(¢). 


We have 
h(t) = — lim arg{e“"(g(r) + cos cr)}. (24) 
u-0+ 


Now, for t = ka/c we have 
e*( g(t) + cos ct) = e”*s(kx/c) > 0. (25) 


Hence, h(kz/c) is some multiple of 27. We can conclude that h(kz/c) 
= 0 from the fact that h is high pass and from the fact that the 
argument (phase) of s(t) is a decreasing function of t, decreasing by + 
at each zero of s (indent into the upper half plane around the zeros). 
There is one, and only one, zero between each extremal of cos ct, and 
h’(t) = —c except when t is a zero of s(t). Thus, 


h(t) = h(O0) — ct + S(t), (26) 
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where J(t) is a nondecreasing staircase function (jump function), 
J(0) = 0, which increases by = at t;, the zeros of s(t); i.e., 


kr (k + 1)x 
— <i << ——— 


; (26a) 
J’(t) = 0, t # tr, (26b) 
and 
J (=) = kr. (26c) 
Thus, we have from (26) and (26c) 
h (*=) = h(0) = 2nz, k=0, +1, +2,---. (27) 


Since h(t) can increase or decrease at most by z in each interval 
kr/eo<t<(k+1)z/c, we have 


—a < h(t) — h(kr/c) <a 
or 
(2n — 1)xr < h(t) < (2n + 1)z. 


But, h(t) is high pass and must change sign, so we must have n = 0, 
and hence, 


—-x1<h(t)<a (—o <t < 0), (28) 


In other words, h(0) = 0 is the only value of arg{e“‘s(t)} at t = 0 
consistent with H(t + iu) = 0(e~™), where we took log{2e"’s(z)} to be 
the principal branch, log{1 + 2e”7g(r)} — 2e""g(r), u — ». Thus, we 
have from (17), (23), and (28) 


log 2| s(t) | = : ia ae dx, |h(x)| <7, (29) 


where h is the function defined in (11). Now, since h is a bounded 
high-pass function, we can obtain an exponential localization of the 
Hilbert transform in (29). 


IV. LOCALIZATION OF THE HILBERT TRANSFORM 


If the spectrum of a bounded function h(t) omits the interval 
(—A, A), there are many equivalent representations for its Hilbert 
transform. We have* 


h(t) = if h(x)K(t — x)dx, (30) 
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where* 


K(t) = eo) (30a) 
at 
and 
f(t) is any bandlimited function whose spectrum is 

confined to [—A, A], satisfying (30b) 
f(0) =1 (30c) 

and 
{ IO] dt < ©, (30d) 

jei>1 | t| 


Although A(t) is not generally bounded, the singular kernel K(t) can 
be truncated to an interval (—T, T) giving K7(t), such that 


es T 
h(t) = h(x)Kr(t — x)dx = i” Kr(x)h(t — x)dx (31) 


is a good approximation to A(t); i.e., 


h(t) — hr(t) = if ix) h(t — x)dx. (32) 
|x|>T WX 
Then 
|A(t) — hr(t)| < M ii Lis) | dx, (33) 
Ixj>7T 7|x| 
where 
M = sup |hA(t) |. 


The bandlimited function f that minimizes the integral in (33), with 
f(0) = 1 and spectrum confined to [—A, A], is (see Ref. 4) 


T ~~ sin A\vt? — T? 








= pee 34 
KM) = aha? Ve — a 
or, equivalently, 
T sinh AVT? - ¢t? 
f(t) = an" aera eee) 


sinh AT VT? — t? 


* Formally, the Fourier transform of K is the convolution of F(w), the Fourier 
transform of /, with —i sgn w, which then agrees with —i sgn w for |w| >A. 


274 TECHNICAL JOURNAL, FEBRUARY 1984 


So, we will set 

1 T__ sinh AVT? = ¢? 

at sinhdT VT? —- 2 ’ 
=0.. ear (35) 

Then, (34) in (33) gives (see Ref. 4) 


‘ 2M Lake 
[ A(t) — hp(t)| < — log ———> 


T 1—e >?" 


Kr(t) = 





lt] <T 


(36) 


So Kz, defined in (35), gives us the desired exponential localization 
of the Hilbert transform. 


V. LOCALIZATION OF THE RECONSTRUCTION OF s(t) 


Now we are prepared to obtain a good approximation to s(t) belong- 
ing to S(b, c), c — 6 = A, from a knowledge of the zeros of s in 
(t — T,t + T). From (29) we have 


s(t) = (1/2){sgn s(t)}exp A(t), (37) 


where / is the Hilbert transform of the high-pass function h defined 
in (11). Now we set 


sr(t) = (1/2){sgn s(t)jexp hr(d), (38) 
where 
hr(t) = { ; h(x)Kr(t — x)dx (39) 
and Kris defined in (35). Since | h(x) | < 7, we have from (36) 
| A(t) — hr(t)| < 2 log ite. (40) 


Owing to the special nature of h(x) (cf. [11]), we can manipulate the 
integral (39) into a more convenient form. We may write 


hr(t) = Z : dh(x)Lr(t — x)dx 


Tv 


= y Ln(t — t) ~~ { _ Exod, oe 


where {t,} are the zeros of s(t) and 


Lr(t)=7 ii K7(x)dx, —o <t<Q, (42) 
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Lr(t) = Lr(-2), (42a) 
Lr(t)=0, |t| >T, (42b) 
and 


* f(x) 
Lr(t) = - a dx, O<t<T. (42c) 
t 
The function L(t) is a sort of truncated log | t/T'|, depending both on 
d and T, although in the notation we have suppressed the dependence 
on \. Using the fact that f(0) = 1, we have 








eae a + | ee lt) <T 
=0, |t/=T (42d) 
and furthermore, since f(x) > 0 for |x| <= T (ef. [42c]), 
Lr(t) < 0, —-~<t<o, (42e) 


Although it may be simpler to compute hr(t) from (41) and then 
compute s7(t) according to (38), it is interesting to express s7(t) in 
product form. So we define 


V7(t) = exp Lr(¢). (43) 


Here, V7(t) is a sort of V-shaped function that we may write using 








(42d), 
t "1 = f(x) 
V(t) = |= ex | | 1-10) art, |t]| <T 
T t x 
ae |t| = T, (43a) 
and we note that 
T 
V7(0+) = = exp { eae iC) dx (44) 
T 0 x 
and 
i) eee 
V(T-) = T smh TT" (44a) 


Now we may write 
s(t) = A{sgn s(é)} i Vr(t — th), (45) 


where, from (41) and (38), 
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A=A(b,c, T) = 5 exp |-2 { Ln(oatt (45a) 


We have 


[trae = = [7 tL 7(t)dt = - | Hoa (46) 


From a table of Fourier transforms,’ we find 


a, a 
a cos wtdt = rIp(TVd? — w”), lw| <A 
="), Jwl >A (47) 
and 
T 
h \vt? — T? 
es aa cos wtdt = ad o(T Vw _- x*); (48) 


where I(x) = Jo(ix) is the modified Bessel function. Setting w = 0, we 
have 


* sin AVt? — T? apse: ” cosh AVt2 — T? dt = nlp(XT) 
Sa! coe ae 2 pe ers [jz — 7? = Cag ? 
from which we find 
i ROE SE ae if ees 
>7 3 3=Vt? — T? do VP TP 
a/2 
= 2 i e~sreneda (49) 
and 
T . n/2 
h AVT? — #? 
= aa dt = rI,(AT) — 2 sf e ATsinédg. (50) 
Thus, 
oo T a/2 
i, Lr(t)dt = Sinh AT? [enon = 2, if esta, (51) 
and 


A =A(b,¢, T) = = exp {s won, (52) 


where \ = c — b, and 
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x/2 
u(x) = eft - —_ = { e*sin’dg ~ Vee 
sinh x sinh x zw Yo T 


x—o, (52a) 


We observe, since V7(t) = 1 for |t| = T, that the infinite product 
in (45) is, in effect, for each t, a finite product. That is, we may write 


sr(t) = Af{sgn Ses Iz Vr(t — tr) (53) 


to emphasize the localization of the reconstruction. 
Now, to estimate the error, we write, using (37) and (38), 


sp(t) — s(t) = (1/2){sgn s(t)}{exp[Ar(t)] — exp[h(¢)]} 











= (1/2)s(t)-{exp[Ar(t) — A(t)] — 1). (54) 
Now, 
F x? x3 
Je*—-1|= xt+oatat 7 
xl? Ix? 
< | |x| EEE y. -| sel#l — 1, 
Setting x = h(t) — h(t) and using the estimate (40) 
P P 1+e%7 
| hr(t) — h(t)| < 2 log [ae 
we have 
1 1+e%7 1+e%? 
Is(t) — sr) | = 5 1-5 - i} te = i} 
Qe>T 
<q eT | s(t) | (55) 


with equality only for t = t,, where s(t,) = 0. 
The error estimate in (55) is conservative since equality in (33) and 
(36) is attained only for (see Ref. 4) 


h(t) = +M(sgn t)-{sgn sin Avt? — T?}, Ree 4 
= 0, t<T. 


Here we have applied the inequality to a high-pass function h(t), the 
“meandering” sawtooth function of (26) with h(0) = 0, which is quite 
different from the above extremal function. However, the exponent 
AT is probably sharp. 

A series for computing L(¢) is given in the appendix. 
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VI. CONCLUSION 

We have seen that signals of a certain class can be reconstructed at 
any point ¢t with a relative error of the order of e~*? from the knowledge 
of their zeros in an interval (t — T, t + JT). The reconstruction 
algorithm can be applied to practical signals of the form 


s(t) = gp(t) + cos ct, 


where g,(t) can be closely approximated by a bandlimited signal, g(t), 
whose spectrum is confined to the interval [—b, b]. We would like the 
zeros of s(t) not to change much when g, is replaced by g. First, we 
require 


le(t)|<M<1 
for some suitable choice of M, say M = 0.75. We also need a condition 


on | gp(t)| and | g‘,(t) | to ensure that s(t) has only one zero between 
each extremal of cos ct, e.g., 


{gp(t)}? + {g>(t)/c}? < 1. 


This condition alone gives | g,(t)| <1, and hence, (—1)*s(kx/c) > 0, 
from which we conclude that s(t) must have an odd number of zeros 
in (kx/c, (k + 1)x/c). It is sufficient to consider (0, z/c). Suppose the 
number of zeros (counting multiplicity) in (0, z/c) is more than 1. 
Then, at some point x; in this interval (say, the second zero), we must 
have s(x,) = 0 and s’(x;,) = 0; i.e., 


8p(Xz) = —COS CXp, 

&p(Xp) ZC Sin cx, 
and then 

{ Sp(xe)}° + {'p(xe)/c}* = 1, 
which is a contradiction. Then, with the two above conditions, we may 
regard the zeros of s(t) as perturbations of the zeros of 
s(t; 6) = g(t) + cos ct, 
assuming that 
| g(t) — gp(t)| < «(b) <1 —M. 


Though the calculation is complicated, estimates can be made for 
the error in reconstructing g,(t) from the zeros of s(t). Of course, the 
relative error in | s(t) — s7(t)| will not be of the order of e~*” but of 
the order of | g(t) — g,(t)| (provided XT is sufficiently large). The 
algorithm has much in common with the reconstruction of almost 
band-limited signals from their samples. That is, if we were given 
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&p(k/c) and knew that 
| g(t) — gp(t)| < €(d), 


we would use the fact that 


= FR «(Ee (-F) 


for any function K(t) of the form 


sin ct 





K(t) = f(t) 

where f,(t) is a bandlimited function whose Fourier transform vanishes 
outside [—-A, A], \ = c — b > 0, and f,(0) = 1. Although the series 
would converge conditionally for f(t) = 1, it would not necessarily 
converge when g(kz/c) is replaced by g,(kz/c). So, it would be far 
better to take, for example, 


is T ~~ sin A\vt? — T? 
Ad) = sinh AT Ve2— T? 


and obtain an exponential localization in the reconstruction of g(t) 
from its samples and then apply the formula to g,(t). The basic 
difference in the problem solved here is that the zeros of s(t) give a 
nonuniform sampling of g(t), but we make use of the fact that the 
fundamental function h(t), given by (11), is high pass to obtain the 
exponential localization of the reconstruction. 

We should note that the validity of the reconstruction could be 
extended to signals of the form (3) with implicit restrictions on g(t) 
such that s(t) has only real simple zeros. This does not gain us anything 
from a practical viewpoint since we do not have a practical way to 
determine all admissible functions g(t), but it does say something 
interesting about the zeros of such functions, viz., that the zeros allow 
the exponential localization of the infinite products. The difference 
now is that we do not have the zeros occurring so regularly and we 
have to replace the upper bound for | h(t) | in (16) by something larger 
than z. This raises an interesting mathematical question: if s(t) is of 
the form 


s(t) = g(t) + cos ct 


and g is bandlimited to [—b, b], 0 < b<c< ©, and such that s(t) has 
only real zeros, what is the longest possible zero-free interval for s(t)? 
It is probably something like z/\, where \ = c — b. Also, what is the 
optimal truncation of the Hilbert transform kernel in order to estimate 
h(t) for the special class of h(t) encountered in this problem? 


280 TECHNICAL JOURNAL, FEBRUARY 1984 


Once a reasonable reconstruction algorithm has been found, one 
can ask for the best function V(t), where V7(t) = 1 for |t| = T and 
corresponding constant A such that 


s(t) _ A{sgn s(t)} ae Vr(t = tr) 


best approximates s(t) of the form 
s(t) = g(t) + cos ct 


for some reasonably smooth class of functions g(t), | g(t)| < M <1. 
For example, one might take typical speech signals g(t) and, with 
suitable c and T, find Vr empirically. 

Some interesting applications may arise when, instead of s(t), one 
transmits p(t)s(t), where p(t) is a low-frequency positive (nonnegative) 
signal. Then, s(t) can be recovered by clipping p(t)s(t) and applying 
the reconstruction algorithm. For example, in single-sideband trans- 
mission of x(t) it is advantageous to transmit the carrier cos ct for 
purposes of demodulation; i.e., s(t) is of the form (lower sideband) 


s(t) = (1 + x(t))cos ct + <(t)sin ct, 


where x(t) is bandpass, [a, b] and [—b, —a], 0 < a < b < c¢, and 
| x(t) | <1. Here, if x(t) is a speech signal, one could take p(t) to be a 
low-frequency (syllabic-frequency) envelope of x(t) and transmit 
p(t)s(t) to reduce the power in the carrier, and then recover s(t) at the 
receiver. The disadvantage is more sensitivity to noise, but this added 
sensitivity occurs only when p(t) is small, so p(t) could be detected at 
the receiver and used as a squelch control to mute the noise portions. 

It appears feasible to construct analog circuitry for effecting the 
reconstruction; the singularity in L(t) is not bad, and for practical 
purposes, L;(t — T) can be approximated by the impulse response of 
a realizable filter. The pulse train }\ 6(¢ — t,) can be ac coupled into 
the filter to obtain h7(t — T’) at the output, since the ac coupling gives 
(for sufficiently large time constant) very nearly (> 6(t — t,) — c/a) as 
input to the filter. The fact that (the approximation to) hr remains 
finite means that exp h(t) is not a good approximation to | 2s(t) | 
near the zeros of s(t); however, this is not important. Also, the problem 
of delaying sgn s(t) for multiplication by exp hr(t — T) can be avoided 
by using a threshold circuit to detect the “near-zeros” (approximately 
t. — T) of exp hr(t — T), and then, pulses derived from the threshold 
circuit can be used to generate sgn s(t — T). 

It may be desirable to use a feedback circuit consisting of an 
operational integrating amplifier and a multiplier connected so that 
the output y(t) of the amplifier satisfies 


y’(t) = hr(t)y(t). 
SIGNAL CLIPPING—I 281 


Then, 
y = M exp hr(t). 
Here, 
up(t) = 7 Y Kr(t — te), 


and one uses a filter whose impulse response approximates K7(t — T), 
except for t very near T. Additional means will probably have to be 
provided for controlling drift in the integrator (or equivalent offset in 
the multiplier). 

Given the state of the art, a specialized digital computer may be 
more practical than an analog computer for effecting the reconstruc- 
tion. The problem encountered here is the quantization of the location 
of the zeros. Some analysis is required, then, to determine admissible 
bounds on the quantization errors. 

For practical purposes, one need only compute s7(t) on a lattice, say 
t = kx/c. This is a good choice since s(k/c) # 0, and so, hr(k7/c) 
will be more easily computed in practice. One then would compute 


&rlkr/c) = sr(kx/c) — cos(kz) 


and then filter the samples to obtain an approximation to g(t). In 
analog computation this method eases the problem of approximating 
h(t) or Lr(t) since one does not require the value of Lr for small 
arguments. 

In connection with this last method, one may be given an unknown 
translate of sgn{s(t)}, so some method is required to determine the 
proper phase for sampling. It can be shown that 


20 
{sgn s(t)}-{sgn sin(ct + 6)} -—, where 5 <6< ; ; 
us 


is high pass, provided @ is small enough to allow the zeros of s(t) and 
sin(ct + @) to interlace. This implies that phase-lock circuitry can be 
used to adjust @ to zero, and thus determine the proper phase for 
sampling, viz., at the zeros of sin ct. 

Alternatively, one could construct the fundamental function h(t) by 
forming an unbiased integral of 


h'(t) =" V(t —t) —c¢, 


where 6(¢) is the Dirac delta function. An unbiased integral of a high- 
pass function is a particular integral that is also high pass. For practical 
purposes, a good approximation to h(t) can be obtained by forming 
a narrow pulse p(t — t,) at each zero t, of s and then ac coupling 
- ¥ p(t — tz) into an RC integrator. The zeros of h that occur between 
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the pulses are then the proper sampling points, the zeros of sin ct (see 
Refs. 11 and 27). 
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APPENDIX 
Series for Computing L;(t) 


Replacing t in (42d) by xT, we have 
Ly(xT) = log|x| + F(x; AT), Ix] <1 











= 0, lx| > 1, (56) 
where 
1 / po. os 7 — 72 5 
The integral in (57) can be evaluated with the aid of the identity 
hpv1—¢t? | 2/Q)n 
ee = 97 SEE toavelB), (8) 
where 
V1/261h/2(B) = = Bp ’ 
V1 /2BI3/2(8) = — = B - — 8 , 


Vi/2BI5/2(8) = (5+ +) sinh 6 — 2 cosh 8, 
and, in general, 
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Vi /2BIn+1p2(B) = (28)"{R[n + 1/2, —Ble’ — (-1)"R[n + 1/2, Ble}, 


where 


1 = (n + 1)! 2 (n + 2)! 3 
R tc +5 a slay OO ona 
a g 2Ml8)™ 
n! 
6 (nt kp). oe 
= Lae Mace 
and 
R{1/2, 6} = 1. 


Alternatively, if we set 


” sinh B 
= V0 4 DM art 59 


we have the recurrence relation, 
(2n + 1) 


n\P). 60 
B n(B) (60) 


dnai(B) = gn-1(B) _ 


In this notation we have 


_sinh BV1 — t* ¢ & (—6t?/2)" on(B) 





(einh @)V1 — 2 61 
(sinh 8) T-f fo» nt do(B)’ (61) 
and 
= (—6/2)" on(B) 
=~ 2 nyt 62 
F(0, 8) 2 Cnn! do(B)’ (62) 
and 
(—6x?/2)" bn(B) 
3 
Fla, 6) = FOB) + 3 rynt dB) ve 
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Signals Designed for Recovery After Clipping— 
Il. Fourier Transform Theory of Recovery: 
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This paper develops a Fourier transform theory for the recovery of signals 
of a certain class from their zeros. The class, denoted by S(b, c), consists of 
real-valued signals of the form s(t) = g(t) + cos ct, where g is bandlimited to 
[—b, b], 0<b<c< ©, and such that (—1)*s(ka/c) > 0, which is satisfied, for 
example, if | g(t) | <1. A very simple method of recovery is given for the case 
c > 3b, and a somewhat more complicated method is given for the case c > 2b. 
The theory also suggests a novel method of effecting amplitude modulation, 
which has advantages in high-power applications. 


I. INTRODUCTION 


Reference 1 shows that there is a practical way to reconstruct signals 
s(t) of a certain class S(b, c) from their zeros. The class S(b, c) consists 
of real-valued signals of the form 


s(t) = g(t) + cos ct, (1) 
where 
g is bandlimited to [—b, b],0<b<c<»o, (1a) 
and such that 
(-1*s(k/c) > 0, k = 0, +1, +2, ---. (1b) 


The alternation condition (1b) ensures that s has only real, simple 
zeros, one between each extremal of cos ct. A practical sufficient 
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condition for (1b) is 
|g(t)| <1, -wo<t<o, (1c) 


However, (1c) is not necessary. 
In Ref. 1 we introduced the so-called “fundamental function” asso- 
ciated with the zeros {t,} of s « S(b, c); viz., 


h(t) = J(t) — ct, (2) 


where J(t) is a jump function increasing by zw at each zero t;, J(0) = 
0. We showed that h(t) is a high-pass function with no spectrum in 
(—A, A), where 


A=c-b>0 (3) 


is the “gap frequency”, i.e., the difference between the top frequency 
of g and the carrier or “bias” frequency. We also established that 


—17r7<h(t)<a,-~<t<»o, (4) 


and derived the important relation 


ie ee | Shoe = i ae (5) 
T ST lami 7 
where / is the Hilbert transform of h. From (5) we obtained, using the 
fact that h is high-pass, an algorithm for approximating s(t) just from 
a knowledge of the zeros of s in the interval (t — T, t + T), the error 
in the approximation being of the order of e~. 

Here we want to consider alternate ways of reconstructing s(t) that 
involve additional facts about the Fourier transform of h(t). Although 
the functions we deal with will not, in general, have Fourier transforms 
in the ordinary sense, we can still attach consistent meaning to the 
statement “the Fourier transform of f vanishes over E,” where E is 
the union of a finite number of disjoint intervals. To do this (see Refs. 
2 and 3) we simply require f to be orthogonal to a suitable class of 
functions that have ordinary Fourier transforms supported on FL. Then 
we can say that “the Fourier transforms of f,; and f. agree over E if 
and only if the Fourier transform of (f,; — fe) vanishes over FE.” We 
will make free use of the results derived in Ref. 2 stemming from these 
basic definitions and also the results in Ref. 3 pertaining especially to 
high-pass functions. There is too much material needed for rigorous 
treatment to collect here. However, the reader should be able to follow 
most of the arguments by supposing the functions to have ordinary 
Fourier transforms or else by supposing the functions to be periodic. 

We need not concern ourselves with the nature of the Fourier 
transform of a bounded function in the sense of Schwartz distributions. 
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We only need the concept of the Fourier transforms of f; and fo 
agreeing over intervals, which has all the utility we require. 

Our first result is that for the case of large gap frequency, c > 30, 
we have 


git) = — = Y (cos ets)Ktt — te) (6) 


where K is in L; and 
ii K(t)e™ = 1, —b<w<=b 


= 0, jw] SAX -—b=c—2b>0b. (7) 


That is, if c > 3b, we can simply sample —cos ct at the zeros t, and 

then filter the impulses to recover g(t). Note that g(t,) = —cos cty. 
We next establish that the Fourier transforms of the fundamental 

function h and the function h, agree over (—(n + 1), (n + 1)A), where 


n (1)\k+t1 
nisin \3 sca fa agile Re oa 6) 
mk J 
and 
v(t) = 2e“g(t) + 7, (8a) 


[In (8), Im denotes the imaginary part.] The spectrum of the function 
y(t) is one-sided, being confined to the interval [\, 4 + 2b] and the 
point 2c. Hence, the Fourier transform of y*(t) vanishes over (—%, 


ky), Rk = 1, 2, --- . So in the case 2\ > ¢, ie., c > 2b, the Fourier 
transforms of h and h; agree over (—2), 2X), where 
h(t) = —2g(t)sin ct — sin 2ct. (9) 


Now low-pass filtering of h(t) is equivalent to filtering h(t), pro- 
vided the cutoff frequency of the filter is no greater than 2X. In case 
2\ > c we can filter h(t) with a sideband filter (upper vestigial) 
to recover a little more than the lower sideband of the signal 
—2g(t)sin ct, rejecting the higher frequency nonlinear components. 
Then standard sideband detection techniques allow us to recover g(¢), 
the details being given later. 

In the case b < c < 2b, nonlinear operations are required to recover 
g(t) or s(t) from h(é), i.e., from the zeros of s. The methods of Analytic 
Modulation Systems (see Ref. 2) are used here. We may take a filtered 
version of the analytic signal H(t), 


H(t) = h(t) + ih(t), (10) 
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namely, 


Hy (t) = fle H(x)K, (t — x)dx, (11) 

where the filter kernel K, belongs to L; and 
ies K(tje dt=1, ASwspr, (11a) 
y>c. _ (11b) 


Both functions here have one-sided spectrum, and since h is high- 
pass, the Fourier transform of H vanishes over (—%, )). By virtue of 
(11a), the Fourier transforms of H* and H agree over (—™, »). 

Next we define 


Gr(t) = exp{-tH7(t)}. (12) 


We then find that the Fourier transforms of G* and G agree over 
(—oo, vy), where 


G(t) = exp{—iH(t)} = 1 + 2e“g(t) + e7*. (13) 


Now, since v > c, we can filter G*(t), using a vestigial sideband filter 
(cutoff frequency v), which is equivalent to filtering G(t), so as to 
recover a little more than the lower sideband of 2g(t)sin ct (taking the 
imaginary part of G*(t)). Then we can proceed as before with the 
standard sideband detection to recover g(t). This method allows the 
minimum bandwidth requirement (y > c) in filtering operations on h, 
and thus eases the problem of obtaining the (unbounded) Hilbert 
transform A(t); i.e., h(t) + iA(t) is replaced by h(t) + ih#(t). Note, 
however, that in practice two filtering operations are required here, 
h — h* and h > hi, to effect the operation indicated in (11). 

Actually, the exponential function in (12) may be replaced by a 
partial sum of its power series, as detailed later. This option may be 
convenient in the case of moderately large gap frequencies where only 
a few terms are required. 

The sideband detection required in the previous method can be 
avoided when (11b) is replaced by vy > 2c. Then G¥(t) can be low-pass 
filtered (cutoff frequency v) to obtain (all of) G(t) = 2e“'s(t). Equiva- 
lently, the Fourier transforms of e“G#(t) and 2s(t) agree over 
(—p, u), where » = v —c >c. Denoting the real part of e“G*(t) by 
2si(t), we have 


s(t) = 5 lcos{h¥(t) + ctl} -explhs(e)). (14) 
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This is to be compared with the basic formula 


s(t) = 5 {sgn s(¢)Jexp{hte), (15) 
which follows from (5). Defining J*(t), analogous to J(t) in (2), by 
hy(t) = Jz(t) — et, (16) 
we may write 
s(t) = 5 {eos{J2(0)]} -exp{AX(e)} (17) 


and then (15), in a similar form, 
s(t) = 5 {eos[J(¢)}}-exp{h(e)}. (18) 


The function s#(t) may be low-pass filtered (cutoff frequency pu) to 
obtain s(t), the resulting formula being a generalization of the basic 
formula (15). 

Finally, a novel method of effecting amplitude modulation is sug- 
gested by the fact that for s in S(b, c), the Fourier transforms of the 
two functions, g(t)sin ct and (7/4) {sgn s(t)}-{sgn sin ct}, agree over 
(—2, 2). The square wave then may be filtered to obtain g(t)sin ct, 
provided 2 >c + borc> 30. 


Il. FILTERING THE INDUCED SAMPLES OF g(t) FOR THE CASE 
c > 3b 


The zeros of a signal s of the form (1), s(t) = g(t) + cos ct, induce a 
certain “oversampling” of g(t). That is, if we denote the zeros of s by 
{tp}, Ra/c < ty < (k + 1)a/c, we have 


g(t,) = —cos ct;, g bandlimited to [—), db]. (19) 


So, in effect, the zeros {t,} give us a nonuniform sampling of g(t) at 
the rate c/z, and we only need a rate slightly larger than b/z to pin 
down g(t). However, in the general case, nonuniform sampling com- 
plicates the recovery problem. 

We have given in Ref. 1 a reconstruction formula for s(t), based on 
(5), which may be regarded as one practical way of solving the non- 
uniform sampling problem (19). We could apply that formula to 
reconstructing only the uniform samples of g(t), say g(ka/c), and then 
use the fact that 


g(t) = - z &(ka/c) K ( ae: “/e), (20) 
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where K is any bandlimited function belonging to L; whose Fourier 
transform satisfies 


{ K(t)e"dt = 1, —b<w<=b 


= 0, |w| = 2c — bd. (20a) 


In engineering terminology, we low-pass filter the uniform samples of 
g. The validity of (20) for bounded g whose Fourier transforms vanish 
outside [—b, b] may be established from the fact (see Ref. 2) that K is 
a reproducing kernel for g, 


g(t) = { g(x)K(t — x)dx, (21) 
and the fact that (for each t) 
i] g(x)K(t — x)e“*dx = 0, |w| = 2c, (22) 


ie., as a function of x, g(x) K(t — x) is a function of L; whose Fourier 
transform vanishes (see Ref. 2) outside (—2c, 2c). Then (20) follows 
by applying the Poisson sum formula to g(x)K(t — x), 


{ g(x)K(t — x)dx = ; y & (kr/c) K (t — kr/c). (23) 


Now we would like to show that the nonuniform samples of g at the 
zeros of s can be treated in the same way as uniform samples in case 
c > 3b. 

We consider the logarithmic derivative of s, 

s‘(r) _ g’(r) — c sin cr 
s(r) g(7) + cos cr 





,THttw, (24) 


which may be written, multiplying numerator and denominator by 
2e*", as 
s'(r) _ 2e**g’(r) + ice" — ic 


= : 24 
s(z) 2e"g(r) +e" +1 eae) 


Then, since g’(7) and g(r) grow no faster than A cos br (see Ref. 4, 
Theorem 6.2.6, p. 83), we have 


a) a re 
Sony a) + ic = 0(e™”), u — ©, (25) 


From the infinite product representation of s(7) we have 


sr) _ $1 


26 
s(t) 21r-te _ 
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where the sum converges conditionally. But if we write 
s‘(t+iu)_ <= 1 
s(t+tiu) 2t-—trt+iu 
= (t — tz) <a u 
So ee ee ee ea Ae 
2(t-bP+w  A2(t-h)? +u? wn 


the second sum converges absolutely (see Ref. 4, p. 86). For any fixed 
u > 0 the second sum is a bounded function of t, since k a/c < t, < (k 
+ 1)z/c. Then (25) implies that for any fixed u > 0, the function of t, 


Py CET ae eee 


“ttre” (28) 


is a high-pass function whose Fourier transform vanishes over 
(—A, A). The integral of each of the pulses in the sum is z, and the 


pulses concentrate around the points ¢, as u > 0. Then, if f belongs to 
L, and its Fourier transform vanishes outside (—A, A), we have 


i; f(t)h’ (t; u)dt = 0, u> 0. (29) 
Then, letting u — 0, we have 


Fi peoax = 2 3 pe. (30) 


Now if \ > 20, i.e., if c > 3b, we can apply (30) to 
f(x) = g(x)K(é — x), 


where K is in L; and 
{ K(t)e"dt = 1, —b<w<b 
= 0, Jw] >A -—b=c— 2>b, (31) 
to obtain 


git) == ¥ g(tK(t - 6) 


ie : y (cos ct,)K(t — th)  (¢ > 30). (32) 


That is, the nonuniform samples of g at the zeros of s, occurring at 
the rate c/7, may be filtered as if they were uniform samples, provided 
c > 3b. The rate of nonuniform sampling must be three times greater 
than the Nyquist rate for g in order to obtain this simple method of 
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recovery. We merely sample cos ct at the zeros of s(t) and then use 
appropriate low-pass filtering to recover g(t). In this connection, if the 
phase of the carrier, or high-frequency bias, cos ct is not known, it can 
in practice be obtained from phase-lock circuitry operating on sgn s(t). 
It has already been noted in Ref. 1 that [cf. eqs. (82) and (83) here] 


{sgn s(t)}-{sgn sin ct} 


is high-pass, which gives a proper condition for phase lock. 

In connection with the necessity of the condition c > 3b, we note 
that no matter how small g is, provided g + 0, the nonuniform 
sampling, though almost uniform, is not perfect in the sense that (32) 
does not hold (exactly) unless c > 3b. We could estimate the error in 
(32) for c > b, which, with proper choice of K, will be small for small 
g (depending, of course, on how large the gap frequency is). Depending 
on the application, (32) may give a suitable approximation to g(t) for 
c only somewhat larger than b. 

We note that in (28), letting u — 0, we may write, formally, 


PADS Peas (33) 


where 6(t) is the Dirac delta function and h’(t) is interpreted as the 
derivative of the fundamental function h(t) defined in (2). So far we 
have described the Fourier transform of h(t) only as vanishing over 
(—A, A). We give a further description next. 


Ill. THE FOURIER TRANSFORM OF THE FUNDAMENTAL FUNCTION 


The complex-valued function, 
H(t) = h(t) + ih(0), (34) 


where h is the fundamental function and h is the Hilbert transform of 
h, has a one-sided spectrum, i.e., its Fourier transform vanishes over 
(—oo, X). Equivalently (see Ref. 2), it extends as a function analytic in 
the upper half-plane of the complex variable 7 = t + iu, satisfying 


H(t + iu) = 0(e~”), u—> 0 (35) 
(—o <t< 0), 
As Ref. 1 establishes, H(t) is related to s(t) by 
H(r) = i log G(r), (36) 
where 
G(r) = 2e" s(r), (37) 


and we take the principal branch, log(1 + z) ~zasz—0. 
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The function G(7) is an entire function. Let us set 
G(r) =1 + y(r) = 1 + Qe” g(r) + e?*, (38) 
where 
y(t + iu) = Ofe~™}, u—> ©, (38a) 


For sufficiently large u, say u > Uo, we will have | y(t + tu)| <1 for 
all t so that the Taylor series for log G converges, i.e., 


H(r) =i Lae 5) + 5G) + | (Im 7 > wo). (39) 


The series does not converge everywhere on the real line because we 
do not have | y(t) | < 1 for —« < t < , However, if we define 


; n (—1)*t! P 
H,(7) Sb tag =, AY (r), n= 1; 2; ry (40) 
pat OR 
we have 
H(r) — H,(r) = 0 fy"""(7)}, =—-y 2 0 
= 0 {ety lu, (41) 


It follows (see Ref. 2) that the Fourier transforms of H(t) and H,(t) 
agree over (—®, (n + 1)A). Also, the Fourier transforms of H(t) and 
H,,(t) vanish over (—®, \). Then for any finite n we may write 


H(t) = Hy(t) + Rasi(t) 
=h,(t) + ihp(t) + Ras(t), (42) 


where the Fourier transform of R,+;(t) vanishes over (—%, (n + 1)dA), 
and h, is the Hilbert transform of hn. 

So, despite the fact that the infinite series in (39) does not converge 
for u = 0, we can still equate the Fourier transforms of H(t) and the 
partial sums H(t) of the series over the frequency interval (—, 
(n + 1)A). Here again we see the importance of a positive gap frequency 
A. Since 


y(t) = 2e“g(t) + e?, (43) 


we have an effective way of describing the Fourier transforms of h(t) 
and A(t) over intervals (—(n + 1)A, (n + 1)dA); i.e., they agree with the 
Fourier transforms of h,(t) and h,(t), respectively, over such intervals. 
These functions in turn are simply related to g through y. For example, 


Hiy(t) = iy(t) 
= Qie™g(t) + ie? (44) 
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1 a 
h(t) = 5 Hy(t) + 3 Hi, (t) 
= —2g(t)sin ct — sin 2ct, (45) 
x 1 1 
h(t) = 3 Ai,(t) - 3 H,(t) 


= 2g(t)cos ct + cos 2ct. (46) 


Thus the Fourier transforms of h(t) and —2g (t) sin ct agree over the 
frequency interval (—2), 2A). Similarly, the Fourier transform of 
h(t) = log 2| s(t) | agrees over (—2X, 2) with the Fourier transform of 
2g(t)cos ct. Since the Fourier transform of g(t)sin ct vanishes outside 
[—b — c, b + c], we see that in case 2\ > b + ¢, i.e., c > 3b, we can 
recover —2g(t)sin ct by appropriate filtering of h(t). Actually, in case 
c > 3b, it is simpler to recover g(t) by the sampling technique of the 
previous section. 

Next we consider the case 2\ > ¢, i.e., c > 2b, in which case there is 
a relatively simple way of recovering g(t) from h(t). 


Iv. A LINEAR RECOVERY METHOD FOR THE CASE c > 2b 


Here we need only be concerned with the case 2b < c S 3b, since for 
c > 3b we have a simpler recovery method. However, in case c > 2b we 
still have a method requiring only linear operations. 

We have seen that the Fourier transform of the fundamental func- 
tion h(t), defined in (2), agrees over (—2A, 2) with that of h(t) in 
(45) and hence, agrees over (—2\, 2A) with that of —2g(t)sin ct. In case 
2X > ¢, ie., c > 2b, the band (—2), 2) includes the lower sideband of 
—2g(t)sin ct and at least a part of the upper sideband. Then, with 
proper filtering of h(t) and subsequent demodulation, we can recover 
g(t) by linear operations. In the first place the filtering must eliminate 
everything outside the band (—2), 2A); i.e., the higher-order nonlinear 
terms h,(t), n = 2, defined in (42). There is nothing in the band 
(—A, A) since the Fourier transform of h; vanishes over (—), A). So let 
kK, denote the kernel (impulse response) of a bandpass filter, where 
we first require 


{ K(t)e“dt = 0, Jw | = 22, (47) 


{ | Ky(t) |dt < ~, (47a) 
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and define 


f(t) = { h(x)Ki(t — x)dx, (48) 


which is equivalent to 


A(t) = - I 2g(x)sin cx Ky(t — x)dx 


=— ii K,(x)2g(t — x)sin {c(t — x)}dx. (48a) 


Now if we multiply f,(t) by —sin ct and band limit the result, we should 
obtain a filtered version of g(t). The spectrum of /f; is confined to the 
intervals [A, 2A] and [—2A, —A]. Then, defining 


f(t) = —fi(t)sin ct, (49) 


we see that the spectrum of f, is confined to the union of the intervals 
[A +c, 2\ +c], [A — c, 2A — c], [-2A — c, -A — c], and [-2A + ¢, 
—r +c]. (Recall that \-= c — 6.) The innermost intervals are the 
second and fourth; viz., [—b, c — 2b] and [—c + 2b, b]. The outer 
intervals are [2c — b, 3c — 2b] and [—3c + 2b, —2c + b]. We are 
interested in the case 2b < c < 3b and want to retain the spectrum in 
the inner intervals and eliminate the spectrum in the outer intervals 
by low-pass filtering. We have ample room for this. So let Ke be a 
kernel (impulse response) of a low-pass filter, satisfying 


i K,(t)e““dt = 0, Jw] = 2c — b> 3b (50) 


{ | Ko(t) |dt < , (50a) 


and then define 


gilt) = B fo(x)Ko(t — x)dx. (51) 


We have 


£i(t) = = | sin cx f,(x)Ko(t — x)dx 


= | sin cx K2(t — x)dx { Ki(y)2g(x — y)sin{e(x — y)}dy. 
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The requirements (47a) and (50a) allow us to interchange the order of 
integration. Then 


{ Kilt — x)g(x — y)2 sin cx sin{c(x — y)}dx 
= i K,(t — x) g(x — y) {cos cy — cos(cx — cy)}dx 


= (cos cy) i} K2(t — x)g(x — y)dx, (51a) 
the other term dropping out since, for each y, the Fourier transform 


of g(x — y)cos(2cx — cy) vanishes over (—2c + b, 2c — b). Thus we 
have 


gi(t) = i i &(x)K3(t — x)dx, (52) 


where 
K3(t) = i) 7 {Ki(y)cos cy}-K2(t — y)dy. (52a) 


> 
So gi(t) is in fact a filtered version of g(t). Denoting by K;(w) the 
Fourier transform of K;(t), the overall filter response is given by 


> 1 [> > > 
K3(w) = 5 {ite +c) + Ki(w - o} Reta (53) 
So, if in addition to the requirements (47) and (50) we require 
1 [> > > 
9 {cw +c) + Ki(w - of Rae) = 1, |w| <b, (54) 


then we will have, with g; defined in (51), 


&i(t) = g(t). (55) 


Conceptually, the simplest way to satisfy (54) is to require K, to be 
a vestigial sideband filter with linear cutoff characteristic, i.e., assum- 
ing 2b<c < 3b, 





> 1 2A -—w 
=. —_ <= — 
Kio) =5, 2 = sle el Se - 2b 
= 1, AX Sw S 2b, (56) 
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and, of course, 
> > 
Ky(w) = Ki(-w) 


> 
K,(w) = 0, Jw] > 2d. 


Then Kz is a low-pass filter with gain of 2 in the passband; i.e., 


> 
=0, Jol =2ce—b>3b. (57) 


We should note that the multiplication of f,(t) by —sin ct in (49) can 
be replaced by multiplication with —7/4 sgn{sin ct}, since the addi- 
tional high-frequency terms drop out in (51a). That is, a switch-type 
multiplier may be used in the demodulation process, thereafter making 
appropriate gain corrections. 


V. A GENERAL RECOVERY METHOD 


In case A is not greater than b, i.e., c is not greater than 2b, nonlinear 
operations are required to recover s(t) or g(t). The general method 
here, which makes further use of the ideas developed in the Theory of 
Analytic Modulation Systems (see Ref. 2), affords an alternative to 
the method in Ref. 1. In effect, the method here allows the Hilbert 
transform in (11) to be replaced by a filtered version. This modification 
must be taken into account.in subsequent operations. 

The method is best explained in terms of the functions with one- 
sided spectrum, i.e., the functions G and H in (36), 


H(t) = i log G(t) = i log[1 + y(t)] = A(t) + ih(t), (58) 
where 
G(t) = 2e%s(t) 
and 
y(t) =1 + Qe“ g(t) + 7". 
We have already noted that the Fourier transforms of H(t) and H,(¢), 
defined by 
(es B ae 
k 
vanish over (—°, \) and agree over (—%, (n + 1)A). 
Let us denote by H(t) any function such that the Fourier transforms 
of H(t) and H(t) agree over (—%, v). In particular, H¥ may be any 
function of the form 


Ay(t) = . A(x) Kt — x)dx, (60) 





H(t) = iD yt), (59) 
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where K,, in the interesting case v > X, satisfies 
{ ; K(t)e“dt =1, ASwsp, (61) 
with 
ix | K,(t) |dt < ©, (61a) 
Or we may take H*(t) to be any function of the form 
Hy) =i Y an", (62) 


where the a; decrease rapidly enough and 
fs | k+1 





a, = forl<k<n, where (n+1)AZ2v. (62a) 
Now let us denote by G*(t) any function such that the Fourier 
transforms of G(t) and G*(t) agree over (—©, v). Since 


G(t) = 1 + 2e“g(t) + er, 


and g(t) is real, we can recover g(t) from any such function G*(t) 
provided v > c, which is to say all we need is a little more than the 
lower sideband of g(t)sin ct or g(t)cos ct. 

We have the inverse relation between G(t) and H(t); viz., 


G(t) = exp{-iH()}, (63) 


and as shown in Ref. 2, the Fourier transforms of G(t) and 
exp{—1H*(t)} agree over (—0, v). So we may take 


Gk(t) = exp{—1H#(t)}, yp>c, (64) 


where H7(t) is any function of the form (50). 

Since the Fourier transforms of H(t) and H(t) vanish over 
(—oo, \), the Fourier transform of the nth power of H¥(t) vanishes over 
(—0o, nA). So we may take, if we prefer, 


n 


Gi(t) = > = [iH (e)}* (nt+l)Xzv>Cc. (65) 
r=0 R! 


This last form may be preferable to (64) for moderately large gap 
frequencies, for example, 3\ >c. For small gap frequencies there would 
appear to be no practical reason for perferring a large-degree polyno- 
mial in (65) to the exponential in (64). In general, if 2\ < c, then 
nonlinear operations are required to obtain a suitable G*(t). We can 
modify this statement in case g(t) itself has a spectral gap about the 
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origin. Then we could recover g(t) from G*(t) when » is slightly less 
than c. For example, in this case if 2\ = c, we could take n = 1 in (65). 

Whether we choose (64) or (65), we have the Fourier transforms of 
—Im{G*(t)} and —2g(t)sin ct agreeing over (—y, v), vy > c. Then the 
recovery problem is essentially that of Section IV, except there the 
Fourier transforms of h(t) and —2g(t)sin ct agree over (—2A, 2A), 
2 > c. So, with obvious modifications of the cutoff frequencies of the 
filters, the recovery process in Section III applies to the general case 
when h(t) is replaced by —Im{G#(t)}, v > c. 

We have 


H¥(t) = he(t) + ih¥(t). (66) 


So the general recovery method allows A(t) to be replaced by a function 
h* (t), whose Fourier transform agrees with that of A(t) only over the 
interval (—v, v), v > c. In particular, h(t) could be bandlimited to a 
larger interval. Alternatively, h*(t) could be obtained from h(t) by a 
transversal filter having the appropriate periodic frequency response. 

In practice, all the required operations can only be approximated. 
As contrasted to the method in Ref. 1, complexity in the recovery 
method here is exchanged for simplicity in approximating the Hilbert 
transform of h(t). 

We now want to specialize to the case vy = 2c and obtain a formula 
for s(t) analogous to that in Ref. 1. 


VI. A GENERALIZATION OF THE BASIC FORMULA 
From (5) we have 


s(t) = 5 fsgn s(¢)}-exp{At)}. (67) 


We call this the basic formula for recovery of s(t) from sgn s(t) or, 
equivalently, from the fundamental function h(t) defined in (2). 

In Ref. 1, the practicality of (67) was demonstrated by using the 
fact that the Fourier transform of h(t) vanishes over (—A, ). Here we 
give a generalization of (67), which may be preferred in analog imple- 
mentations. 

We take H(t) to be any function of the form (60) where v > 2c, and 
set 


Gr(t) = exp{-iH;(é)}, (68) 


as in (64). Now we have the Fourier transforms of G(t) and G*(t) 
agreeing over (—%, v) where vy > 2c. Recall that 


G(t) = 2e*s(t), 


and, therefore, its Fourier transform vanishes outside [0, 2c]. So we 
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can multiply G*(t) by e~“ and then, with appropriate filtering, recover 
2s(t). 
Let us define then 


Qsi(t) = Refe““GF(t)}. (69) 


Since s(t) is real, the Fourier transforms of s(t) and s*(t) agree over 
the interval (—p, ») where p» =v —c>c. 
In (68) we set 


Hi(t) = h(t) + ihz(t), 
and then write (69) as 


st(t) = 5 {eosfht(t) + ctl} -exptht(e)}, (70) 


and note the similarity to (67). To carry this a bit further, we recall 
the definition (2) of h, 


h(t) = J(t) — ct, (71) 


where J(t) is a jump function increasing by z at each zero of s(t), 
J(0) = 0. We may define J*(t) in an analogous way by setting 


h*(t) = J*(t) — ct. (72) 
Since h*(t) is a filtered version of h(t), their Fourier transforms 
agreeing over (—», v), we may think of Ji*(t) as a filtered version of 
J(t), although the corresponding convolution with J(t) alone need not 
make sense. 
Now we may write 
1 ‘3 
si(t) = 3 {cos[J7(t)]}-exp{hr(t)}, (73) 
and observing that {sgn s(t)} = cos[J(t)], 
1 ‘ 
s(t) = 3 {cos[J(t)]}-exp{h(t)}. (74) 


This is the analogy we seek. Now s(t) may be obtained by bandlimiting 
sa(t); Le., 


s(t) = { s*(x)K.,(t — x)dx, (75) 


where K,,,(t) belongs to L; and satisfies 
i, K,, (tle “dt =1, |w|se 


=0 lw] pre. (76) 
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We may regard (74) as a limiting case (u — ©) of (75). That is, (75) is 
a generalization of (74). 


VII. NOTE ON HIGH-POWER AMPLITUDE MODULATION 


Compared to conventional modulators, switching-type modulators 
have the advantage that little power is dissipated in the switches (e.g., 
transistors). Of course they also have the advantage of easily effecting 
linear modulation. The disadvantage is the introduction of high- 
frequency components, which must be removed by filtering. However, 
if the filter has low losses, the high-frequency energy is, for the most 
part, stored in the reactive elements of the filter, resulting in an 
efficient modulation system. 

If g(t) is bandlimited to [—b, b] and we wish to modulate sin ct with 
g(t), then a straightforward method to do this with a switching-type 
modulator is to form 


f(t) = g(t) {sgn sin ct} 


= - ecsin ct + 5 &(t)sin 3ct + -- |. (77) 


and then, assuming c > b, an output filter may be designed to deliver 
4/x g(t)sin ct to the load. 

In high-power applications, the disadvantage of this method is that 
the input g(t) to the switches must furnish the power to the load. That 
is, high-power baseband amplifiers are required. An interesting alter- 
native is suggested by the results in Section III. 

There we showed that the Fourier transforms of h(t) and 
—2g(t)sin ct agree over (—2A, 2X) where \ = c — b. 

Recall that 


h(t) = J(t) — et, (78) 


where J(t) is a jump function with J(0) = 0, which increases by z at 
each zero of the function s(t), 


s(t) = g(t) + cos ct. (79) 
Now denote by a(t) the periodic sawtooth function, 


o(t) = 5 — et, 0<t<za/ec 


a ( = | = a(t). (80) 
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This function, associated with log(1 — e’*“), has the Fourier series 


oo 1 . 
o(t) = ¥ — sin 2nct. (81) 
n=1 NV 
Since s(t) has one and only one zero between each extremal of 
cos ct, we have 


J(kx/c) = kr, k=0, +1, +2,.--. 
Now the function ¢(t), defined by 
p(t) = h(t) — o(t) = J(t) — ct — oft), (82) 


is a pure step function increasing by a at the zeros {t,} of s(t) and 
decreasing by z at the zeros of sin ct. We see that ¢(0+) = —7/2 and 
since the zeros of s(t) and sin ct interlace, it is easy to see that we can 
express ¢(t) in the form 


g(t) = — ; {sgn s(t)}-{sgn sin cfé}. (83) 


Clearly, from (81) and (82), the Fourier transforms of ¢({t) and h(¢) 
agree over (—2c, 2c), and hence, the Fourier transforms of ¢(¢) and 
—2g(t)sin ct agree over (—2), 2A), since 2c > 2\ = 2(c — b). 

In order to filter out —2g(t) sin ct from ¢(t) we require c + b < 2d 
or c > 3b. Thus, if c > 3b we can use (83) as a keying or switching 
voltage and generate 


f(t) = M {sgn s(t)}-{sgn sin ct} 
= ee g(t)sin ct + hi-freq, (84) 


where M represents a constant voltage that furnishes power to the 
load and “hi-freq” represents components of f(t) outside the band 
(—2d, 2) that may be removed by filtering. 

Thus, if we add cos ct to g(t), and simply require | g(t) | = 1, with 
c > 3b, where b is the “top-frequency” of g (i.e., add a high-frequency 
bias larger than | g(t)|) and then clip the result to obtain {sgn s(t)} 
and then multiply by {sgn sin ct}, we obtain a square wave, which we 
may use to switch transistors alternately between +M, the output 
having the desired “in-band” component (4M/z) g(t)sin ct. 

Alternatively, we may multiply s(t) by {sgn sin ct} and then clip the 
result to obtain the keying or switching voltage. 

This method clearly has an advantage over the straightforward 
method of (77) in high-power applications. The disadvantage of both 
methods is that filters capable of handling high power are required. 
Conceivably, in some applications (audio) the filtering may not be 
required. 
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VII. DISCUSSION AND CONCLUSION 


Here we have, in essence, detailed the remarkable properties of the 
zeros {t,} of any function s(t) of the form (1). We may regard the 
description of the Fourier transform of the fundamental function 
h(t) = J(t) — ct as summarizing the properties of the zeros {¢;,}. 

The rude property is that the Fourier transform of h vanishes over 
(—A, A), where X is the gap frequency in (1). This means that 


: h’(t) = 2 5(t — t,) — c/x, 


is high-pass, and so provides a quadrature formula for functions f in 
L, and bandlimited to (—A, A), i.e., 


{ _a(t)dt = - 5 f(ts). 


This formula is the basis of the simple recovery procedure for c > 3b, 
which is not an uninteresting case, e.g., in recording applications, 
where c is the bias frequency. 

In the first stage of refinement, the zeros {t,} have the property that, 
insofar as regards spectrum in the interval (—2A, 2d), the transfor- 
mation s(t) — h(t) or s(t) — }*.. 6(t — ty) is linear. Then, when 2) > 
c, i.e., c > 2b, we have a linear recovery method that we would use in 
case c were not larger than 3b. What is involved here is essentially 
vestigial sideband detection. 

In further stages of refinement, the zeros {t,} have the property that 
we can equate the Fourier transform of h(t) (or the pulse train) over 
frequency intervals (—(n + 1)A, (n + 1)A) to the Fourier transform of 
a finite sum involving no more than the nth power of g(t). For recovery 
of g(t) in case c < 26 the theory of analytic modulation and detection 
after bandlimiting as developed in Ref. 2 comes into play. The theory 
given there is specialized here to detection of analytic signals z(t) of 
the special form z(t) = 2e%'s(t) = G(t), where s(t) is a real-valued 
bandlimited signal of still more special form. The analytic “modulation 
law” f(z) in this case is i log z. We are given 


h(t) = Re[z log{z(t)}], 
and then can obtain z(t) from 
2(t) = exp{—-iH(t)}, H(t) = h(t) + ih(t). 


This is essentially the basic formula in Ref. 1. However, to allow us as 
much latitude as possible in approximating h(t), we assume, in effect, 
that we are given hx(t), a function whose Fourier transform agrees 
with that of h(t) over the frequency interval (—»v, v). We only require 
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v >c, half the bandwidth of z(t), because of the special form of z(t). 
The general recovery method assumes no more than py > c, which 
requires the vestigial sideband detection techniques in recovery of g(t) 
from 


2*(t) = Gk(t) = exp{-iH¥(t)}, Hi (t) = he(t) + ihk(t). 


The general recovery method is interesting from a mathematical 
viewpoint, since we assume no more than is necessary. However, the 
simplification obtained by requiring v > 2c, as in Section VI, appears 
to be well worth the price paid in filter requirements, i.e., the filters 
are required to have prescribed characteristics over the band (—», v). 
The resulting recovery formula may be preferred to the basic formula 
given in Ref. 1, especially for analog implementation. 


REFERENCES 


1. B. F. Logan, Jr., “Signals Designed for Recovery after Clipping—I. Localization of 
Infinite Products,” AT & T Bell Lab. Tech. J., this issue. 

2. B. F. Logan, Jr., “Theory of Analytic Modulation Systems,” B.S.T.J., 57, No. 3 
(March 1978), pp. 491-576. 

3. B. F. Logan, Jr., “Properties of High-Pass Signals,” doctoral thesis, Elec. Eng. 
Dept., Columbia Univ., New York, 1965. 

4, R. P. Boas, Jr., Entire Functions, New York: Academic Press, 1954. 


AUTHOR 


Benjamin F. Logan, Jr., B.S. (Electrical Engineering), 1946, Texas Tech- 
nological College; M.S., 1951, Massachusetts Institute of Technology; 
Eng.D.Sc. (Electrical Engineering), 1965, Columbia University; AT&T Bell 
Laboratories, 1956—. While at MIT, Mr. Logan was a research assistant in 
the Research Laboratory of Electronics, investigating characteristics of high- 
power electrical discharge lamps. Also at MIT he engaged in analog computer 
development at the Dynamic Analysis and Control Laboratory. From 1955 to 
1956 he worked for Hycon-Eastern, Inc., where he was concerned with the 
design of airborne power supplies. He joined AT&T Bell Laboratories as a 
member of the Visual and Acoustics Research Department, where he was 
concerned with the processing of speech signals. Currently, he is a member of 
the Mathematical Research Department. Member, Sigma Xi, Tau Beta Pi. 


306 TECHNICAL JOURNAL, FEBRUARY 1984 


AT&T BELL LABORATORIES TECHNICAL JOURNAL 
Vol. 63, No. 2, February 1984 
Printed in U.S.A. 


Integrated Voice/Data Services on Fasnet 


By J. W. MARK* and J. O. LIMB? 
(Manuscript received July 8, 1983) 


Fasnet is a high-speed local area network with a pair of unidirectional 
communication channels. While individual stations control their own access 
to the medium, global scheduling of the medium for various types of traffic is 
controlled centrally. Centralized scheduling permits a slot to be endowed with 
a type, so it is feasible to integrate synchronous and asynchronous traffic 
services on Fasnet. In this paper we describe a policy for admitting new voice 
calls into the system and an approximate analysis of the queueing behaviour 
of data packets in an integrated voice and data services environment. Queueing 
behaviour can also be studied by means of simulation, which is, however, an 
expensive undertaking, particularly at high transmission rates, for example, 
at 100 Mb/s. The approximate method described in this paper allows the 
performance evaluation for high-speed local area networks. 


I. INTRODUCTION 


Packet-switching technology has enhanced the potential for com- 
plete integration of different types of traffic on a common transmission 
facility, more particularly within a local environment where signal 
propagation delays are relatively small. Over the past decade a number 
of papers have been published on the integration of voice and data.’® 

We will address the problem of integrated voice and data services 
using a specific local area network (LAN), Fasnet,’ as the backbone 
support. The techniques explored and evaluated here would carry over 
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to other related systems (e.g., certain implicit token-passing systems 
and some ring networks) with small modifications. 

Fasnet is a dual-channel, packet-switched LAN. Channel time is 
partitioned into contiguous slots, each of a duration long enough to 
transmit one packet of data plus the necessary overhead to allow for 
access control, addressing, error checking, etc. A voice source, engaged 
in a conversation, may be silent (either listening to the other party in 
the conversation or in a short pause between speech spurts) or active 
(producing a speech spurt). An active voice source is partitioned into 
segments, each of length T' seconds. Each segment of voice is mapped 
into one information packet. The length of an information packet is 
thus determined by the code rate of the incoming voice source and the 
transmission rate of Fasnet. To maintain voice continuity, it is nec- 
essary to transmit one voice packet of the same source every T seconds. 
From the LAN’s point of view, to maintain an approximate continuous 
service, it must offer service to voice once every T' seconds. The initial 
portion of the T-second interval will be devoted to serving voice calls. 
The interval T between consecutive starts of voice service will be 
referred to as a service cycle. The service cycle period can be expressed 
in terms of the number of slots, and it is conceivable that the number 
of voice sources connected to Fasnet is larger than the number of slots 
in one service cycle. Obviously, the system cannot accommodate more 
simultaneously “active” voice source than the number of slots available 
in one service cycle without degrading the signal in some way; some 
kind of admission control needs to be exercised. 

In this paper we assume that a station may be handling voice, data, 
or both. If we consider voice only, a station that is not in a voice 
session will be referred to as idle. A nonidle station producing a speech 
spurt (talk spurt) will be referred to as active; it will be referred to as 
silent if its voice session is in a silent condition. 

From an operational point of view, calls already in progress, whether 
they be in a talk spurt or in a silent interval, must be guaranteed 
channel time when needed. Since data can be buffered and voice needs 
continuity, voice will be given a higher priority than data. New calls 
that are not admitted will be blocked. In the event that a subscriber 
is blocked, an appropriate tone (short busy) is returned to the sub- 
scriber and he/she would try again at a later time. (Of course a 
satisfactory system design would provide sufficient channel capacity 
that such an event would occur with no more than a prespecified 
probability.) On the other hand, unused voice capacity is available to 
data traffic until there is an increase in voice traffic. In this paper we 
propose an admission policy for new calls under the constraint that 
ongoing sessions are guaranteed channel access and that the average 
number of active voice calls is controlled to within a small neighbour- 
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hood of the target operating value, N;. (Due to statistical variations 
in the number of active voice calls and the fact that ongoing voice 
sessions cannot be dropped, the mean number of active voice calls can 
only be controlled to within a neighbourhood of N;.) We evaluate the 
performance of data under this admission policy. An alternative strat- 
egy has been explored in Ref. 10, in which data traffic is guaranteed a 
minimum capacity. Voice stations returning from the silent state to 
the active state are not guaranteed a voice slot, and clipping of a 
speech spurt can occur as experienced in an overloaded Time Assign- 
ment Speech Interpolation (TASI) system."! The maximum voice load 
is thus determined by the amount of clipping that may be tolerated. 
The study showed that, in the distributed environment of the LAN, 
performance approached very closely that of an ideal TASI system. 

The priority admission mechanism and the data queueing behaviour 
in Fasnet are similar to those in Welnet.® The difference lies in the 
media access control, the network architecture, and the method of 
switching from one traffic type to another. In Fasnet, a “dead time”, 
referred to as a sojourn time, usually accompanies switching from one 
type of service to another. In Welnet, provided the packet transmission 
time is greater than the propagation delay, there is no dead time due 
to switching of traffic type. Otherwise, the service disciplines of Fasnet 
and Welnet are similar. 

Allowing voice sessions to capture the channel when needed will 
cause the service rate for data to fall below the data arrival rate during 
some epochs when the system is heavily loaded with voice and data. 
When such an event takes place, the traffic intensity, defined as the 
ratio of data arrival rate to service rate, will temporarily exceed unity. 
In this mode of operation, the aggregate data queue will experience a 
sustained growth with a nonstationary distribution, and the queueing 
process will be temporarily unstable. An approximate analysis method, 
based on a periodic switching single-server model and the fluid ap- 
proximation method,” was developed in Ref. 8 to analyze the data 
queueing behaviour for integrated voice and data services on Welnet. 
The approximate analysis method involves a decomposition of the 
queueing process into stationary and nonstationary queues and a 
determination of the mean queue contribution from each of the oper- 
ational modes. With a suitable modification to account for the sojourn 
time, the approximate analysis method developed in Ref. 8 is used to 
analyze the data queueing behaviour for integrated voice and data 
services on Fasnet. 

This paper is organized as follows: Section II briefly describes the 
Fasnet architecture, which includes a description of the frame format 
of a Fasnet slot and the concept of a voice cycle. Section III discusses 
the evolution of the state equation for voice, which leads to the 
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proposal of an ad hoc admission policy for new calls in Section IV. 
Fasnet provides centralized control of classes of traffic, where a class 
may be a traffic type, a priority, or some combination. However, access 
by individual stations is a distributed function, being determined by 
the station itself. It is important to assure that the centralized control 
is implemented in such a way that service to each traffic type remains 
fair. (We define a fair system as one in which all active stations have 
an equal chance to access the common channel.) In Section V we 
evaluate the performance of data for a “fair” Fasnet, under the new 
voice admission policy introduced in Section IV and the constraint 
that ongoing voice calls are guaranteed service. 


Il. FASNET ARCHITECTURE 


Fasnet uses two channels with unidirectional signal propagation in 
each. The first station (the head station in Fig. 1a) in each channel 
exercises centralized access control. In the event the first station 
malfunctions, the second station senses the lack of timing information 
and assumes the role of the first station. Since there is complete 
symmetry in the two channels, we will consider the information flow 
in one of them only. 

Each channel has M taps, each of which supports a signal processor 
that handles both transmit and receive functions. The processor will 
be referred to as a station. Each station will support a number of user 
machines or processes. The first station permits access to each class 
by issuing a slot of a given type. A new class of traffic may be admitted 


FORWARD PATH (CONTROL IN S AND B SUBFIELDS) 


Re 


HEAD ‘TH END 
STATION STATION STATION 


4) 


REVERSE PATH (CONTROL IN E SUBFIELD) 





(a) 


BED hee 
ACCESS | 
eee FIELD ~~?" —-—-——" INFORMATION FIELD-————— . 


(b) 
Fig. 1—Fasnet: (a) topology; (b) frame format. 
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by sending a request to the first station, which, in turn, upgrades its 
traffic classification to accommodate the new class. 


2.1 Fasnet frame format 


Channel time is divided into slots, each of a duration sufficiently 
long to transmit an information packet plus control, addressing, and 
cyclic redundancy checking information. A Fasnet slot comprises an 
access field and an information field in a manner depicted in Fig. 1. 
The access field is endowed with an S (Start), an E (End), and a B 
(Busy) subfield. The S subfield contains a code that designates the ~ 
start of a sequence of slots of a particular traffic type. The B subfield 
is a single-bit subfield that indicates whether the information field is 
full or empty. Individual stations have the capability to read and write 
into the B bit. Guard times on either side of the B subfield allow for 
the read and write deliberations. As shown in Fig. 1, the contents of 
the information field are under the control of the individual stations. 
A slot is initiated with an empty information field and the content of 
the B subfield is a “0”. If B = “1”, the slot is busy; if B = “0”, the slot 
is free. Stations can only write into free slots of the correct type. 


2.2 Channel occupancy 


Let r in b/s be the code rate of a voice source, R in b/s be the 
network transmission rate, and H be the number of overhead bits. H 
includes the access field, the source/destination addresses, control, 
and Cyclic Redundancy Check (CRC). Then, the slot length in bits is 


given by 
F =H +rT bits/slot. (1) 


The number of bits transmitted by the system in T seconds is RT, 
and the number of bits in an information packet is rT. The voice cycle 
length, L, in slots, is then 





fh 
L= : ‘ slots. (2) 


Let C = 1 be the normalized channel capacity, i.e., the channel can 
be utilized to transmit at most one packet per slot. Since a slot is 
longer than the time necessary to transmit one information packet, 
the maximum channel occupany, Umax, is 


rT 
nix yy. A lots. 
U; Hatt slots (3) 
The capacity required to coordinate packet transmission, C,, is 
s = C at Oras (4) 
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2.3 Code designations of priority scheduling 


A voice cycle repeats every L slots. All admitted voice sessions are 
guaranteed one slot per voice cycle. It is expected that the number of 
slots taken by voice calls during any voice cycle will be less than L 
slots so that there will be slots available for data transmissions on a 
per-cycle basis. The concept of a voice cycle is depicted in Fig. 2, with 
a single slot expanded to emphasize the characteristics of the access 
field. The access field of each slot has three subfields, namely, S, E, 
and B. As indicated in Fig. 1a, the first (or head) station issues slots 
in the forward path, and the last (or end) station echoes acknowledg- 
ment in the reverse path. The S subfield contains a code word that 
designates the traffic type for that slot. The B subfield contains a 
token bit to indicate whether the slot is full or empty. The code word 
in the S subfield is set by the head station when the slot is issued. The 
B bit is to be set by the transmitting station. The end station observes 
whether a slot, after traversing the entire channel, is full or empty. If 
the slot is empty, then every station with the designated type of traffic 
has transmitted the permitted number of packets. The status of an 
empty slot is indicated (by the end station) in the E subfield. The 
occurrence of an empty slot signals the end of one type of service and 
the start of a new service type. The dead time associated with the 
circulation of an empty slot, as indicated to the head station by the 
marked E subfield, has been referred to earlier as the sojourn time. 

The S subfield, which takes on one of the code words H, G, V, D, 
and C, determines the mode of the slot. To ensure “fairness”, admission 


START OF START OF 
/TH CYCLE (j+ 1)TH CYCLE 
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CYCLE CYCLE 


Fig. 2—Service cycle. 
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of new voice is carried out in two stages. Upon arrival, a new call 
enters the wait state. When the new voice in the wait state sees an H 
in the S subfield, it transits to the hold state and ultimately transmits 
when a visiting slot is free. The station model for new voice admission 
is depicted in Fig. 3. The code words H, G, V, D, and C have the 
following connotation: 

H is a type designation that allows new voice in the wait state to 

transit to the hold state and to start competing for a slot 
G allows new voice in the hold state to start competing for a slot 
V allows old “active” voice to start competing for a slot 


NEW VOICE 


MOVE TO 
WAIT 






TRANSMIT 
PACKET 


BECOMES OLD VOICE 





Fig. 3—New voice admission algorithm. 
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D allows “active” data stations to transit from the wait state to the 
hold state and to start competing for a slot 

C allows active data stations in the hold state to start competing for 
a slot. 


A data subcycle differs from a voice subcycle. Since a data subcycle 
may be interrupted by a voice subcycle, a data subcycle may span 
more than one service cycle. The access control code, D, starts a new 
data subcycle after all data stations in the wait state have been served 
so as to maintain fairness: amongst the data stations. There are 
basically three types of services, namely, new voice calls, old voice 
calls, and data. Both H and G admit new voice calls, and both D and 
C allow data stations to transmit. However, each H (D) slot is issued 
at most once per voice (data) subcycle. That is, an H slot is followed 
by a string of G slots and a D slot is followed by a string of C slots. 
Fig. 4 shows the flowchart for the issuance of slots by the head station 
using feedback information, denoted by an asterisk, that arrives from 
the end station, after one channel propagation delay. The quantity Q,, 
determined in Section IV, is the number of new voice calls to be 
admitted in the jth admission cycle. 


Ill. EVOLUTION OF STATE EQUATION 


Consider just the voice subcycle. A station can be either in a voice 
session or idle. A voice session can be in a talk spurt or silent. Figure 
5 shows the state diagram of a voice source. 

Figure 5 depicts the state transitions at the start of the jth service 
cycle [or at the completion of the (j — 1)th cycle]. The symbols used 
have the following notations: 


nj-1 denotes the number of active voice stations in the (j — 1)th 
cycle 

dj;-1 denotes the number of stations moving to the silent state after 
the completion of the (j — 1)th cycle 

D;-. denotes the number of stations terminating their voice ses- 
sions at the completion of the (j — 1)th cycle 

b; | denotes the number of stations moving from the silent to the 
active state during the (j — 1)th cycle 

B; denotes the number of stations with new calls that are admit- 
ted in the jth voice cycle 

G; denotes the number of new calls arriving in the jth voice cycle. 


Then, the number of calls not accepted in the jth voice cycle is G; 
— B;. The state equation at the beginning of the jth voice cycle is 
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Fig. 4—Flow chart for slot issuance by head station with feedback from end station... 
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Fig. 5—State transition diagram for voice at start of jth service cycle, assuming 
deaths occur while in active state. 


Nj = nj-1 — Aj — Dj-1 + b; + B;. (5) 


Under the access control mechanism described in Section 2.1, nj-; is 
the only known quantity in the above equation. (It is possible to obtain 
information about D;-; and dj-; by the incorporation of a more 
complex access control mechanism.) We will describe, in the next 
section, an admission policy for new calls under the Fasnet structure 
described in Section II and the state development in eq. (5). 


IV. NEW VOICE ADMISSION CONTROL 


It is desirable to maintain an equilibrium operating state such that 
the mean number of active voice sessions (in talk spurt) will be 7 = 
Ni, where N, is the target number of active voice sessions. Let N be 
the total number of voice sessions to be supported and p be the fraction 
of voice sessions that are active. Then, 7 = pN. In this paper, we let p 
= 0.4, ie., 40 percent of the time a (nonidle) voice source is in talk 
spurt. 

Suppose we call each run of identically designated slots, i.e., H, G, 
V, D, and C, a sequence denoted by S,, x € {H, G, V, D, C}. The 
termination of a sequence is detected by the last station whenever it 
sees the first empty slot. Before a new sequence of slots can be 
initiated, a sojourn time, [,.;, elapses. I;.; equals one slot time plus one 
round trip propagation delay. (The round trip propagation delay equals 
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Fig. 6—Concept of sojourn time in jth cycle. 


2. [5] where 6, measured in slots, is the propagation delay in one 
direction, if the starts occur at the same instant in each channel, and 
the symbol [-] denotes the integer greater than or equal to its 
argument.) 

At the start of a new service cycle, the first station issues a sequence, 
Se, of 2; H- or G-designated slots (see Fig. 4). The number actually 
admitted is 


(6) 


p= G; if Gj <Q; 
Q; if G; = Q;. 


The format of a voice cycle is shown in Fig. 6. A sequence Sy will 
follow Sg. If G; = Q;, there will be no sojourn time separating Sg and 
Sy. Otherwise, there will be a sojourn time <I,,; slots separating Sg 
and Sy. Also note from Fig. 6 that there is at least one sojourn time 
between an Sy and an Sc sequence within a service cycle, where Sc 
denotes a sequence of D or C starts. 

Admission of new voice calls is controlled by the design parameter 
N,. Since a cycle time is quite short, it is only necessary to admit new 
voice calls once every P cycles. Let A; be the accumulated number of 
calls that left the system since the last admission cycle. Then A; may 
be estimated by 
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A z nj T 

a ogee pT’ 7) 
where T is the cycle length in seconds, T}, is the holding time of a 
voice session, and p is the probability that a voice session is in talk 
spurt. The value n,/p is an estimate of the total number of voice 
sessions in the system in the jth cycle. On the assumption that, at 
steady state, the net movement from silent to talk spurt and vice versa 
is zero, the maximum number of active voice in the jth admission 
cycle can be approximated as 


nj = Nj-p + B; = A;. (8) 


The objective is to drive the system to operate at a point such that 
the average number of active voice approximately equals the design 
parameter N;,. We propose to compute a running average ni; by the 
following recursive formula: 


nj = yy + (1 — y)nj, (9) 


where y, in the range 0 < y <1, is a constant factor that determines 
the relative weighting of the last running average and the current 
number of active voice. The constant y also fades the memory of the 
recursive process and hence affects the size of the moving average 
window. ‘To have a relatively large window, y should only be slightly 
less than unity, e.g., y = 0.98. Thus, during start-up, 7; will be small. 
By the same token, the number admitted is expected to be considerably 
larger than the number that leaves the system. Then, during start-up, 
the number of new voice already admitted should play an important 
role in controlling the admission process. When the system attains an 
equilibrium condition, the running average will approximately equal 
N, and should assume primary control of the admission process. That 
is, at equilibrium, the number admitted should balance the number 
that left the system. 

Let c; be the parameter that accounts for the difference between the 
number admitted and the number that left. The number admitted in 
the past admission cycles should only be remembered until the system 
attains steady state, ie., until 7; attains a value close to N,. It is 
proposed to compute the control variable c; by a recursive equation of 
the form 


cj = nej-1 + B; — Aj, (10) 


where 0 < 7 < 1 is chosen to fade the memory at an appropriate rate. 
n should be related to the time constant with which the system attains 
steady state from start-up. In fact, the combination of (9) and (10) 
forms a robust iterative algorithm so that an optimum choice of 7 is 


318 TECHNICAL JOURNAL, FEBRUARY 1984 


not critically important. The strategy is to admit new voice calls 
whenever ij + c; < N;. The number to be admitted in the jth cycle is 
governed by the equation 


Q; _ max(0, N igs nj — c;). (11) 


We have run numerous tests of the admission control algorithm 
described in this section. In all cases, the system quickly locks onto a 
small neighbourhood of the target operating point, N;. It is noted that 
admission control takes effect only when the total number of active 
voice exceeds N,. If the total number of voice calls is fewer than N, so 
that the total number of active voice in the system is fewer than Nz, 
no admission control would be needed, and all calls would be admitted 
into the system. 


V. QUEUEING ANALYSIS OF DATA PACKETS 


Although the topology and the media access protocol of Fasnet?® 
differ from those of Welnet,® the operational features of the two 
systems are similar. That is, each of these two systems has a deter- 
ministic media access protocol, which allows for priority scheduling 
and a bounded packet delivery time. With appropriate modifications, 
the periodic switching model and the mean value approximation 
method used in Ref. 8 to analyze the data buffer queueing behaviour 
are applicable here. The queueing analysis described below assumes 
that the procedure for admitting new voice calls is able to maintain 
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Fig. 7—(a) Single-server model. (b) Periodic switching model. 
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the voice service operating at the target value Nj, i.e., the mean number 
of active voice in the system approximately equals N,. 

The periodic switching model shown in Fig. 7 opens every T seconds 
and closes t; seconds after the opening of the switch at the start of the 
jth service cycle. Here, t; is the time needed to serve n; active voice 
plus the total (i.e., including also those sojourns occurring in the data 
subcycles) sojourn time incurred in the jth service cycle. We have 
hitherto used n; to represent the number of active voice calls. Since 
each active voice occupies one slot per cycle, in what follows, n; will 
also denote the number of slots occupied by voice in the jth cycle. 

Let m; denote the number of sojourns in the jth cycle. We then 
have 


t; = at Riel , in seconds, (12) 
where T, L, and [,.; are as previously defined. The time available for 
servicing data packets in the jth cycle is then T — t; seconds or (L — 
nj — Mjlzo3) slots. Since the transmission of each packet requires one 
slot time (L — n; — mjJ,.;), data packets can be transmitted in the jth 
cycle. We define 


_L-(nj=k) — mj; Le 
rs, 


to be the service rate available to data in the jth cycle. The mean 
service rate is then 


packets/s (13) 


a= Lo =a UES) packets/s. (14) 


Let \q be the aggregate data arrival rate in packet/s. The data utili- 
zation, pa, is defined by 
da 


pa = — 
Ma 


hal 

=e (15) 
For stable operation, it is required that pg< 1. For a given total number 
of voice calls in the system, N, there is a small probability (caused by 
the statistical fluctuation in the number of active voice calls in the jth 
cycle) that, in a certain jth cycle, \; > [(Z — (n; = k) — mjl,.;)|/T and 
the system is said to experience temporary overload, where ),; is the 
number of arrivals in the jth cycle. Under the fluid approximation 
model to be discussed below, we approximate ; by \a, the mean arrival 
rate, and define the data utilization in the jth cycle by 
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Fig. 8—Fluid approximation model for the jth cycle. 
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Depending on the value of k, p, can exceed unity. 

The mean value, or fluid, approximation method” is depicted in 
Fig. 8, where N,(v) is the critical number of active voice at which the 
service available to data in the fluid approximation model equals the 
number of data packets that have arrived in the same cycle. The jth 
cycle is on the verge of overload when the combination of voice, 
sojourn time, and data is such that the queue at the end is the same 
as at the beginning of the jth cycle, i.e., Qj11 = Q;. At this operating 
point, the number of active voice is N,(v). Ifn; > N.(v), there is a net 
growth in the queue and Qj, > Q;. The queueing process temporarily 
becomes unstable and the system experiences overload. Since the 
number of voice in talk spurt obeys a binomial distribution and since 
the mean time between transitions from silent to talk spurt, and vice 
versa, is long compared to a cycle time (e.g., 1.50 seconds compared to 
10 ms), once the system enters an overload state, there is every 
likelihood that it will stay in the overload state for many consecutive 
service cycles. 

To analyze the queueing behaviour, we separate the service cycles 
into two sets, a stable set in which the queueing process is stationary 
and an overload set in which the queueing process is nonstationary. 
In addition, during the interval when the system is servicing voice 
only, new data arrivals will have to be scheduled for transmission. Let 


(16) 
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Q, and Q, be the sets of stable and overload states, respectively. The 
average queue in the system is given by the sum of the contributions 
from the stable and overload states and the temporary queueing due 
to scheduling: 


Q = Qs) + Q(0) + Qecr, (17a) 
where 
Q(s) = E[Q|k € 2.]-Pr(k € 2) (17b) 
is the contribution from the stable states, 
Q(0) = E[Q|k € .]-Pr(k € 24) (17c) 


is the contribution from the overload state, and Q.n is mean queue 
due to scheduling of the data packets during a service cycle. 

An approximate queueing analysis method, based on a decomposi- 
tion of stable and overload operational modes, was developed in Ref. 
8 for the queueing analysis of integrated voice/data services on Welnet. 
Functionally, the media access protocols of Fasnet and Welnet differ 
only in that Fasnet may incur a sojourn time whenever a change in 
the type of traffic to be served takes place. With the incorporation of 
the sojourn time into the analysis model, the method developed in 
Ref. 8 can be used to analyze the data queueing behaviour on Fasnet. 

When the number of active data stations is small (at light loads), 
there will be many data subcycles within one service cycle and consec- 
utive data subcycles will be separated by a sojourn time. While the 
existence of sojourn time represents a loss of channel capacity, not 
every occurrence of a sojourn affects the data queueing process. The 
ultimate limit occurs when every data station has packets to be 
transmitted (as is likely to happen in an overload state). We define 
the mean number of sojourn times per service cycle, which represents 
a loss in channel capacity, by 


L- n(o) ed T+ To; 
M:- Prax 


where the 1 in the right-hand side of (18a) represents the sojourn 
between the end of the voice subcycle and the beginning of the data 
subcycle, 7(o) is the mean number of active voice in the overload state, 
M is the number of data stations, and Prax = 1 is the maximum 
number of packets that an active data station is permitted to send per 
data subcycle. A value of Pax > 1 is beneficial for small M, so that 
the mean number of sojourn times in the overload state is a minimum. 
For M > L — N,(v), a value of Pmax = 1 suffices. The expression for 
m can be rearranged to yield 


m= + 1, (18a) 
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_ L—nlo) + M-Praax 
= —___—___—_. 18b 
_ MEP as 4, ( i 
It can be shown that the mean number of voice while in the overload 
state, 7i(o), is given by Ref. 8: 


(n; = k)-Pr(n; = k) 
k=N,(v)+1 


N 
y% = Pr(n; = k) 


k=N,(v) +1 


ni(o) = (19) 


From Fig. 8, we deduce that 
hal = L — Nv) — Mg. (20) 


That is, when the number of active voice n; = N,(v), the number of 
packets arrived in the jth cycle equals the number of packets served 
in the same cycle. The above equation can be rearranged to yield 


N(v) = L — aT — rilegj- 


Substituting \aT' by pa(L — ni — mI,.;) from (20) and rearranging, we 
have 


N,(v) = L(1 — pa) + pa: — M-Tgoj(1 — pa). (21) 


The parameters m and N,(v) are related by a set of transcendental 
equations. m and N,(v) can be computed by the following algorithm: 


1. Substitute N.(v) for n(o) in (18b) and then combine the result 
with (21) to obtain 


[Li ==. pa) a7 npa|(M- Pax + Tsoj) 
_ T5o5(1 = pa)(L + M: Prnax) 
T503(1 _ pa)(L — M:- Prax) 
. Calculate N,(v). 
. Using this value of N.(v), calculate n(o) by eq. (19). 


. Calculate m using (18b). 
. Recalculate N,(v) using (21). 


N,(v) = 


ot fe cc bo 


The approximate analysis model developed in Ref. 8 is based on the 
fact that the number of active voice obeys a binomial distribution and 
that the system (for voice) is at equilibrium when the number of active 
voice equals the mean, n = pN, where p is the probability that a voice 
call is in talk spurt and N is the total number of voice calls in the 
system. When the number of active voice calls deviates from n, there 
is a tendency for the system to move towards the equilibrium point. 
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Define d, to be the net drift toward the equilibrium point. Then, at 
any operating point n; = k, the net drift, d,, is given by Ref. 8: 


d, = a-T.(nj =k) — B-T-(N — (nj = k))|kR EQ, slots/cycle, (22) 


where a7! and 6” are, respectively, the mean holding times of a talk 
spurt and a silent event and T, as previously defined, is the cycle 
length in seconds. With probability Pr[overload] = Pr[k € Q,] = 
Pr[n; > N.(v)], the system is in the overload state. Once it enters the 
overload state, the mean time that the system spends in the overload 
state is a function of the net drift d, given by (22). The mean and 
second moment of the duration that the system spends in the overload 
state is given by Ref. 8: 


. (N — (n; = k))Pr(n; = k) 


__ k=N,(v)+1 d, 
Orme N 
» (N- (nj = k))Pr(nj; = k) 
k=N,(v)+1 
x (N — (n; = k))Pr(n; = k) 


S (N= (n; = k))Pr(n; = k) 
k=N,(v)+1 
and 
5 (N — (n; = k))Pr(n; = k) 
= — ~:. 7: —_ = 2 
gp = MNwH ee (n, = k) — B-T-(N — (n; = k))] (24) 
>» (N - (n; = k))Pr(n; = k) 
k=NAv)+1 


The coefficient of variation of the duration in the overload state is 
then given by 


So 
[ro]? 


It turns out that C, < 1. This is a consequence of the fact that the 
increase in the number of active voice, which drives the data queueing 
process into an overload condition, is a slowly fluctuating process. 
Thus, the second moment of the overload duration has negligible effect 
on the queueing behaviour when the system is in the overload state. 
As observed in Ref. 8 and depicted in Fig. 9, a transient state always 
follows an overload state. The chain of service cycles depicted in Fig. 
9 will be referred to as an “overload/transient” chain. The contribution 
to the average queue in the system from an overload/transient chain 
is given by the queue in a typical overload/transient chain multiplied 





C= 1. (25) 
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Fig. 9—Chaining of overload/transient cycles. 


by the probability of the occurrence of an overload/transient chain. 
The mean duration of the overload portion of the chain depicted in 
Fig. 9 is t,. The queue growth in the overload state is given by the 
difference of the queue size at the ends of consecutive service cycles, 
ie., AQ; = gjui(v) — Q;(v). The average queue growth per cycle will be 
denoted by E[AQ]. The average queue growth in an overload state is 
then given by an accumulation of the growth per cycle over the 
overload duration, i.e., 
[ro] 
2 iB(AQ), 

where [x] denotes the integer part of x + 0.5, plus AgN.(v)P/2R, 
which is the average growth in the cycle when the system enters the 
transient state. The queue shrinkage will be slightly faster than the 
queue growth so that the mean duration in the transient state is 
nominally shorter than 7,. This is due to the fact that the number of 
active voice while the system is in the overload state is larger than 
that when the system is in the transient state. The average number of 
slots available for data per overload cycle is L — n(o) — m- Iq; that 
per transient cycle is L — n(t) — m-I,.;, where n(o) and n(t) are the 
mean number of active voice in the overload and transient states, 
respectively. ri(o) is given by (19), and n(t) <= N.(v). For convenience, 
we approximate 7i(t) by N.(v) and define 


L — A(o) — m-Lno3 


h= < 
L — N.(v) — m-Tooj 


(26) 


to be the fractional contribution from the transient portion to the 
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queue size in an overload/transient chain. Then, the number of packets 

in an average overload/transient chain is given by Ref. 8: 

E[Q|k € 24] oa 

[To | -( [ To | a2 1) 
2 


E[AQ] + 


Na: n(o) -P 
2R : 


=a+m{ 


where 


E[AQ] = NaT — {L — Ao) — m-Lee3 (28) 


is the average queue growth per cycle when the system is in an overload 
state. The average queue contribution from the overload/transient 
chain is then given by 


Q(o) = E[Q|k € Q,]-Prik € Oo]. (29) 


The stationary contribution can be written as 
N,(v)-1 


Q(s) = 2 E[Q|nj = k]Pr(n; = k), (30) 


where E[Q|n,; = k] is the mean number of packets in the system given 
the number of voice in the jth cycle is less than N,(v) and, for Poisson 
arrivals, is given by the M/G/1 formula: 


px(1 + C3) 
2(1 — pr) © 


px is the utilization defined by (16) and Cj is the squared coefficient 
of variation of service time. Since voice changes slowly compared to a 
cycle time, Cj < 1. 

The Pollaczek-Khinchin M/G/1 formula describes the queueing 
effect at the imbedded points of the service process which, in our case, 
correspond to the ends of service cycles. During the initial part of the 
service cycle when the system is serving voice stations, the data 
stations experience a temporary queueing. The average temporary 
queue per cycle corresponds to a delay during which a typical data 
packet is scheduled for channel access. Depending on the state of the 
system, the number served in a given cycle may be fewer than the 
number that has arrived during the same cycle. The queueing of those 
new arrivals that are not served in the same cycle is accounted for by 
the M/G/1 formula. To compute the queueing due to scheduling, it is 
only necessary to consider new arrivals when the system is in the 
stable mode. Consider the fluid approximation model depicted in Fig. 
10. The total number of packets that need to be queued temporarily 


E[Q|nj = k] = pet (31) 
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Fig. 10—Fluid approximation model for the evaluation of mean queue due to sched- 
uling. 


in the jth cycle is given by the area under the triangle; the average 
number of packets in the temporary queue is given by the area divided 
by the base. Thus, the average number of packets in the temporary 
queue due to scheduling during one cycle is given by 


=~ da NW 
Qsch = age > [(n; = k) + m+ Isoj]-Pr(n; = k). (32) 
2L  x=0 


The mean number of packets in the system is then given, as indicated 
in (17a), by the sum of Q(o), Q(s), and Q,ch, which are prescribed 
respectively by (29), (30), and (32). 


VI. NUMERICAL RESULTS 


Recognizing the difficulty with an exact analysis of the data 
queueing behaviour in an integrated voice and data service environ- 
ment, we have introduced an approximate analysis in Section V. 
Another approach to evaluating the system performance is by means 
of computer simulation or emulation. Emulation portrays all the steps 
involved in the actual system. The drawback in an emulation process 
is that it is extremely time-consuming and that it is not feasible to 
emulate a high-rate system. We compare the approximate analysis 
results with those obtained by means of computer emulation. Based 
on the closeness of the approximate analysis results to those obtained 
from emulation on a 10-Mb/s Fasnet, we conjecture that the perform- 
ance results obtained for a high-rate Fasnet (e.g., 100 Mb/s) using the 
approximate analysis method are reasonable portrayals of system 
performance. 


VOICE AND DATA ON FASNET 327 


6.1 Simulation 


The simulation consisted of an emulation at the packet level. We 
assume the system operates at 10 Mb/s. The length of the line was 
assumed to be 2.5 km, resulting in normalized propagation and signal 
processing delay for a single line of 0.2 slot. Data and voice stations 
were randomly positioned along the line. We assumed that 142 voice 
station pairs were in conversation with an additional 20 active data 
stations. The other parameters used in the simulation are r = 64 
kb/s, T = 10 ms, and P = 700 b/packet. 

An ongoing voice conversation is modeled as a two-state Markov 
process. It spends 40 percent of the time in the active state and 60 
percent in the silent state (activity factor = 0.4). The average talk 
spurt duration is 1.5 seconds and the average silent interval is 2.25 
seconds, both with an exponential length distribution. All data stations 
generated traffic at the same mean rate and each station transmits a 
maximum of three packets from its queue during one data subcycle 
(Pmax = 3). Data packets arrive at each station according to a Poisson 
distribution. 

Each simulation run was for 500,000 packets, corresponding to a 
duration of 35 seconds. The simulation started with voice conversa- 
tions assigned to talk and silent stations in the same proportion as 
the steady-state distribution. To ensure that the steady state was 
reached before data were recorded, statistics were not gathered until 
after the first 20,000 packets. 

Performance is measured in terms of the average number of packets 
in the system. In a given service cycle, there is temporary queueing 
when the system is serving the voice stations. The queue at the end of 
each cycle is smaller than it is during the cycle. The simulation results 
shown in Fig. 11 reflect the effects of temporary queueing due to 
scheduling during one service cycle. 

For each data point shown in Fig. 11, the results for between five 
and twenty independent runs were combined. At lower data rates, the 
variance is small and few runs were required to obtain accurate 
measures of the mean data queue. At higher data rates, the variance 
was significantly larger. The simulation did not guarantee that the 
average amount of voice traffic would be equal to the value used to 
calculate the theoretical results. This affects the measurement of the 
average queue size for data since the channel capacity available for 
data depends directly on the capacity taken by the voice stations. To 
improve the estimate of the mean queue for data, the average number 
of active voice conversations (which was measured for each run) was 
regressed against mean queue and, from this regression, the value at 
56.8 (0.4 X 142) active calls (the value used in the calculations) was 
estimated. An example of this is shown in Fig. 12, where the experi- 
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Fig. 11—Data queueing characteristics. 


mental points, regression line, and 95-percent confidence interval are 
shown. 


6.2 Calculation 


The same parameters used in the simulation were employed to 
calculate the average queue contributions from the stable and overload 
modes and channel access scheduling, using the equations developed 
in Section V. The calculated results are tabulated in Table I. For 
comparison purposes, the mean number of packets in the system is 
plotted in Fig. 11 for the cases with and without the mean scheduling 
queue per cycle. Inspection of Fig. 11 shows that the calculated results 
closely approximate those obtained from simulation. 

The data arrival rate in packets/slot is expressible as \g7/L. Then, 
by Little’s formula,’* the mean packet delay is given by 
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Fig. 12—Regression analysis of mean data queue vs. mean number of active voice 
with 95-percent confidence. 


Table |—Mean number of packets in the system as a function of the 
arrival rate Ag in packets/s* 


Pr[over- = a Z 2 

pa da N.(v)  m load] Qscn QS) Q(0) Q 
0.50 4091.80 97.72 1.68  2.6449e-12 5.19 1.01 — 6.21 
0.60 4891.12 89.43 1.82 1.6032e-08 6.86 1.38 — 8.24 
0.70 5691.21 81.19 1.95 1.4707e-05 8.78 1.96 — 10.73 
0.80 6483.12 73.01 2.08  2.3096e-03 11.02 3.25 0.23 14.51 
0.85 6877.15 68.94 2.15 2.3249e-02 12.01 4.71 4.31 21.03 
0.90 7269.90 64.88 2.21 9.4180e-02 12.16 6.86 30.56 49.58 
0.95 7895.67 60.83 2.28  2.6206e-01 10.54 8.64 228.68 247.86 
0.98 7895.67 58.41 2.32 3.8351e-01 8.02 6.30 1512.02 1526.34 


* R = 10 Mb/s, r = 64 kb/s, H = 60 bits, P = 700 b/packet, a“! = 1.50 seconds, 67 
= 2.25 seconds, L = 142 slots, N = L, M = 20, Prax = 3, Ino; = 2. 


d = QL/ dal’. 


The analysis and simulation results of the delay-throughput charac- 
teristics are shown in Fig. 13. These curves are an alternative pres- 
entation of the results shown in Fig. 11. 

To demonstrate the effects on data queueing when the number of 
voice calls is reduced, we compare in Fig. 14 the cases for N = 128, 
136, and 142. It is observed that at a relatively high data rate, the 
average queue decreases significantly when the number of voice calls 
is reduced by approximately 5 percent, ie., from 142 to 136. The 
improvement in queueing comes from the fact that a reduction in the 
number of voice calls is accompanied by a relatively large reduction in 
the probability that the system is in the overload state. 
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Fig. 13—Mean data packet delay characteristics. 


As the voice arrival rate r decreases or the network transmission 
rate increases, the number of slots in a cycle increases. For example, 
with R = 100 Mb/s, r = 64 kb/s, P = 700 b/packet, and T = 0.01 
second, the number of slots per cycle is approximately 1408. For r = 
32 kb/s and P = 380 b/packet (H = 60 bits in each case), the number 
of slots per cycle is approximately 2624. Letting the number of voice 
calls N = L, the average number of data packets in the system as a 
function of the traffic intensity pq is plotted in Fig. 15. It is noted that 
the mean contribution to the queue due to scheduling is the same for 
the two cases, since the packet length for the r = 64 kb/s case is 1.8421 
times longer than that of the r = 32 kb/s case. The longer cycle length 
(due to a shorter packet length) case exhibits better long-term 
queueing performance (excluding scheduling queueing). Consider for 
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Fig. 14—Effects on queueing with a decrease in the total number of voice stations. 


example, at pa = 0.95, the average number in the system (excluding 
scheduling queueing) is 18.21 packets for the r = 64 kb/s case and 
12.72 packets for r = 32 kb/s. Since the relative packet length ratio is 
700/380, the effective mean queue for the 64 kb/s case is 2.64 times 
that of the shorter packet case. Thus, the mean queue for the case 
where P = 700 bits is equivalent to 33.58 380-bit packets at pg = 0.95 
compared to 12.72 380-bit packets for the 32 kb/s case at the same 
utilization. 

For comparison purposes, the M/D/1 curve is also plotted in Figs. 
11 and 15. The M/D/1 result is based strictly on a stationary distri- 
bution, i.e., ignoring the statistical variations associated with the 
speech spurts. At light loads where the system always operates in a 
stable mode, the M/D/1 represents a reasonable lower bound. At 
heavy loads where the system experiences a reasonable percentage of 
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Fig. 15— Queueing characteristics for different packet lengths. 


overload, the M/D/1 curve serves no useful comparison (see Fig. 11). 
Of significance are the numerical results shown in Fig. 15, in which 
the perfect scheduling Fasnet results deviate from the M/D/1 curve 
by almost a constant amount up to 95-percent data utilization. This 
is a consequente of the fact that, at the higher speed (R = 100 Mb/s), 
overload does not become significant until the data utilization ap- 
proaches unity. The discrepancy between the M/D/1 and the perfect 
scheduling Fasnet curves can be attributed to the presence of sojourn 
time in the Fasnet access discipline. 


VII. DISCUSSIONS AND CONCLUSIONS 


In previous work, Limb and Flamm” described a scheme for accom- 
modating voice and data that differs from the present method in that 
calls returning from idle may get clipped. New calls are accepted until 
an unacceptable level of clipping is approached, whereupon they are 
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blocked. Unused voice capacity may be utilized by data, although a 
sojourn time is required to switch between the voice subcycle and the 
data subcycle unless the voice subcycle is full. 

A direct comparison between the two schemes is difficult because of 
the different way in which the voice traffic is treated—the new calls 
admission policy and the clipping of talk spurts in one case and the 
absence of clipping in the other. However, a couple of observations 
can be made: Both systems have comparable performance since they 
both utilize the available slots efficiently when traffic is heavy; clipping 
may be undesirable in certain instances, for example, where a number 
of Fasnet links are connected in tandem. A speech spurt could be 
clipped more than once, leading to a noticeable degradation. In this 
instance, the present method would be preferable. 

The global access control that is enabled by the endowment of each 
slot with a type designation renders Fasnet particularly suitable for 
integrated services. We have described the ramifications for integrated 
voice and data services on Fasnet and a policy for admitting new voice 
into the system. Once admitted, a voice station is guaranteed service 
when needed. Voice stations are served on a periodic basis, so that 
each voice source sees almost synchronous service, i.e., the service is 
synchronous within the statistical variations dictated by the binomial 
distribution of the silent and talk spurt intervals in a voice process. 
The mean interpacket delay can be expected to be negligibly small. 
Note that the service discipline described will not clip any voice 
packets from a talk spurt. 

Since data traffic can be buffered, it is not necessary to maintain 
continuous service for data. During intervals when the system is 
servicing voice stations, data stations are temporarily locked out and 
have to be buffered. Depending on the data input load relative to the 
service time available for data during one service cycle, the temporary 
queues may be cleared entirely at the end of the service cycle, or some 
residual packets will carry over from one cycle to the next. If the latter 
phenomenon occurs over many consecutive cycles, the system will 
sustain a net growth in its queues. Under this condition the system 
experiences overloading, and the data utilization during these cycles 
will exceed unity for a short period and the data queues will experience 
significant growth. 

The conventional approach to analyzing the queueing behaviour is 
the generating function analysis method. However, the size of the 
state space in a practical network renders the generating function 
method impractical. An approximate analysis method, based on a 
merge of a periodic service single-server model and a fluid approxi- 
mation method, has been developed in Ref. 8 to analyze the queueing 
behaviour for integrated voice and data services on Welnet, a local 
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area network in which the service discipline is functionally similar to 
that of Fasnet. With a suitable modification to account for the sojourn 
time that accompanies each switch in the type of service by the Fasnet 
access control, the approximate analysis method introduced in Ref. 8 
was used to analyze the data buffer queueing behaviour for integrated 
voice and data services on Fasnet. 

System performance may also be evaluated by means of computer 
simulation. Simulation is time-consuming and it is impractical to 
simulate high transmission rate systems. In this study, we simulated 
a 10-Mb/s Fasnet and measured performance in terms of the mean 
queue in the system as a function of the data input load. The analytic 
results closely approximate those obtained from simulation. It is 
conjectured that the analysis method yields reasonably accurate re- 
sults. Analytical results for a 100-Mb/s Fasnet have been presented. 
An interesting observation is that the larger bandwidth has the effect 
of reducing the system overload probability. 

In this paper we used an approximate model to analyze the network 
performance. We verified the reasonableness of the approximate 
model. An alternative approach is to approximate the correlated voice 
source by a memoryless process, e.g., a birth-and-death process. The 
validity of such an approximation needs to be shown: An analysis of 
the backlog in the buffers for a birth-and-death process is discussed 
in Ref. 14. 


REFERENCES 


1. K. Kummerle, “Multiplex Performance of Integrated Line and Packet-Switched 
fae Proc. 2nd Int. Conf. Comp. Commun., Stockholm, Sweden, August 1974, 
pp. 517-23. 

. M. J. Fischer and T. C. Harris, “A Model for Evaluating the Performance of an 
Integrated Circuit-and Packet-Switched Multiplex Structure,” IEEE Trans. Com- 
mun., COM-24 (February 1976), pp. 195-202. 

3. C. J. Weinstein, M. K. Malpass, and M. J. Fischer, “Data Traffic Performance of 
an Integrated Circuit-and Packet-Switched Multiplex Structure,” IEEE Trans. 
Commun., COM-28 (June 1980), pp. 873-7. 

. B. Maglaris and M. Schwartz, “Performance Evaluation of a Variable Frame 
Multiplexer for Integrated Switched Networks,” IKEE Trans. Commun., COM- 
29 (June 1981), pp. 800-7. 

. I. Gitman, W-N. Hsieh, and B. J. Occhiogross, “Analysis and Design of Hybrid 
Silene Networks,” IEEE Trans. Commun., COM-29 (September 1981), pp. 
1290-1300. 

. B. R. Iyer and J. B. Sinclair, “Analysis of Data Delays in an Integrated Data-Voice 
Link,” Technical Report 8204, Electrical Engineering Department, Rice Univer- 
sity, Houston, Texas, April 1982. 

. A. Leon-Garcia, R. H. Kwong, and G. F. Williams, “Performance Evaluation 
Methods for an Integrated Voice/Data Link,” IEEE Trans. Commun., COM-30 
(August 1982), pp. 1848-58. 

8. J. W. Mark, “Priority Scheduling for Integrated Voice/Data Services in Local Area 
Networks,” CCNG Report E-111, Computer Communication Networks Group, 
University of Waterloo, Waterloo, Ont., Canada, March 1983. 

9. J. O. Limb and C. Flores, “Description of Fasnet—A Unidirectional Local Area 
Communications Network,” B.S.T.J., 61, No. 7, Part 1 (September 1982), pp. 
1413-40. 


bo 


AN 


on 


jor) 


~] 


VOICE AND DATA ON FASNET 335 


10. J. O. Limb and L. E. Flamm, “A Distributed Local Area Network Packet Protocol 
for Combined Voice and Data Transmission,” GLOBECOM ’82 Conf. Rec., 
November 29-December 2, 1982, Miami, Florida, pp. B2.5.1-6. 

11. C. J. Weinstein, “Fractional Speech Loss and Talker Activity Model for TASI 
pe eens Speech,” JEEE Trans. Commun., COM-26 No. 8 (August 1978), 
pp. 1253-7. 

12. G. F. Newell, Applications of Queueing Theory, London: Chapman and Hall Ltd., 
1972, Chapter 2. 

13. J. D. C. Little, “A Proof of the Queueing Formula L = \W,” Oper. Res., 9 (May 
1961), pp. 383-7. 

14. S. Halfin, “The Backlog of Data in Buffers with Variable Input and Output Rates,” 
Proc. of the Second International Symposium on Performance of Computer- 
Communications Systems, Zurich, March 21-23, 1984. 


AUTHORS 


John O. Limb, B.E.E., 1963, Ph.D. (Electrical Engineering), 1967, University 
of Western Australia, Nedlands; Research Laboratory, Australian Post Office, 
1966-1967; AT&T Bell Laboratories, 1967-1983. Present affiliation Bell Com- 
munications Research, Inc. At AT&T Bell Laboratories, Mr. Limb worked for 
a number of years on the coding of color and monochrome picture signals to 
reduce channel capacity requirements, and has published widely in this area. 
He is now manager of the Communication Research Department in Bell 
Communications Research, Inc. His current research interest is wideband, 
local area networks. Fellow, Institute of Electrical and Electronics Engineers; 
member, Association for Computing Machinery, Association for Research in 
Vision and Opthalmology, Optical Society of America. 


Jon W. Mark, B.A. Sc., 1962, University of Toronto; M.Eng., 1968, Ph.D., 
1970, McMaster University; Canadian Westinghouse, 1962-1968; IBM, 
Thomas Watson Research Center, 1976-1977; AT&T Bell Laboratories, 1982- 
1983; University of Waterloo, 1970—. At the IBM Research Center Mr. Mark 
was a visiting scientist. At AT&T Bell Laboratories, he was a consultant to 
the Computer Technology Research Laboratory. He worked for a number of 
years in the areas of radar and sonar signal processing, equalization, data 
compression, and spread spectrum communication. At present, Mr. Mark is a 
Professor of Electrical Engineering at the University of Waterloo. His current 
research interests are network modeling and performance evaluation, satellite 
networks, and local area networks. Senior Member, IEEE. 


336 TECHNICAL JOURNAL, FEBRUARY 1984 


AT&T BELL LABORATORIES TECHNICAL JOURNAL 
Vol. 63, No. 2, February 1984 
Printed in U.S.A. 


Experimental lon-Implanted Bubble Memory 
Device With 16-um7? Cell 


By T. J. NELSON,* D. J. MUEHLNER,* V. V. S. RANA,* 
B. J. ROMAN,* and G. P. VELLA-COLEIRO* 


(Manuscript received April 1, 1983) 


This paper describes a 4-~m-period (16 ym?/bit) magnetic bubble device 
based on ion-implanted propagation patterns. Single-layer LaSmLuBiGa-IG 
films were implanted with Ne* and He* ions using densified photoresist 
masking patterns. Coarse (three eighths of a period) minimum features were 
used in the Permalloy and Al-Cu levels. A 21-Oe overlap bias range was 
obtained for generation, propagation, transfer, and detection at 55-Oe drive. 
Generation and detection were demonstrated to have a 39-Oe bias range at 
44-Oe drive, and a 39-Oe bias range for minor-loop propagation was obtained 
at =33-Oe drive. 


I. INTRODUCTION 


Using proposed scaling laws for Permalloy}-based magnetic bubble 
devices, Kyrder’ has concluded that 6 X 10° b/cm?, 1-um bubble devices 
should be achievable. Fontanta et al.” have studied devices using 4.75 
x 5.25-um? cells and 0.75-um gaps fabricated with 10X reduction, 
direct step on wafer projection printing. Quasistatic gate operation 
was achieved at 70-Oe drive field. Recently, Yanase et al.? demon- 
strated propagation on 4-um “club foot” circuits with 1-um gaps at 


* AT&T Bell Laboratories. 
t Permalloy is a trademark of the Allegheny Ludlum Steel Corporation. 
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and abstract, but no other portions, of this paper may be copied or distributed royalty 
free by computer-based and other information-service systems without further permis- 
sion. Permission to reproduce or republish any other portion of this paper must be 
obtained from the Editor. 


337 


drive fields above 45 Oe. It appears, therefore, that Permalloy-based 
bubble devices are already near their anticipated density limit. 

One-micrometer bubble devices with 30-ym? cell area have been 
operated‘ at 150 kHz, using Ion-Implanted Propagation Patterns? 
(I2P2s) instead of Permalloy. And Komenou et al.® have reported 
promising results using 16-ym? I2P2 cells. Because of their extreme 
simplicity (a sine-wave-shaped track will propagate, for example), 
I2P2 devices can be resolved with three-eighths of a period minimum 
features, which, together with a 0.75-um optical lithography limit, 
imply a possible density of 2.5(10)" b/cm?. 

This paper describes magnetic bubble material and test-circuit 
design and characterizes generation, propagation, transfer, and detec- 
tion in 4-~m-period devices based on ion-implanted propagation pat- 
terns (see Fig. 1).* With one exception, the minimum features in the 
Permalloy and Al-Cu levels were coarsened to three-eighths of a period. 
We have continued to use single-layer garnet films, and the layers 
processed on the chip were unchanged, except that the Permalloy 
thickness was reduced from 500 AU to 375 AU. 


Il. MATERIAL FOR 4-y~m-PERIOD BUBBLE DEVICES 


Films of the nominal composition (LaSmLu)3(FeGa)50,2, with and 
without Bi, were used in the present study. Propagation at 4-~m period 
in films of (LaSmLuCa)3 (FeGe);0i2 composition was previously re- 
ported.” However, this composition required a high magnetic moment 
(4rM, = 784G) to produce 1.1-um-sized bubbles. The saturation 
magnetization (47M,) is given by 4%M, = V32xAQ/Z’, where Z is the 
material length parameter, Q is the quality factor, and A is the 
exchange constant. Hence, a similar bubble size (same 7, and Q) can 
be obtained at a lower 47M, by going to a composition having a lower 
exchange constant. Since A scales with the Curie temperature,’ the 
required 47M, can be lowered by going to a composition having a lower 
Curie temperature. In the present work we have used Ga substitution, 
instead of CaGe, to reduce the 47M, and, possibly, the minimum drive 
field. The desired anisotropy was obtained from the SmLu site selec- 
tively. The bubble mobility was improved by making use of the 
Lanthanum effect;”* i.e., the Sm content was kept low and the uniaxial 
anisotropy was increased by increasing the Lu content and using La 
to match the film lattice parameter to that of the substrate, La having 
no effect on the uniaxial anisotropy. 


* The contents of this paper were presented by T. J. Nelson in the Fourth Interna- 
tional Conference on Magnetic Bubbles held September 24-27, 1980, in Tokyo. 
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Fig. 1—Minor-loop ion-implanted propagation patterns and orientation with respect 
to wafer, preferred positions, and rotating field conventions. 


Films were grown from PbO—B,0; flux by the Liquid Phase Epitaxy 
(LPE) technique on 50-mm diameter (111) GdsGa;O0,2 substrates by 
horizontal dipping with a unidirectional rotation of 100 rpm. The 
growth temperatures ranged from 872° to 880°C. Table I lists the 
properties of typical films. The film thickness and demagnetized stripe 
width intentionally exceed the 0.9 um desired for 4-um period to allow 
for the thickness of the implant. H, is the collapse field. The domain 
wall mobility (u) was calculated from the resonance line width using 
standard Ferromagnetic Resonance (FMR) technique. Initially, films 
were grown without bismuth (type A). The bubble visibility in these 
films was poor because of the low Faraday rotation (6). The addition 
of Bi was made (film type B) to produce an enhanced (and reversed) 
Faraday rotation to improve the bubble visibility. 

The anisotropy of the implanted films was determined from micro- 
wave resonance using the method described in Ref. 9. Measurements 
were made at 11.85 GHz on 50-mm diameter wafers after removing 
the film on the back side by etching in phosphoric/sulfuric acid. Figure 
2 shows the derivative microwave absorption vs. field of a film (type 
A). This film was implanted with 80/Ne/1E14, 270/Ne/2E14, and 130/ 
He/6E15 followed by a 300°C anneal. (The implantations are referred 
to in a shorthand form representing energy in keV/ion species/dose 
in ions/cm?.) The anisotropy fields of both the film and its implant 
can be read from Fig. 2. For this film, 

770 


Q= 1+ ay = 2.16. 


A sufficiently large change in anisotropy of about 6(10)* ergs/cm?® was 
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Table |—Properties of (LaSmLu)3(FeGa)5O 2 films 
Property 


OF 

(at 

» _-5000A) 
th sw 2 H, 4c7M, Hx (cm/ (degrees/ 

Type (um) (um) (um) (Oe) (G) (Oe) s/Oe) cm) 


A Without Bi 1.39 108 011 378 664 1485 485 1400 
B With Bi 1.22 1.07 O11 365 655 1650 400 $-—4500 


caused by the implantation. Figure 3 shows the visually obtained 
minor-loop bias margins vs. drive field on another film (type B). This 
film was implanted with 80/Ne/1E14 (no pattern) and 270/Ne/2E14 
+ 130/He/4E15 patterned implants. This sample was subjected to full 
processing, which included vacuum anneal at about 350°C for 1 hour 
prior to the Permalloy evaporation. The demagnetized strip width in 
the implanted region after processing was 0.85 yum, sufficiently close 
to the desired value. 


Ill. PROCESSING OF 4-ym-PERIOD BUBBLE DEVICES 


This study employed contact printing using chrome masters on 
glass. The test devices were originally encoded with 0.5-um address 
steps at 8-um period. For this study they were shrunk to 0.25-um steps 
at 4-um period using an Electron Beam Exposure System (EBES)?° 
option. Scanning Electron Microscope (SEM) photographs of the 
propagate patterns are shown in Fig. 4a as printed in AZ1350J pho- 
toresist in the ordinary way and Fig. 4b as transferred to densified 
AZ1350J by trilevel'’ method. The latter lithographic technique per- 
mits the transfer of a high-resolution image, which is more easily 
achieved in a thin resist layer, into a thick polymerized layer that is 
dimensionally stable under extreme implantation conditions. Plasma- 
deposited SiO» separates the thin and thick organic layers and, after 
plasma etching through the primary resist pattern, serves as a reactive 
ion-etch mask to complete the image transfer process. Because of the 
highly anisotropic etching, trilevel patterns have vertical edge profiles 
and were used exclusively to guarantee edge integrity for the implanted 
patterns in this study. 

The processing after implantation consisted of 2000 AU of plasma- 
deposited SiO, (P-SiO2), 375 AU of Permalloy, and an additional 2000 
AU of P-SiO, and 5000 AU of Al-Cu. Via holes were plasma etched 
into the second SiO, layer to connect the Permalloy sensors to their 
Al-Cu pads. The devices received no additional passivation. The 
Permalloy and Al-Cu layers were patterned with contact printing and 
ion-beam milling. As will be discussed below, the minimum features 
in the metal layers, except in the generator, were coded to be 1.5 um. 
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Fig. 2—FMR absorption (derivative) as a function of field of a (LaSmLu)3(GaFe);O1. 
film implanted with 80/Ne/1E14, 270/Ne/2E14, and 130/He/6E15. 
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Fig. 3—Minor-loop propagation margins vs. drive field at 40°C using 80/Ne/1E14 
(unpatterned) + 270/Ne/2E14 + 130/He/4E15. Open circles refer to two bubbles in 
adjacent positions. 


IV. LAYOUT AND OPERATION OF 4-y-PERIOD BUBBLE DEVICE 


Figure 5 shows the components and layout used in this study. The 
generator consists of a simple Al-Cu hairpin, much like that used in 
Permalloy devices. Between the generator and detector, a line of 
clover-leaf-shaped features was interpolated. These features were sep- 
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(b) 


Fig. 4—-SEM photographs of 4-um-period resist patterns (a) with contact-printed 
AZ1350J resist and (b) using trilevel resist. 


arated by 1-um merge gaps.’” The clover-leaf features have threefold 
symmetry, except for small EBES address effects, and so provide 
strong-side propagation all around. The transfer-enable was designed 
to frustrate propagation of bubbles by the merge gap and send them 
through the gaps and around the cover-leaf features. In this way, a 
block of bubbles coming from the generator may, when properly 
positioned, be brought to the transfer ports for insertion into the 
minor loops. The trapping-transfer gates’* operate by confining bub- 
bles within the lower slot of the conductor and between the two 
implant boundaries. Bubbles transfer from one path to the other after 
being trapped for approximately 180° of field rotation. When a group 
of bubbles is transferred out from the minor loops to the tops of the 
clover-leaf features, each bubble propagates around through the merge 
gap to the path between the generator and detector. Once the block is 
thus assembled, the bubbles propagate towards the detector with the 
same bit order as originally generated. The detector consists of an Al- 
Cu conductor hairpin to stretch the bubbles and a Permalloy stripe to 
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Fig. 5—Test-circuit layout used in design study at 4-um period. Minimum feature 
was 1 ym in the implant level and 1.5 um in the Permalloy and Al-Cu levels. 


sense them. The bubbles are collapsed by a second current pulse to 
restore the detector for the next bubble. 

The minimum coded feature used in the implant level was 1 ym, 
which was the distance between implanted regions through the nar- 
rowest place on the minor loops or the width of the implanted region 
in the gaps between clover-leaf features. Except for the generator, the 
minimum feature in the metal levels was coded to be 1.5 wm, which 
was the line width of the Permalloy sensor and the gap width in the 
Al-Cu hairpins. The Al-Cu line width in the transfer gates was made 
2.5 um, which is the maximum permitted by the 8+yum interloop spacing 
and the 1.5-um gap width. 

The experimental 4-~m-period devices (film type B) were operated 
on a microscope-based test set under microprocessor control. The 
temperature of the chips was unregulated, but as the duty factor for 
the drive coils was small and the bias magnet was water cooled, it was 
less than 25°C. All testing was carried out at 50 kHz with sine-wave 
currents in the drive coils. Figure 5 shows positions 1 through 3 where 
bubbles were stopped during the tests. A typical test consisted of (1) 
generation and propagation to 1; (2) propagation with transfer-enable 
and transfer-in to 2; (3) propagation around the minor loop back to 2; 
(4) propagation and transfer-out to 3; and, finally, (5) propagation and 
detection. The test was automatically repeated at 3-Oe bias intervals. 
The data pattern generated was compared with that actually detected, 
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and the exclusive OR of the two was printed for each test. The bias 
ranging could be restricted to a single component of the test, in which 
case the other components were performed at a fixed bias. All current 
pulsers were triggered at 0° in the appropriate cycles by the computer, 
and delays, widths, and amplitudes were manually adjusted on the 
pulsers. 


V. DEVICE COMPONENTS AT 4-um PERIOD 
5.1 Generation 


The generator design, shown in more detail in Fig. 6, was included 
in a serial test circuit also present on the chip. Ordinary ungapped 
propagation patterns were employed in the serial test circuit, which 
was used to obtain the generator and detector data given in Figs. 7, 9, 
and 10. A reduction in upper bias margin with increased pulse width 
was observed. At 100-ns pulse width, the full bias range was obtained 
with generator current varied almost by a factor of 2, from 110 to 200 
mA. 


5.2 Detection 


Figure 8 shows the detector design employed. The Al-Cu hairpin 
was pulsed to remove the bubble from the strong-side track and 
simultaneously stretch it along the Permalloy at a drive-field angle of 
90°. At this time the magnetization in the Permalloy switches from 
up to down. A decreased resistance is observed when a stripped-out 
bubble is present, indicating increased cross-magnetization. Figure 9 
shows the stretch-current amplitude dependence of the bias margins 
for the approximately 50-um-long detector in the serial test circuit. 
Good bias margins were obtained between 50 mA and 90 mA. Figure 





level. 
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Fig. 7—Bias range at 55-Oe drive and 25°C for generation at 0° as a function of 
current. 





Fig. 8—Destructive Read Out detector design using 1.5-~m Permalloy line width and 
2.5-um Al-Cu line width at 4-um period. 


10 shows the bias-range dependence on collapse current amplitude. 
The collapse pulse was applied typically quite late, which stroboscopic 
observation (at 8-um period) suggests gives the strip time to contract 
to a bubble before the collapse pulse is applied. Therefore, the collapse 
pulse is thought to collapse only a bubble, not a strip. Figure 11 shows 
the bias range vs. drive field for generation, strong-side propagation 
and detection. The bias range has opened up at 44-Oe drive, though 
operation was possible down to 33 Oe with some loss of low bias 
margin. The data presented in Figs. 7, 9, and 10 were taken at 55-Oe 
drive to be compatible with the transfer results given below. 

Figure 12 shows the differential voltage between detector and 
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Fig. 9—Bias range vs. detector stretch current at 55-Oe drive field with pulse active 
between 72° and 126°. 
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Fig. 10—Bias range vs. collapse pulse current at 55-Oe drive field with pulse timing 
252° to 288°. 


dummy, amplified by 10. This waveform was obtained at 55-Oe drive 
with the circuit for which transfer and minor-loop margins are given 
below. Reliable detection with this longer (100 wm compared to 50 wm 
in the serial test circuit) detector required earlier phasing for the 
stretch pulse, about 54° instead of 72°. In this particular circuit a 
higher collapse current was also needed. The 500 detector gave about 
6 mV at 2.5 mA, corresponding to a figure of merit of 6V/V = 0.5 
percent. 


5.3 Transfer-enable 


The transfer-enable and trapping-transfer ports are drawn to larger 
scale in Fig. 13. The major line is composed of the bottoms of the 
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Fig. 11—Bias range for generation, strong-side propagation, and detection at 25°C 
vs. drive field. 





Fig. 12—Differential detector output with bubble and no bubble traces superimposed. 
Gain is 10 for the signal and 1 mV/mA for the stretch and collapse pulses. Positive 
peak of the sine wave is at 0°. 
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Fig. 13—Transfer-enable and trapping-transfer designs used to communicate be- 
tween minor loops and gapped major line. 


clover-leaf features. While bubbles propagated past the gaps at lower 
drive, at least 50 Oe was needed to propagate bubbles away from the 
clover-leaf cusps. The data reported below were obtained with 55-Oe 
sine-wave drive. Transfer-enable bias range vs. pulse current is shown 
in Fig. 14. Four bubbles, aligned with four ports, were made to pass 
through the merge gaps by a pulse active near a rotating field phase 
of 0°. The bubbles started from the cusps on the major line, passed 
through the gaps at the end of the first cycle, and propagated three 
steps around the clover leaf. A 33-Oe bias range was obtained for 25 
through 45 mA pulses active from —72° to 36°. Some variation has 
been seen with other chips of the same wafer. In other tests, for 
example, the low bias limit increased with increasing current instead 
of being independent of current as in Fig. 14. 


5.4 Trapping transfer 


Figure 13 also shows trapping-transfer gates linking the tops of the 
clover-leaf features with the bottoms of the minor loops. For the data 
presented in Figs. 15 through 20, a single bubble was loaded into each 
of three adjacent 73-step minor loops. Note that the trapping-transfer 
conductor slot opens away from the minor loops. This means that the 
pulse increases the bias on neighboring positions in the minor loops. 
In this work, the effect of the pulse current on neighboring positions 
was investigated by pulsing one cycle before and one cycle after the 
bubbles propagated by the gates within the minor loops. The open 
circles in Figs. 15 through 20 represent, with these extra pulses, both 
transfer-in and transfer-out. The slot orientation was chosen to avoid 
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ae ata range vs. transfer-enable current at 55-Oe drive with pulse active from 
—72° to 36°. 
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Fig. 15—Bias range vs. transfer-in leading-edge phase at 55 Oe. Trailing edge was 
held constant at 144° and the current was 35 mA. 


the bad conductor crossings that developed after passivation in 8-~m 
trapping-transfer gates. In the present design, the bad crossing would 
be expected if the conductor slot extended past the implant boundary 
of the bottom of the minor loop. In the 8-um-period gates, the crossings 
were still good after passivation when the slot stopped short, as in Fig. 
13, and left continuous metal over the bubble path. 
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Figure 15 shows bias range vs. leading-edge phase of the transfer-in 
pulse. The pulse could be started as late as +36°; thereafter, the upper 
bias margin dropped abruptly. The bias range was fairly constant for 
earlier phases. Transfer-in trailing-edge phase was varied to produce 
Fig. 16. Figure 17 shows the bias range vs. transfer-in current, with 
the lowest current for reliable operation at 30 mA. The high bias effect 
of the pulse on neighboring positions in the minor loop was most 
noticeable in this plot. The overall margins were limited by transfer- 
in at low bias, which improved as the current was lowered. Figure 18 
shows the bias range vs. transfer-out leading-edge phase. Figure 19 
shows that a higher upper bias limit exists for trailing-edge phases 
near 360°. Such phasing, however, loses high bias margin when the 
effect on bubbles in adjacent cycles is considered. This high bias 
failure mode is also clearly evident in Fig. 20, which shows the bias 
range vs. transfer-out current. The overall margins were in fact limited 
at high bias by transfer-out. 


VI. OVERALL BIAS RANGE AND CONCLUSIONS 


Figure 21 shows the bias ranges for individual device components 
at 55-Oe drive. These results were obtained with the chip also used to 
produce the data shown in Figs. 14 through 20. The enable, transfer, 
and detector designs have been discussed above. The generator in the 
major-minor test circuit had 1-um gap and 1.75-um coded line widths. 
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Fig. 16—Bias range vs. transfer-in trailing-edge phase at 55 Oe with leading edge at 
0° and the current at 35 mA. 
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Fig. 17—Bias range vs. transfer-in current at 55 Oe with pulse active from 0° to 144°. 
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Fig. 18—Bias range vs. transfer-out leading-edge phase at 55 Oe. Pulse current was 
45 mA and lasted until 306°. 


In other work we have found 1.5-um gaps to be satisfactory for 
generation at 4-um period. A generator pulse of 105 mA, lasting 80 ns, 
worked well on the chip of Figs. 14 through 21. In other nominally 
identical chips, higher currents, over 200 mA, were required for gen- 
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Fig. 19—Bias range vs. transfer-out trailing-edge phase at 55 Oe. The leading-edge 
phase was 180° and the current was 45 mA. 
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Fig. 20—Bias range vs. transfer-out current at 55 Oe with pulse active from 180° to 
306°. 


eration. The wide minor-loop margins indicate that low end failure 
due to stripping or hopping between loops is not the limiting problem 
at 55 Oe. However, Fig. 3 shows that the low bias limit for minor-loop 
propagation would increase if the drive could be decreased. Transfer- 
out limited the overall range at high bias. Figure 20 shows that the 
high end loss in transfer-out limited is due to interference with bubbles 
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Fig. 21—Bias range for test components at 55 Oe. The generator pulse was 105 mA 
and 80 ns wide and occurred at 0°. Transfer-enable lasted from —72° to +36° at 35 mA. 
Transfer-in was 35 mA from 0° to 144°. Transfer-out was 45 mA from 180° to 306°. 


in adjacent minor-loop positions. The next limit at high bias was due 
to detection. The low bias limit was caused by transfer-in. Since the 
low bias limit increases with increasing transfer-in current, some sort 
of stripping is indicated. The low drive limit, 50 Oe, was caused by 
bubbles sticking in the cusps of the clover-leaf-shaped major path 
features. 

In summary, an experimental 4-um-period functional major-minor 
bubble device has been operated at 55-Oe drive with 21-Oe bias range. 
Generation and detection were demonstrated to have a good bias range 
down to 44-Oe drive, and good minor-loop propagation was obtained 
as low as 33 Oe. 
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Theoretical Design of Single-Layer Antireflection 
Coatings on Laser Facets 


By G. EISENSTEIN* 
(Manuscript received November 19, 1982) 


Single-layer antireflection coating films are used to transform semiconduc- 
tor injection lasers into different kinds of active devices such as superlumi- 
nescent diodes or optical amplifiers. In this paper, optimum film parameters 
(thickness h and index of refraction n) are established for a wide range of 
InGaAsP lasers emitting at 1.3 um. Optimum film parameters are different 
for the TE and TM polarizations. The minimum theoretical power reflectivity 
is higher for the TM polarization than the minumum reflectivity for the TE 
polarization by a factor of about three. Both are very low, on the order of 10°. 
Tolerances in film parameters for a power reflectivity R < 10° (which is 
acceptable for most practical applications) are calculated for a typical laser 
having a spot size a = 0.5 ym. The tolerances are Ah = +50A and An/n = +3 
percent for the TE polarization and Ah = +43A and An/n = +2.4 percent for 
the TM polarization. Processing of high-quality antireflection coating films 
on InGaAsP devices is possible according to these tolerances by using sputtered 
SisN,, which allows a very slow deposition rate (on the order of 75A/min) and 
the tailoring of the film index by adjusting the nitrogen pressure in the plasma. 


I. INTRODUCTION 


Single-layer antireflection (AR) coating films can be very useful in 
transforming semiconductor injection lasers into different kinds of 
active devices. For example, a laser whose emitting facet reflectivities 
(one or both) are reduced to zero is transformed into a superlumines- 
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cent diode’ or an amplifier-modulator,”? respectively. Such matching 
layers are analyzed in this paper. The purpose of the analysis is to 

1. Establish the optimum parameters (thickness and index of re- 
fraction) of a matching layer for a given laser 

2. Estimate the minimum theoretical reflection from a matching 
layer 

3. Calculate the allowed tolerances of the film parameters for an 
acceptable reflectivity 

4, Compare the reflectivities of TE and TM polarized light for 
possible explanation of the observed preferential polarization in prac- 
tical devices.' 


Il. THEORY 


The device shown in Fig. 1 is assumed to have the properties 
described in the following paragraphs. 

The laser diode’s index of refraction is n, = 3.52, which is that of 
the InGaAsP active region in which the field is largely confined. We 
ignore the step in the index due to the cladding material to simplify 
the calculations. Kaplan and Deimel have shown that in most cases 
this simplification has a very small effect on the resulting reflectivity.® 
In some cases, however, this index step cannot be ignored and one 
should use a weighted average of the indices or resort to the rigorous 
calculation presented in Ref. 3. 

The emitted beam is assumed to be Gaussian in both transverse 
directions. A Gaussian functional form is a good approximation to the 
emitted beam of practical devices’ and enables us to characterize the 
laser in terms of a measurable parameter, the spot size. The spot size 
measured is that of the emitted beam and is slightly different from 
the spot size of the field inside the device. This is because of the 
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Fig. 1—Geometry of the laser and AR film. 


358 TECHNICAL JOURNAL, FEBRUARY 1984 


angular-dependent reflection coefficient of the AR film. The difference 
between the two is negligible however. 

The spot size in the direction perpendicular to the junction (X 
direction) is assumed to be much smaller than the corresponding spot 
size in the Y direction. Therefore, it is reasonable to analyze a two- 
dimensional model in the XZ plane (see Fig. 2). This model is valid 
for most lasers with the possible exception of buried structure lasers, 
which emit a nearly circular beam. The electric field is assumed to be 
linearly polarized along the Y, or along the X directions (TE or TM 
polarization, respectively.) The problem at hand is approached in a 
manner similar to the problem of calculating the modal reflectivity of 
a laser.*° There, a narrow guided beam is impinged on a semiconduc- 
tor-air interface, while here, the same guided beam impinges on an 
AR matching layer. 

The technique used is that of the angular plane wave spectrum 
representation.® The calculation is an extension of two previous papers 
by R. H. Clarke,’* and details are given in the appendix. 


Ill. RESULTS 


The power reflectivity R of an optimum AR layer was calculated for 
a wide range of laser spot sizes. Detailed results for a 0.5-um spot size, 
typical of practical devices, are shown in Figs. 3, 4, and 5. (The spot 
size is defined as the electric field radius at its 1/e point.) 

The reflected power in the case of TE polarization is shown either 
in Fig. 3a as a function of film thickness assuming optimum film index 
Nz = 1.8346, or in Fig. 3b as a function of film index assuming optimum 
normalized thickness hn2/\ = 0.2606. Similar results are shown in 
Figs. 4a and b for the TM polarization. The optimum index and 
thickness are nz = 1.9361 and hn2/d\ = 0.2511, respectively. In these 
expressions XQ is the free-space wavelength = 1.3 um. 

The minimum reflectivity under the theoretically optimum condi- 
tions is very low (R < 10°°). This extremely low reflectivity is of little 
importance, however, since achieving it in practice would require 
tolerances too stringent for the film parameters. Assuming, for ex- 


Fig. 2—Definition of the coordinate system and the angle of incidence. 
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Fig. 3—Reflectivity: (a) as a function of film thickness (with optimum index), and 
(b) as a function of index (with optimum thickness) in the TE case for spot size a = 0.5 
pum and ) = 1.3 wm. 
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Fig. 4—Reflectivity as a function of: (a) thickness (with optimum index), and (b) as 
a function of index (with optimum thickness) in the TM case for spot size a = 0.5 um 
and \ = 1.3 um. 
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Fig. 5—Tolerances in thickness and index of refraction for R = 10-3 for TE and TM 
polarizations. 


ample, that a reflectivity R < 107° is acceptable, we deduce from Fig. 
5 the tolerances in the film parameters to achieve this reflectivity. 
They are assuming nz = 1.8346 and hn2/d\ = 0.2606, Ahn2/dX = +0.00744 
(Ah = +50A) and An2/n2 = +3 percent for the TE polarization, and 
assuming ng = 1.9861 and hno/X = 0.2511, Ahn2/dA = 0.2511, 
Ahn2/d = +0.0064 (Ah = +43A) and An2/nz = +2.4 percent for the 
TM polarization. The tolerances in film thickness and index can be 
traded for each other according to the contours described in Fig. 5. 

Similar calculations were performed for a wide range of laser spot 
sizes. Optimum film parameters as a function of laser spot size are 
shown in Figs. 6a and b for the TE and TM polarizations, respectively. 
In both cases spot sizes a > 3.5 um (for \ = 1.3 ym) approach the 
condition of a plane wave, i.e., hn2/A = 0.25 and nz = Vn,. Small spot 
sizes, however, deviate significantly from these asymptotic conditions. 
Figure 6 should be used as a design curve for optimum AR films for a 
given laser. 


IV. DISCUSSION 


The reflectivity of AR matching layers on the emitting facet of an 
injection laser was calculated for a wide range of laser spot sizes. The 
calculation assumes that the optical field is Gaussian and that the 
spot size in the direction perpendicular to the junction is much smaller 
than the corresponding spot size in the direction parallel to the 
junction. This allows for a two-dimensional model that is valid for 
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Fig. 6-—Optimum film parameters as a function of spot size for: (a) TE polarization 
and (b) TM polarization. 


most lasers. For lasers with an essentially circular spot the results 
presented here can be used only qualitatively. 

Optimum film index and thickness were calculated, as well as the 
allowed tolerances, to ensure a reflectivity R = 107°, which is an 
acceptable reflectivity for most practical applications. Optimum film 
parameters for TE and TM polarizations were found to differ slightly. 
The minimum reflectivity for TM polarization is slightly higher than 
for TE polarization. Moreover, the tolerances for the TE case are 
higher. Therefore, we can assume that in practical devices, the reflec- 
tivity for the TM component of the optical field is higher than for the 
TE component. However, a practical superluminescent diode has most 
of its power polarized TE (more than 75 percent).’ We have to 
conclude, therefore, that this observed preferential polarization is not 
due to the reflectivity but rather due to differences in internal gain 
and losses. 

The reflectivity of practical AR films on superluminescent diodes 
was found to be on the order of 107~* to 10™? (see Refs. 9 and 10). This 
reflectivity is quite realistic according to the tolerances calculated. 

Finally, the tolerances calculated tell us that the use of sputtered 
SisN, on InGaAsP devices can result in high-quality AR coatings. 
Indeed, the slow deposition rate (on the order of 75A/min) and the 
possibility of tailoring the film index by adjusting the nitrogen pressure 
in the plasma allow the processing of films with the required toler- 
ances. 
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APPENDIX 
Calculation of Modal Reflectivity 


We assume a Gaussian beam (with spot size a) to the left of the 
laser facet (Z = 0) 
E(X, 0) = exp(—X?/a’). (1) 


A time-variation exp(iwt) is assumed throughout. This field can be 
represented as an angular spectrum of plane waves incident at different 
angles © (see Fig. 2): 








F;(s) = < { E(x, 0) exp (i mm | dx, (2) 
where s = sin® and A) = free space wavelength = 1.3 um. Equation 2 
becomes 
Vran 
F;(s) = exp [> Py \* (3) 
Tan, 


The amplitude reflection coefficient r(s) for each plane wave inci- 
dent on the AR film (thickness h, index of refraction n2), followed by 
a semi-infinite region (n3 = 1) is" 


= Nig + re3exp(—278) 


= 1+ r2123eXp(—2i8) : (4) 
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where 
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and rio(s), re3(s) are the well-known Fresnel reflection coefficients at 
the first and second boundaries, respectively. (Note that the Fresnel 
coefficients for the TE or TM polarizations are different from one 
another.) 

The reflected wave F,(s) is easily found to be 


F,(s) = r(s)Fi(s). (6) 


The power reflection coefficient of the film R is found by calculating 
the portion of the reflected wave that is coupled back into the laser 
mode and the proper normalization, i.e., 
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where Q = cos@ in the TE case and Q = (cos® — sin?0)/cos® in the 
TM case.® The integration range is from s = —» to s = , The 
corresponding angle is real for |s| < 1 and becomes complex for 
|s| > 1; its value is 0 = +7/2 + iy. W varies from 0 to © as |s| varies 
from 1 to «©. The integrals in eq. (7) are solved numerically and the 
results are plotted in Figs. 3, 4, 5, and 6. 
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Coherent optical communications, integrated optics, and interferometric 
sensors require that single-mode fibers maintain a stable state of polarization. 
With a broadband source, wavelength averaging of the power in either orthog- 
onally polarized HE,, mode permits a simple measurement of the power 
transfer to the cross-polarized state (the h value). Polarization-holding quality 
was measured on several single-mode fibers with strain birefringence that had 
been induced by a new preform deformation method or a gas-phase etching 
technique. The fibers exhibited h values from 2.3 x 10~° m™ (average polari- 
zation-holding to -24 dB) to 8.9 x 10~* m™ with beat lengths (at \ = 632.8 
nm) from 1.7 to 10 mm when the power is wavelength averaged over the 800- 
to 1100-nm range. 


I. INTRODUCTION 


Coherent optical communications,’ integrated optics, and interfer- 
ometric sensors require that single-mode fibers maintain a stable state 
of polarization over long propagation lengths. Environmental pertur- 
bations are one cause of cross-polarization coupling of the two orthog- 
onally polarized HE, modes?“ resulting in an unstable state of polar- 
ization. The birefringent properties of polarization-maintaining fibers 
should be designed to be insensitive to environmental changes. 

An important characteristic of single-polarization fibers is the av- 
erage rate at which power is coupled into the cross-polarization state. 
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With a broadband source, wavelength averaging of the power in either 
mode can allow a simple measurement of the power transfer to the 
cross-polarized mode.® Following that approach, a statistically deter- 
mined polarization-holding parameter, h, and beat lengths (at \ = 
632.8 nm) were measured for several single-mode fibers. The strain 
birefringence was induced by a preform deformation method’ or a gas- 
phase etching’ technique. The power is wavelength averaged over the 
800- to 1100-nm spectral range. 


Il. RATE OF TRANSFER OF CROSS-POLARIZATION 


One parameter commonly used in describing birefringent fibers is 
the mode beat length defined by 


_ 247 _ dX 
~ 66 B’ 
where 66 is the difference between the propagation constants for the 
slow- and fast-polarization eigenmodes, ) is the wavelength of light, 
and B is the modal or normalized linear birefringence. The modal 
birefringence can be separated into a stress anisotropy component, B,, 
and a geometrical, or shape, component, Bg. 

Kaminow’ obtained an analytical expression for the ensemble av- 
erages of the relative powers in the two polarization modes. Rashleigh® 
et al. extended the random-coupling theory to incorporate a finite 
range of wave numbers on a single fiber since the random-coupling 
theory requires measurements of a large number of fibers. The rate of 
transfer of cross-polarization in a single-mode fiber is defined by 


£ = {P,J/{P} = 1/2[1 — exp(—2hL)], (2) 


where P, indicates the power in the cross-polarization mode, P = P, 
+ P, is the total power, P, is the power in the excited mode, and L is 
the fiber length. The inverse of the parameter h describes the char- 
acteristic distance for power transfer to the cross-polarization state. 
The ensemble of N statistically equivalent fibers is directly propor- 
tional to the bandwidth 6) of a broad bandwidth light source;? i.e., 


N = oLid ® L/ASZ (3) 


where A 7 represents one member of an ensemble of equivalent fibers 
and 


Ly (1) 





= 1 HB) 1 E — a8) 


c dk, clk dk, (4) 


is the polarization-mode dispersion, the group delay difference between 
orthogonally polarized modes. The relative standard deviation of the 
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rate of power transfer to the cross-polarization state is of the order 
NW. 


Ill. EXPERIMENTAL TECHNIQUE 


Figure 1 shows a schematic for the measurement of the rate of 
power transfer to the cross-polarization state. The source for these 
experiments is white light from a xenon arc lamp that is chopped and 
then filtered by Corning glass filters 5-56 and 2-58. The combined 
filters produce a bandwidth of 800 to 1100 nm. The filtered light is 
linearly polarized and focused into the single-mode test fiber. Clad- 
ding-mode strippers are used to ensure that the light being examined 
is guided along the fiber core. The polarized fiber output light is 
focused through a linear analyzer and detected by a silicon avalanche 
photodiode. The detected signal is input to a lock-in amplifier and 
referenced to a synchronizing signal from the chopper. The lock-in 
function-output voltage is digitized and sent to a computer where P, 
and P, are averaged in time. 

Stolen and Turner® showed that the best polarization-holding azi- 
muthal angle in a fiber can be determined by rotating the polarizer 
and analyzer until the detected power is minimum. At the minimum 
power, the measured output is in the cross-polarization state. To 
examine the power in the excited mode, the linear analyzer is rotated 
90 degrees. 


IV. RESULTS AND DISCUSSION 


The rate of power transfer to the cross-polarization state and beat 
lengths (at \ = 632.8 nm) were measured and are presented in Table 
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Fig. 1—Measurement system for the rate of power transfer to the cross-polarization 
state. 
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I. Table I also shows the average polarization holding in decibels 
defined as 10-logio (h-fiber length). The measured fibers were all 
single mode for wavelengths longer than 800 nm. 

Fibers 2 through 4 were fabricated with a cladding refractive index 
below that of the silica support tube (depressed index or W-type fiber 
design). The ellipticity of the stress-applying region (shown in Fig. 2) 
immediately surrounding the core gives rise to the birefringence ob- 
served. The stress resulting from this structure has been analyzed in 
detail by Stolen.® Fibers 2 through 4 contain an Al,03-SiO; core, B203- 
P.05-SiO2 stress region, and a fluorine-doped SiO, cladding. Two 
modes were present in these fibers at a wavelength of 0.63 um, and a 
high fundamental mode tunneling loss was observed in the 1.0-um 
wavelength region. 

Fiber 1 (shown in Fig. 3) was made by a gas-phase etching tech- 
nique.® It contained a GeO2-SiOz core, a fluorine-doped silica stress 
region, and a SiO, cladding. This fiber has a core diameter of 3.6 um, 
and was multimode at 0.63 wm. The long wavelength, leaky-mode loss 
edge was beyond 1.5 um. 

The h values for the fibers listed in Table I vary from 2.3 x 107° m™ 
to 8.9 x 10°* m7! with the beat lengths varying from 1.7 mm to 10 
mm. There was no correlation between the changes in the h values 
and the beat lengths. This agrees with other results in the litera- 


Table |—Average polarization-holding quality over 800- to 1100-nm 
spectral range 


Fiber Polarization Core 
Fiber Length h Lp Holding Diameter 
Number (m) (m7!) (mm) (dB) (um) 
1 96 6.2 x 107+ 10 -12.2 3.6 
2 50 8.9 x 107* 2.6 -13.5 4.3 X11 
3 450 1.5 x 1074 3 -11.7 5.8 
4 174 2.3 x 107° 4.3 -24 6 


-——-- SiOz SUPPORT TUBE 





———-F-Si02 OUTER CLADDING 


Fig. 2—Single-mode birefringent fiber cross section fabricated by a preform defor- 
mation. 


368 TECHNICAL JOURNAL, FEBRUARY 1984 






———Ge07-Si07 CORE 
——-Si0. SUPPORT TUBE 


—~-/-—— F-Si0. STRESS REGIONS 


Fig. 3—Single-mode birefringent fiber cross section fabricated by a gas-phase etching 
technique. 


ture.°-’ The fiber-core diameters are listed in Table J. Fiber 2 has 
shape-induced birefringence due to the elliptical core, but this effect 
is small compared to the stress-induced birefringence. 

The coherence length, A 4 for lightwaves in a fiber is 

Age (5) 
Oo 
where Ar = \2/céX is the coherence time of the source and , is the 
mean wavelength for the spectral width 6). Using the definition for o 
from eq. (4) and neglecting the second term (small for stress-induced 
birefringence), eq. (5) becomes 
No 
Af= Sy Lp. (6) 
As an example, fiber 1 has a beat wavelength, L, = 10 mm, which 
results in a coherence length, A 7 = 3.2 cm, for a source whose spectral 
width is 6\ = 300 nm centered about \, = 950 nm. For fiber 1 in Table 
I, this measurement represents an average of N ~ 3000 independent 
samples with a relative standard deviation of N~’? ~0.018. 

It has been shown that birefringent multimode fibers preserve the 
state of polarization.’° The observation was that in fibers supporting 
a moderate number of modes, all modes maintain linear polarization 
with the exception of one mode near cutoff. Multimode polarization 
holding measurements were made in the wavelength region from 500 
to 900 nm. The multimode results for fiber 1 indicate polarization 
holding to -27 dB, corresponding to an h value of 2.2 x 107° m™. 

A 1-km length of high-birefringence fiber fabricated by Hitachi 
Cable Ltd. was reported in the literature’! to have an h-value of 5.5 X 
10-° m™ for an average polarization holding to —-22 dB. These values 
are in the same range as those measured here. However, low-loss high- 
birefringence fibers have been reported to hold polarization to better 
than -30 dB.” 
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V. CONCLUSION 


In conclusion, the polarization quality of high-birefringence single- 
mode fibers has been evaluated using a broadband source. The bire- 
fringent fibers were fabricated by a preform deformation method and 
a gas-phase etching technique. A statistical determination of the rate 
of power transfer to the cross-polarized state was made by wavelength 
averaging. The h values measured in the single-mode regime of the 
fibers varied from 2.3 X 10-° m7 to 8.9 x 10° m“! with beat lengths 
from 1.7 mm to 10 mm. Polarization preservation was found in the 
multimode regime supporting a moderate number of modes, but it is 
not yet understood how mode mixing affects polarization performance. 


VI. ACKNOWLEDGMENTS 


We would like to thank K. L. Walker and R. Csencsits for providing 
the etched-preform-type fiber and acknowledge valuable discussions 
with L. G. Cohen and R. H. Stolen. 


REFERENCES 


1. Y. Yamamoto and T. Kimura, “Coherent Optical Fiber Transmission System,” 
IEEE J. Quant. Electr., QE-17 (June 1981), pp. 919-35. 

. I. P. Kaminow, “Polarization in Optical Fibers,” IEEE J. Quant. Electr., QE-17 

(January 1981), pp. 15-22. 

. T. Okoshi, “Single-Polarization Single-Mode Optical Fibers,” IEEE J. Quant. 
Electr., QE-17 (June 1981), pp. 879-84. 

D. N. Payne, A. J. Barlow, and J. J. Ramskov Hansen, “Development of Low- and 
High-Birefringence Optical Fibers,” IEEE J. Quant. Electr., QE-18 (April 1982), 
pp. 477-88. 

. S. C. Rashleigh et al., “Polarization Holding in Birefringent Single-Mode Fibers,” 

Optics Letters, 7 (January 1982), pp. 40-2. 

. R. H. Stolen, W. Pleibel, and J. R. Simpson, unpublished work. 

. K. L. Walker and R. Csencsits, unpublished work; R. D. Birch, D. N. Payne, and 
M. P. Vardham, “Fabrication of Polarization Maintaining Fibres Using Gas- 
Phase Etching,” Electron. Lett., 18 (November 1982), pp. 1036-8. 

8. R. H. Stolen and E. H. Turner, “Faraday Rotation in Highly Birefringent Optical 
Fibers,” Appl. Optics, 19 (March 1980), pp. 842-5. 
9. R. H. Stolen, “Calculation of Stress Birefringence in Fibers by a Infintesimal 
Element Method,” J. Lightwave Tech., LT-1, No. 2 (June 1983), pp. 297-301. 
10. R. E. Wagner, R. H. Stolen, and W. Pleibel, “Polarization Preservation in Multi- 
mode Fibres,” Electron. Lett., 17 (March 1981), pp. 177-8. 
11. S. C. Rashleigh and M. J. Marrone, “Polarization-Holding in a High-Birefringence 
Fibre,” Electron. Lett., 18 (April 1982), pp. 326-7. 
12. T. Katsuyama, H. Matsumura, and T. Suganuma, “Low Loss Single-Polarization 
Fibre,” Electron. Lett., 17 (June 1981), pp. 473-4. 


~~ win 


10 oOo 


AUTHORS 


Frederick M. Sears, S.B. and S.M. (Mechanical Engineering), 1975, Mas- 
sachusetts Institute of Technology; Ph.D. (Mechanical Engineering), 1980, 
University of California at Berkeley; Lawrence Livermore National Labora- 
tory, 1978-1980; AT&T Bell Laboratories, 1980—. At the Lawrence Livermore 
National Laboratory, Mr. Sears developed experimental and theoretical tech- 
niques for digital ultrasonic testing of geologic materials to examine their 
anelastic properties. At AT&T Bell Laboratories, he is currently pursuing 


370 TECHNICAL JOURNAL, FEBRUARY 1984 


research and development on optical-fiber characterization techniques as a 
Member of Technical Staff. Recent investigations include fiber interferometry, 
single-polarization fibers, polarization-maintaining fibers, polarization mode 
dispersion, and optical time-domain reflectometry. Member, Optical Society 
of America, the American Geophysical Union, and the National Technical 
Association. 


J. R. Simpson, B.S. (Glass Science), 1972, Alfred University, Alfred, NY; 
M.S. (Ceramic Science), 1978, Rutgers University, AT&T Bell Laboratories, 
1976—. At AT&T Bell Laboratories, Mr. Simpson’s research has included the 
characterization of optical fibers for fiber preform development and methods 
of increasing deposition rates in the Modified Chemical Vapor Deposition 
(MCVD) process. As a Member of Technical Staff, he has pursued the study 
of alternate dopants for the MCVD process and for making single-mode single- 
polarization optical fiber. 


OPTICAL FIBERS 371 


PAPERS BY BELL LABORATORIES AUTHORS 


COMPUTING/MATHEMATICS 


Dodd S. L., Shape-Preserving Spline Interpolation for Specifying Bivariate 
Functions on Grids. IEEE Comp G 3(6):70-79, 1983. 

Fishburn P. C., Packing Graphs With Odd and Even Trees. J Graph Th 7(3):369- 
383, 1983. 

McConkey J., The 2nd-Annual-Symposium on Small Computers in the Arts 
(Editorial). Comput Mus 7(3):25-30, 1983. 


ENGINEERING 


Alles D. S., Direct Electron-Beam Lithography (Editorial). Sol St Tech 
26(9):125, 1983. 

Beni G., Hackwood S., Hornak L. A., Jackel J. L., Dynamic Sensing for Robots— 
An Analysis and Implementation. Int J Rob R 2(2):51-61, 1983. 

Bank R. E., Rose D. J., Fichtner W., Numerical-Methods for Semiconductor- 
Device Simulation. IEEE Device 30(9):1031-1041, 1983. 

Chen C. Y., Alavi K., Cho A. Y., Garbinski P. A., Pang Y. M. New Gao.47 Ino.53As 
Sheet-Charge Field-Effect Transistor for Long-Wavelength Optoelectronic 
Integration. Electr Lett 19(19):791-792, 1983. 

Chu T.S., An Imaging Beam Waveguide Feed. IEEE Antenn 31(4):614-619, 1983. 
Coughran W. M., Grosse E., Rose D. J.,. CAZM—A Circuit Analyzer With Macro- 
modeling. IEEE Device 30(9):1207-1213, 1983. 

Cox D. C., Time-Division Adaptive Retransmission for Reducing Signal Im- 
pairments in Portable Radiotelephones. IEEE Veh T 32(3):230-238, 1983. 
Farrow R. C. A Model for Microcomputer Control of the Scanning Electron- 
Microscope Used in Basic Research. Scan Elec M1983(P2):479-488, 1983. 
Fichtner W., Rose D. J., Special Issue on Numerical-Simulation of VLSI Devices- 
Foreword (Editorial). IEEE Device 30(9):967, 1983. 

Fichtner W., Rose D. J., Bank R. E., Semiconductor-Device Simulation. IEEE 
Device 30(9):1018-1030, 1983. 

Goodman D. J., Sundberg C. E., Quality of Service and Bandwidth Efficiency of 
Cellular Mobile Radio With Variable Bit-Rate Speech Transmission. IEEE 
Veh T 32(3):211-217, 1983. 

Hackwood S., Beni G., Hornak L. A., Wolfe R., Nelson T. J.. A Torque-Sensitive 
Tactile Array for Robotics. Int J Rob R 2(2):46-50, 1983. 

Iannuzzi M., Bias Humidity Performance and Failure Mechanisms of Nonher- 
metic Aluminum SICs in an Environment Contaminated With Clz. IEEE 
Compon 6(2):191-201, 1983. 

Lamorte M. F., Abbott D. H., AlGaAs/GaAs Cascade Solar-Cell Computer Mod- 
eling Under High Solar Concentration. Solar Cells 9(4):311-326, 1983. 

Liou K. Y., Single-Longitudinal-Mode Operation of Injection-Laser Coupled to 
a Grinrod External Cavity. Electr Lett 19(19):750-751, 1983. 

Mayo J. S., Design Automation—Lessons of the Past, Challenges for the Fu- 
ture. IEEE Comp G 3(6):13-18, 1983. 

Penumalli B. R., A Comprehensive (2)-Dimensional VLSI Process Simulation 
Program, Biceps. IEEE Device 30(9):986-992, 1983. 

Petersen W. P., Fichtner W., Grosse E. H., Vectorized Monte-Carlo Calculation 
for the Transport of Ions in Amorphous Targets. IEEE Device 30(9):1011-1017, 
1983. 

Pini R. et al., Ultraviolet Four-Photon Mixing in a Multimode Silica Fiber 
Raman Amplifier. Opt Commun 47(3):226-229, 1983. 

Shashaty A. J., Nonlinear and Hysteretic Twisting Effects in Ocean Cable 
Laying. J Energ Res 105(3):341-345, 1983. 


373 


MANAGEMENT/ECONOMICS 


Luss H., An Inspection Policy Model for Production Facilities. Manag Sci. 
29(9):1102-1109, 1983. 


PHYSICAL SCIENCES 


Ashkin A., Gordon J. P., Stability of Radiation-Pressure Particle Traps—An 
Optical Earnshaw Theorem. Optics Lett 8(10):511-513, 1983. 

Aur S., Kofalt D., Waseda Y., Egami T., Wang R., Chen H. S., Teo B. K., Local- 
Structure of Amorphous MosoNiso Determined by Anomalous X- Ray-Scatter- 
ing Using Synchrotron Radiation. Sol St Comm 48(2):111-115, 1983. 

Auston D. H., Smith P. R., Generation and Detection of Millimeter Waves by 
Picosecond Photoconductivity. Appl Phys L 43(7):631-633, 1983. 

Bjorkholm J. E., Ballman A. A., Localized Wet-Chemical Etching of InP Induced 
by Laser-Heating. Appl Phys L 43(6):574-576, 1983. 

Cage M. E., Dziuba R. F., Field B. F., Williams E. R., Girvin S. M., Gossard A. C., Tsui 
D. C., Wagner R. J., Dissipation and Dynamic Nonlinear Behavior in the Quan- 
tum Hall Regime. Phys Rev L 51(15):1374-1377, 1983. 

Carroll. P. J., Patterson G. D., Cullerton S. A., Light-Scattering Photography of 
Polymers Near the Glass-Transition. J Pol Sc PP 21(10):1889-1896, 1983. _ 
Chraplyvy A. R., Stone J., Burrus C. A., Optical Gain Exceeding 35 dB at 1.56 
um Due to Stimulated Raman-Scattering by Molecular Dz in a Solid Silica 
Optical Fiber. Optics Lett 8(7):415-417, 1983. 

Feibelman P. J.. Hamann D. R., Electronic-Structure of Metal Overlayers on 
Rhodium. Phys Rev B 28(6):3092-3099, 1983. 

Goodby J. W., Relationships Between Miscibility Grouping and Smectic Meso- 
phase Structure. Molec Cryst 92(6):171-177, 1983. 

Johnson P. D., Smith N. V., Production of Circularly Polarized-Light From 
Synchrotron Radiation in the Vacuum Ultraviolet Region. Nucl Instru 214(2- 
3):505-508, 1983. 

Khanarian G., Tonelli A. E., Non-Linear Electrooptic and Dielectric Properties 
of Flexible Polymers (Review). ACS Symp S'1983(233):235-244, 1983. 

Lovinger A. J., Ferroelectric Transition in a Co-Polymer of Vinylidene Fluoride 
and Tetrafluoroethylene. Macromolec 16(9):1529-1534, 1983. 

Luther L. C., Kometani T. Y., Chemical-Analysis of Bi-Containing Magnetic- 
Bubble Garnet-Films. J Am Ceram 66(9):619-622, 1983. 

Ostermayer F. W., Kohl P. A., Burton R. H., Photo-Electrochemical Etching of 
Integral Lenses on InGaAsP/InP Light-Emitting-Diodes. Appl Phys L 
43(7):642-644, 1983. 

Ota Y., Silicon Molecular-Beam Epitaxy (Review). Thin Sol Fi 106(1-2):3-136, 
1983. 

Pang Y. M., Chen C. Y., Gardinski P. A., 1.5 GHz Operation of an Al,Ga;-xAs 
GaAs Modulation- Doped Photoconductive Detector. Electr Lett 19(18):716-717, 
1983. 

Patterson G. D., Carroll P. J.. Photon-Correlation Spectroscopy of Polymers in 
the Glassy State. J Pol Sc PP 21(10):1897-1902, 1983. 

Pulito V. L., Miller D. L., Sassa S., Yamane T., DNA Fragments in Friend- 
Erythroleukemia Cells Induced by Dimethyl-Sulfoxide. P Nas Biol 80(19): 
5912-5915, 1983. 

Rowe M. W., Gallagher P. K., Gyorgy E. M., Establishing the Absence of Any 
Influence by an External Magnetic-Field Upon the Intrinsic Rate of Reduction 
of Magnetite by Hydrogen. J Chem Phys 79(7):3534-3536, 1983. 

Schneider M. V., Cho A. Y., Kollberg E., Zipath H., Characteristics of Schottky 
Diodes with Microcluster Interface. Appl Phys L 43(6):558-560, 1983. 

Tell B., Bjorkohm J. E., Beebe E. D., Subnanosecond Pulsed Laser Annealing of 
Se- implanted InP. Appl Phys L 43(7): 655-657, 1983. 

Tsang W. T., Olsson N. A., Enhanced Frequency-Modulation in Cleaved-Cou- 


374 TECHNICAL JOURNAL, FEBRUARY 1984 


pled-Cavity Semiconductor-Lasers With Reduced Spurious Intensity Modula- 
tion. Appl Phys L 43(6):527-529, 1983. 

Valley J. F., Slusher R. E., Acoustic-Wave Calibration for Co2-Laser Scattering 
Experiments. Rev Sci Ins 54(9):1157-1162, 1983. 

Wertheim G. K., Disalvo F. J., Buchanan D. N., Site Inequivalence in 
Fe,,,Nb3-;Seio. Phys Rev B 28(6):3335-3338, 1983. 

White J. C., Henderson D., Anti-Stokes Raman Laser-Emission at 149-nm in 
Atomic Bromine. Optics Lett 8(10):520-522, 1983. 

Yamada M., Tien P. K., Martin R. J., Nahory R. E., Ballman A. A., Double Zinc 
Diffusion Fronts in InP—Theory and Experiment. Appl Phys L 43(6):594-596, 
1983. 

Yan M. F., Cannon R. M., Bowen H. K., Chowdhry U., Effect of Grain-Size 
Distribution on Sintered Density. Mater Sci E 60(3):275-281, 1983. 


Zucker J. E., Pinczuk A., Chemla D. S., Gossard A., Wiegmann W., Raman-Scattering 
Resonant With Quasi-Two-Dimensional Excitons in Semiconductor Quantum 
Wells. Phys Rev L 51(14):1293-1296, 1983. 


SOCIAL AND LIFE SCIENCES 


Frauenthal J. C., Swick K. E., Limit-Cycle Oscillations of the Human-Popula- 
tion. Demography 20(3):285-298, 1983. 

Kiger J. I., Glass A. L., The Facilitation of Lexical Decisions by a Prime Occur- 
ring After the Target. Mem Cognit 11(4):356-365, 1983. 


Ortony A., Turner T. J., Antos S. J., A Puzzle About Affect and Recognition 
Memory. J Exp Psy L 9(4):725-729, 1983. 


PAPERS BY BELL LABORATORIES AUTHORS = 375 


CONTENTS, MARCH 1984 


Signals Designed for Recovery After Clipping—III: Generalizations 
B. F. Logan, Jr. 


Click Modulation 
B. F. Logan, Jr. 


Design Diagrams for Depressed Cladding Single-Mode Fibers 
W. T. Anderson and P. F. Glodis 


Single GRIN-Lens Directional Couplers 
F. H. Levinson and S. W. Granlund 


The DTWP: An LPC-Based Dynamic Time-Warping Processor for 
Isolated Word Recognition | 
M. K. Brown, R. T. Thorkildsen, Y. H. Oh, and S. S. Ali 


On the Application of Embedded Training to Connected-Letter 
Recognition for Directory Listing Retrieval 
L. R. Rabiner, J. G. Wilpon, and S. G. Terrace 


An Improved Word-Detection Algorithm for Telephone-Quality 
Speech Incorporating Both Syntactic and Semantic Constraints 
J. G. Wilpon, L. R. Rabiner, and T. Martin 


Performance of Cross-Polarized M-ary QAM Signals Over 
Nondispersive Fading Channels 
M. Kavehrad 


377 











AT&T BELL LABORATORIES TECHNICAL JOURNAL is abstracted or indexed by 
Abstract Journal in Earthquake Engineering, Applied Mechanics Review, Applied 
Science & Technology Index, Chemical Abstracts, Computer Abstracts, Current 
Contents/Engineering, Technology & Applied Sciences, Current Index to Statistics, 
Current Papers in Electrical & Electronic Engineering, Current Papers on Computers & 
Control, Electronics & Communications Abstracts Journal, The Engineering Index, 
International Aerospace Abstacts, Journal of Current Laser Abstracts, Language and 
Language Behavior Abstracts, Mathematical Reviews, Science Abstracts (Series A, 
Physics Abstracts; Series B, Llectrical and Electronic Abstracts; and Series C, Computer 
& Control Abstracts), Science Citation Index, Sociological Abstracts, Social Welfare, 
Social Planning and Social Development, and Solid State Abstracts Journal. 
Reproductions of the Journal by years are available in microform from University 
Microfilms, 300 N. Zeeb Road, Ann Arbor, Michigan 48106. 


== Atel 


Bell Laboratories 


